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Lift on Inclined Bodies of Revolution in 


Hypersonic Flow 


G. GRIMMINGER,* E. P. WILLIAMS,?# anv G. B. W. YOUNG? 
The RAND Corporation 


SUMMARY 


The importance of body lift lies in the fact that at moderate 
angles of attack and high Mach Number it can constitute an 
appreciable part of the total lift of a winged missile. In this 
paper an attempt has been made to analyze body lift in hyper 
sonic flow by an approximate method and, together with a cor- 
relation of existing experimental data, to indicate the probable 
variation of body lift over a wide range of Mach Numbers extend- 
ing from low supersonic to hypersonic. The method of analysis 
of hypersonic flow over inclined bodies of revolution employed 
herein has been denoted as the hypersonic approximation. It is 
an improvement on the Newtonian corpuscular theory of aero 
dynamics, since it considers the centrifugal forces resulting from 
the curved paths of the air particles in addition to the impact 
(Newtonian) forces. 


SYMBOLS 


= cross-section area 
unit vector tangent to a circular parallel 
normal force coefficient = N/gorr,? 
local normal force coefficient = dCy/d(x/dg) 
(Pp — po)/go = pressure coefficient 
axial force coefficient = X/gorrp? 
local axial force coefficient = dCy/d(x/dz) 
= diameter 
length 
free-stream Mach Number 
rate of mass flow 
= normal force—that is, force normal to longitudinal axis 
of body 
= unit vector normal to surface, positive inward 
Static pressure 
(1/2) V? = dynamic pressure 
= radius of normal curvature of a stream line 
radius 
surface area 
Presented at the Aerodynamics Session, Eighteenth Annual 
Meeting, I.A.S., New York, January 23-26, 1950. Revised 
and received April, 1950. 
* Physical Scientist. Now with the Department of Defense, 
Washington, D.C., as Scientific Warfare Adviser. 
t Aerodynamics Engineer. 
t Research Engineer 


= unit vector tangent toa meridian 

= velocity 
free-stream velocity vector 
velocity of flow over the surface of the body 

that is, force in direction of longitudinal 
axis of body 

distance along longitudinal axis of body 


axial force 


angle of attack 

angular position of a point on the surface of the body 
upper limit of integration around a circular parallel 
angle between body stream line and a meridian 
thickness of body layer 

angle between Vi andn 

angle between ¢t and the longitudinal axis of the body 
semivertex angle of cone 

mass density 

density in body layer 


Subscripts 


= based on maximum body dimensions 
centrifugal force values 
converging or boat-tail section 
nose or diverging section 
refers to expansion flow 
impact values 
straight or cylindrical section 
based on free-stream conditions 


(I) INTRODUCTION 


i» THE FIELD OF GUIDED MISSILES, investigations of 
possible performance will inevitably lead to the con- 
sideration of flight at higher and higher Mach Num- 
bers. Although for some missiles the Mach Number 
My = 5 might be considered as the upper end of the 
speed range, for missiles of a different category the 
high-speed range of flight can conceivably extend to 
My = 20 or 25. As far as wing aerodynamics is con- 
cerned, for supersonic speeds up to My = 5 the wing 
lift may be obtained with satisfactory accuracy by 
means of the linearized supersonic wing theories, and at 
higher Mach Numbers satisfactory results are ob- 
tained on the basis of two-dimensional gas dynamics. 
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LIFT ON INCLINED 


body of revolution. The resulting forces on a cone and 
cylinder are then given. Centrifugal force effects re- 
duce the cylinder normal force resulting from the New- 
tonian analysis by approximately 10 per cent. Corre- 
sponding effects for slender cones and ogives are less 
jor the angle-of-attack range of general interest (a < 
94,); therefore, for practical purposes, the Newtonian 
analysis needs no modification for predicting the lift on 
the nose of a body of revolution at extremely high Mach 
Numbers. A qualitative discussion of the pressures on 
cone and cylinder areas situated in regions of expansion 
flow at hypersonic speeds is then presented. It serves 
as a guide for extending the results of the following 
correlation of experimental lift data through the hyper- 
sonic region to the hypersonic approximation values. 
A more detailed discussion and analysis of centrifugal 
force effects are given in the Appendix. 
(11) THe NEWTONIAN (IMPACT) AERODYNAMIC 
FoRCES ON A YAWED Bopy OF REVOLUTION 


In view of the general lack of exact gas dynamic re- 
sults for three-dimensional flow, it is important to real- 
ize that valuable results concerning the lift on a pointed 
body of revolution can be obtained from the relatively 
simple theory of Newtonian aerodynamics. At high 
supersonic Mach Numbers, particularly when the angle 
of attack is appreciable, the gas pressure forces on a 
body may be approximated in a simple manner on the 
basis of the concept of Newtonian flow.*~> In Newton- 
ian flow it is assumed that the gas stream maintains its 
speed and direction unchanged until it strikes the solid 
surface exposed to the flow, whereupon it loses the com- 
ponent of momentum normal to the surface and moves 
along the surface with the tangential component of 
momentum unchanged. Thus, in this concept the 
shock wave is assumed to lie on, or follow, the surface of 
the body. The Newtonian approximation does not 
specify the pressure on surfaces that do not “‘see’’ the 
flow—that is, surfaces on which gas dynamics would pre- 
dict expansion flow. For a flat plate inclined at an 
angle to the flow, the Newtonian pressure coefficient on 


the lower surface is 


2 sin? a (1) 


C, = (pb — po)/G 


where p) denotes free-stream pressure and go = poVo?/2 


is the free-stream dynamic pressure. 


The concept of Newtonian flow can also be ap- 
proached from the exact two-dimensional gas dynamical 
equations by letting M, —~ ©. Considering a flat plate 
inclined to the flow, there will be shock flow over the 
lower surface and expansion flow over the upper sur- 
As the Mach Number increases, the shock wave 
approaches closer and the lower surface 
(leading to increasing pressures) and the amount of ex- 


face. 
closer to 


pansion increases on the upper surface (leading to de- 
creasing pressures). In the limit when co, both 
the pressure and the pressure coefficient on the upper 
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TRANSVERSE SECTION 
(VIEWED FROM STERN) 
r= RADIUS VECTOR IN yz- 
PLANE TO A POINT ON 
THE SURFACE 





ELEMENT OF 
SURFACE dS 


V, = FREE-STREAM VELOCITY VECTOR 

n = UNIT VECTOR NORMAL TO A SURFACE ELEMENT 

t =UNIT VECTOR TANGENT TO SURFACE ELEMENT AND 
LYING IN A PLANE CONTAINING THE x-AXIS 

d, = MAXIMUM DIAMETER OF BODY 

£, = TOTAL LENGTH OF BODY 

#. = LENGTH OF DIVERGING PORTION OF BODY 

@, = SEMI-VERTEX ANGLE AT THE NOSE 

Fic. 3. Body of revolution inclined at an angle to Newtonian 

flow. 


surface become zero, and the pressure coefficient on the 
lower surface becomes® 


= (2) 


6a (y + 1) sin? a Cy 


where y is the ratio of the specific heats. Since the 
upper surface pressure coefficient is zero, Eq. (2) is also 
the expression for the normal force coefficient Cy. It is 
indicated in reference 4 that y ~ 1 as My ~ ©, which 
brings Eq. (2) into agreement with the Newtonian re- 
sult, Eq. (1). It is worth pointing out that for a given 
Mach Number and angle of attack the Newtonian theory 
gives considerably better results for a three-dimensional 
body than for a two-dimensional body. Thus, the 
Newtonian hypersonic approximation for a pointed body 
of revolution (such as a cone, for example) is surprisingly 
good and is much better than for a two-dimensional 


flat plate. 


(A) Arbitrary Body of Revolution 


To derive the Newtonian pressure forces on a body of 
revolution of general shape, consider the body shown in 
Fig. 3 for which the longitudinal axis of symmetry is the 
positive x-axis. The angle of attack a is the angle 
in the xz-plane between the free-stream velocity vector 
V, and the positive x-axis. The y-axis is perpendicular 
to the xz-plane, forming a right-handed system of co- 
ordinates. Consider a differential element of surface 
area dS at an arbitrary point O’ on the surface of the 
body, and let x’, y’, z’ denote a local right-handed sys- 
tem of coordinates at the point O’, such that x’, y’, 2’ are 
parallel respectively to x, y, z. Let r, 8 be the polar 
coordinates in the yz-plane of the point O’ on the sur- 
face of the body. Let n be a unit vector that is normal 
to the surface element dS and positive in the inward 
direction. Let ¢ be a unit vector that is tangent to the 
element dS, forming an angle @ with the positive x’-axis. 
Let 6 denote a unit vector that, together with n and 
t, forms a right-handed coordinate system, such that 
b = n Xt. For a body of revolution, which is the 
only case to be considered here, ¢t is tangent to a me- 
ridian and therefore lies in the plane formed by the line 
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y’ 
Fic. 4. Diagram of directions and components at a local element 


of surface area. 


O’x’ and Ox, and 6 is a tangent to a circular parallel 
and therefore lies in the y’z’-plane. The angular posi- 
tion of the point O’ (and the vector n) is given by the 
angle 8, which is measured positive counterclockwise 
from the positive y’-axis. The angle @ between O’x’ 
and ¢ is considered positive in the sense of a counter- 
clockwise rotation about 0. 

The relations between the n, t, b coordinate system 
and the x’, y’, z’ system (Fig. 4) are given in Table 1. 


TABLE | 
Direction Cosines 
x y’ Z 
sin 6sin B 


, 


t cos 6 sin 9 cos B 


b O — sin B cos B 
n sin 6 — cos @cos 8 — cos @ sin B 
Vo COs a@ O sin a 


The free-stream velocity vector Vy lies in the x’z’- 
plane and makes an angle a (angle of attack) with the 
positive x’-axis. The angle » between the velocity 
vector Vp and the normal n is given by 


cos 7 = cos (Vo, n) = cos a sin @ — 


sin acos@sin B (3) 


The condition for Newtonian flow is imposed by 
specifying that the gas stream upon striking the sur- 
face loses all of its momentum in the direction normal 
to the surface. Since the component of Vo normal to 
the surface element dS is Vp cos n and the rate of mass 
flow striking the element is po Vo cos n dS, the rate of 
change of momentum on the surface element in the 
direction of its normal is 


Vo cos n X poVo cos n dS = poVo? cos? n dS 


Thus, the excess local pressure force df produced by 
the momentum change is 


AL 
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df = (p — po) dS = poVo? cos? n dS 
20 cos? n d§ 
and the local pressure coefficient is 


cos? 7 2(cos @ sin § — 


C, = (pb — fo)/g = 2 
sin a@ Cos @ sin 8)? 


With respect to body axes, the forces on the body n; 
be separated into a normal force N in the 2-directi 
and an axial force X in the x-direction (Fig. 3). For, 
element of area which “‘sees’’ the flow, the force cop 
ponents are 


dN = —q@C, sin B cos 9dS = —qoC,r sin B dB dy 
since dS cos 06 = r dB dx and 
dX = gC, sin 0dS = qC,r tan 6 dB dx 


The total force is obtained by integration over the 
face of the body. In general, C, will be a function 


both x and ~. 


Consideration must now be given to those portions 
the body surface which are inclined away from the fre 
stream direction Vj and which may therefore be thoug 
of as lying in the “‘shadow”’ of the free stream. Ti 
situation is illustrated in Fig. 5. As the flow procee 
over the body, there will exist a boundary ac determin 
by the condition C, O (that is, b = py»). All 
the body surface situated upstream from ac is expos 
to the oncoming flow that, upon striking the surfa 
undergoes compression according to Eq. (5). All 
the body surface downstream from ac is in a region (| 


expansion flow for which Eq. (5) has no meanix 
Along the boundary ac, the tangent vector ¢ and tl 
free-stream vector V) both lie in the same plane—t 
tangent plane—and, consequently, nm and Vy are pe! 
pendicular. This condition defines ac and, from F 


(3), leads to the relation 
cos n, = cos a sin 0, — sin a cos 6, sin B, = 0 


( 


r 
sin 6, = (tan 6,/tan a) y 


where the subscript ( ), refers to conditions along t! 











B=B,(x) 


SHIELDED PORTION OF SURFACE 











Diagram to illustrate shielded portion of body in Net 
tonian flow. 
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boundary ac, which defines the limit of the compression 
nt as | 


flow area. 

For all transverse sections (sections normal to the 
y-axis) from the nose back to the section ab, the limits 
integration for 6 are from —7/2 to +7/2. Down- 


of 1 re 
stream from the section ab, the uppef limit for 6, 8, 
must correspond to the points lying on the boundary 


ac and will be a function of x. The point c, desig- 
nated by x = 2X, is the last pr int on the body to inter- 
cept aly of the and at this point 
= —m/2 (Fig. 9). 
a pointed body may be assumed to be conical over a 
Let the semivertex angle of the cone 
When a < 6,, the point a is situ- 
m/2) at the point 


flow, 
The extreme forward tip of 


free-stream 


short distance. 
tip be denoted by @,. 
ated along the top meridian (8 = 
where @ = a. When a > 6,, the point a is situated at 
the beginning of the body (x = 0). 
), two.cases must be distinguished: 


Thus, in integrat- 
ing Eqs. (6) and 
(1) a < 6,, some transverse sections are completely 
exposed to the flow with 6, = 7/2; (2) a > 8@,, the 
transverse sections are only partially exposed to the 
flow with 8, = sin! (tan @/tan a) 

Concerning the pressures on the shadowed or shielded 
portions of surface lying in expansion flow, which may 
be denoted by p,, little can be said except that 0 < 
pb. < po. If it should be assumed that the flow is com- 
pletely separated over the shielded regions, it would then 
0 for 
general relation for the 


po, which gives C,, 
The 


be appropriate to use p, 
the pressure coefficient. 
pressure coefficient is* 

(9) 


Coe = (2/y¥ Mo") [(Pe/Po) — 1] 


* The value y = 1.4 should be used for expansion flow. 


_ ra. lz a /2 
N= —2q0 / | Cor sin B dp dx — 240 f J 

0 —n/2 0 B, 
Cpr tan 6 dB dx + 2qQ f° '¥ C,.r tan 0dB dx 


coon ff 
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OF REVOLUTION 679 
If it is assumed that separation does not occur, p,/ Po 
becomes small compared to unity as Mp) increases, and, 
when My = both p,/po) and C,, are zero. Since 
Cy, = 0 for either of these possible extremes, the 
shielded portions of the surface would contribute noth- 
ing to the integrals of Eqs. (6) and (7). 

Thus, for true Newtonian conditions corresponding 
to Wy = ©, the total normal and axial force on the 


body obtained from Eqs. (6) and (7) aret 


N = —2q t i . C,r sin B dg dx (10) 
uo B,, 
X = 24 J / Cpr tan 0 dB dx 


where C, is given by Eq. (5). The Newtonian results, 
Eqs. (10) and (11), can be used as hypersonic approxi- 
mations provided the Mach Number is sufficiently large 
When separation does not 


(11) 


(Mo = 15, for example). 
occur, C,, increases with decreasing Mach Number, and, 
as lower Mach Numbers are considered (My < 10 or 15, 
for example), it would be necessary to allow for values 
of C,, different from zero because the pressures on the 
expansion areas begin to have an appreciable effect on 
the lift. On the for decreasing Mach 
Numbers, the value of C, given by Eq. (5) for the pres- 
sure on compression areas is less than a gas dynamic 
value and can therefore still be used as a conservative 
estimate for the flow. In this Mach 
Number range, an approximation to the aerodynamic 


other hand, 


compression 


forces could be written in the form 


+ Although no mention has been made of the base pressure 


coefficient, when My = © it vanishes in the same fashion as 


Cr 


C,,r sin B dB dx (12) 


(13) 


where it is understood that C,, may be put equal to zero without introducing appreciable error when the Mach 


Number is high enough or when separation occurs. 


Before obtaining the integrals (10) and 


eral expressions for the local forces on a transverse section. 
2r,z is the maximum diameter of the body. 
Considering only the first term in Eq. (10) corresponding to the surfaces 


defined by Cy = N, gotr yz”, where dg = 


x-axis will be expressed in units of x/d,. 


(11) for a particular body shape, it is convenient to derive the gen- 


The normal force coefficient Cy for the entire body is 
Axial distance along the 


in compression flow, the local normal force coefficient per unit length, Cy’, is 


,_. dy l dN 
d(x dp) Gomr 2” d(x dz) 
where r is a function of x only. Introducing 
C, = 


Dp 


Irom Eq. (5), the integration yields 


C, sin B dB (14) 


4rfr [* 
T?R SJ 2/2 


2(cos a sin 6 — sin a cos 6 sin 8)? 
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CC,’ = : 4 s 
N , © Sjmilai 
7 Ip I Simil 


] w\ . , . 1, ; ’ a - | 
; : By, + 5) sin 2a sin 26 + cos B, | 2 cos? a sin? 6 — , sin 2asin 20sin 8, + 3 sin” a cos? 6 (sin? B, + 2) 


2 2 2 f SR qtisin 
Introducing an axial force coefficient Cy = X/qorr,”, the local axial force coefficient per unit length is Writt: 
~ , “2 my sal 
dCx l dX 4) u For a 
Cc,’ = : = = tan 0 C, dg I 
d(x/dp) Got? zp” d(x/dz) Tp x/2 “F result 
extrem 
where 6 is a function of x only. The integration yields it sim 
4 r T pletels 
Cy = tan 6 («, oa ) (2 cos? asin? @ + sin? a cos? 6) + cos B, (sin 2a cos 20 — sin? a cos? 6 sin By) 1 (Dy Cas 
Trp 2 
Distinction is now made between the two cases a < 6, and a > 6,. 
Casel. a<8,: Forthiscase 8 = /2and Eqs. (15) and (17) reduce to 
. ; “ae fr 
Cy’ = 2(r/ryz) sin 2a@ sin 26 r and 
and 
Cx’ = 4(r/rg) tan 0[2 sin? @ + sin? a(l — 3 sin? @)] 
¥ 6 ° ~ ) TY 
Case?. a>86,: Inthis case, Eq. (8), Th 
takin 
B, = sin! (tan 6/tan a) tance 
distat 
and 
cos 8B, = V1 — (tan? 0/tan? a) 
al 
and Eqs. (15) and (17) reduce to 
r ; [8, + (m/2)] I ‘ ) _ Usiny 
Cy = 4 cos” @ sin 2a pat tan 6+ — cos 6,(cot a tan’ 6 + 2 tan aye (2 
'rR l T Bye 
: r (8. + (x/2)] — — oe _ 3 ; { o . 
Cx’ = 4— tan 0 = [2 sin? 6 + sin? a(1 — 3 sin? 6)] + cos B, sin 2a sin 26; (2g (C 
'r | v Arr 7 
To 
Eqs. (18) through (21) are the general expressions for the normal and axial forces on local transverse sections cylin 
the body, assuming C,, = 0 on the shielded surfaces and neglecting centrifugal forces in the flow. The totd| whic! 


normal and axial forces are obtained by integrating these local values over the length of the body. In ordert| of th 
apply the equations to a particular body of revolution, the profile of the body shape is introduced by specifying 'S 
6 as a function of x. Of particular interest in missile aerodynamics is the body consisting of cone plus cylinder sutla 
According to the Newtonian approximation that has been outlined, the pressure forces on a body are determined) there 


entirely by the local value of #6; consequently, the force on any portion of a body may be evaluated separately, ani) = ¢ 

the total force may be obtained by addition. der) 
(B) Cone 

: ‘ ; ” 
For a right circular cone of length /p, base radius 7p, and semivertex angle 6,, we have I Is 
infin 

’p r r x 
My, =lyb = ; 6 = const. = 6,; == = 

tan 6, a ™ i. 


As before, two cases must be distinguished corresponding to a < 6, and a > 8,. 


Casel. a <0,—8, = 4/2: From Eq. (18), it follows that whet 
The 

: In/4p , ' , : (dC, 

Cy = Cy’ d(x/dp) = (lp/dp) sin 20, sin 2a W 

by t 

Since /p/dp = 1/(2 tan @,), this may be written imps 
cylir 


Cy = cos? 6, sin 2a (22) 


| Thu 
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LIFT ON 


Similarly, from Eq. (19), the axial ccefficient is 


Cy = 2sin? 6, + sin? a(1 — 3 sin? @,) 
It is interesting to examine the value of the initial lift curve slope for the cone. 
(dC,/da),-0 = (dCy/da),-9 = 2 cos? 8, 


For a slender cone (6, —~ 0) this reduces to the slender body result, dC,/da = 2 
sult that at small angles of attack the lift coefficient for a slender diverging body has nearly the same value at 


re 


BODIES 


extremely high Mach Numbers as at extremely low supersonic speeds. 
at small angles of attack the lift coefficient is essentially independent of Mach Number. 
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(23) 
It is simply 


(24) 


2. This gives the rather surprising 


This indicates that for a body of this type 
This conclusion is com- 


pletely borne out by the Stone-Kopal values for cone lift at small angle of attack shown in Fig. 1. 


Case 2. a > 0,—f8, = sin! (tan @,/tan a): 


By, 2) 
sla. + (x / 


Cy = cos? 6, sin 2a 


T oT 


and from Eq. (21) that 


= [Buct (4 2)] 


Tv 


Cx 


The center of pressure on the cone may be found by 
taking moments about the vertex. If @ is the dis- 
tance from the vertex to the center of pressure, this 


distance is determined by the relation 
x 
C N 7 d ( ) 
dp 


“lp/dp x 
dp P 
moment 0 dp 


i = Zé) 


normal force Cy 


Using either Eq. (22) or (25), it is found that 


a= (2/3)lp (28) 


C) Cylinder 

To determine the Newtonian pressures on a circular 
cylinder, consider an infinite circular cylinder in a flow 
which is inclined at an angle a@ to the longitudinal axis 
of the cylinder. In the case of a circular cylinder, it 
is evident that only the lower half of the cylindrical 
surface is exposed to the free stream (Case 2) and, 
therefore, that 6, = 0. Also, since 0 = O andr = rz 
= const. = rs (the subscript S is used to denote cylin- 
der) for a circular cylinder, it follows from Eq. (20) that 


Cy’ = (16/37) sin? a (29) 
If /s is used to denote the length of any portion of the 
infinite cylinder, the normal force coefficient is 

5.33 Ls 


; sin? a = sin® a 
37 dg Tv ds 


: N 16 Ls 
( \ =>= = pe. 
Jot’ s 


(30) 


where V is the normal force on the cylinder length /s. 
The axial force on the cylinder, of course, is zero. Also 
(dCy/da),-9 is zero. 

When the effect of centrifugal forces as influenced 
by the boundary layer is considered in addition to the 
impact (Newtonian) forces, the normal force on the 
cylinder is reduced by 10 per cent (see Appendix). 
Thus Eqs. (29) and (30) become 


[2 sin? 6, + sin? a(l 


For this case, it follows from Eq. (20) that 


{ ed 
cos B, (tan 0, cot a + 2 cot 6, tan a)¢ (25) 
le 3 ; , = 
3 sin? 6,)] + cos 6, sin 26, sin 2a (26) 

dr 
Cy’ = (4.8/2) sin? a (31) 

and 

Cy = (4.8ls/mr ds)sin® a (32) 


The analysis, including centrifugal forces, indicates 
that such effects are smaller for conventional slender 
noses, such as cones and ogives at moderate angles of 
attack, and, thus, the pressure forces on such slender 
noses are satisfactorily approximated by the New- 
tonian (impact force) Consequently, the 
aerodynamic characteristics of a cone as given by Eqs. 
(22) through (28) are not modified. 


method. 


(D) The Cone and Cylinder Areas Situated in Expansion 

Flow 

All of the expressions (14) through (32) can be inter- 
preted in two ways. For truly Newtonian conditions 
corresponding to My = ©, it was pointed out previ 
ously that C,, = 0 and p,/p) = 0. In this case, Eqs. 
(14) through (32) give the total forces on the entire 
body, including the surfaces situated in expansion flow. 
If Eqs. (14) through (32) are considered as hyper- 
sonic approximations, they may be thought of as 
applying only to the forces on the surfaces in compres- 
sion flow. In this case—as discussed in connection 
with Eqs. (12) and (13) 
a nonzero value for the pressure coefficient C,,, where 


Coe = (2/y¥Mo?) [(pe/Po) — 1] 


corresponding to the surface areas situated in expansion 
flow. Although C,, is certainly negligible at extremely 
high Mach Numbers, it would have to be given con- 
sideration in any attempt to extrapolate the hypersonic 
approximations down to lower Mach Numbers. In 
view of the extremely approximate nature of such a pro- 
cedure, it would probably be sufficient to use an aver- 
age value p,, which is independent of 6 but which may 
be different on the cone and on the cylinder. 


it is appropriate to consider 
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flow on the lift of a yawed body of revolution. 
16 , ; 1 
c, - NORMAL FORGE 
. Qo As 
A,= CROSS-SECTIONAL AREA OF CYLINDER 
14 -— “* FREE-STREAM DYNAMIC PRESSURE 
© DATA FROM REFERENCE 7 
12 SE ——— —— 
z 
P | 
° a:6) 
a 0 ——_—-——— 
° Pa 
rs 4 
uw 7 } 
w Po 
° 
7 
° ( — 
a Fa | 
a] -4ae 
=z a=4 
= .06 | --e} ——_+—-—- | — 
= ot 
4 P| | 
= av aif 
- l [a 6" 
04 : | _ 
a tm” | 
-- | 
een HYPERSONIC | 
~~ 
on eco | | APPROXIMATION). e 
| 
| | 2:2") 
re) | ib l <a 
4 18 &2 2.6 3.0 3.4 3.8 
MACH NUMBER, M, 
Fic. 7. Normal force coefficient for a circular cylinder situated 
behind an ogive—the German A4 missile. 
(1) Cone.—If an average pressure is used, from 


Eq. (12) the normal force V, on the expansion flow areas 
of the cone would be written 


4qo SD koe 
N, = q - (1 «te yf / r sin B dB dx 
7M? Po B 


where p,, is the average pressure on the shielded area 
of the cone and 6, = sin~! (tan @,/tan a). The corre- 
sponding normal force coefficient Cy, is 
s “ - 
—_ Ne _ 2[1 — (Pep/Po)] 
Ne — = 


= V cot? 6, — cot? a 
got’ pp” 


Ty M,? 
(33) 
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Similarly, the axial coefficient is found to be data f 
2[(Pen/ po) — 1] (x _ 
ryM? 2 Dre it 
(2) Cylinder.—For the cylinder (since 8, = 0, 7. hes “ 
const., and it is assumed that p = p, = const.), the, _ 
sult is mee 
enoug. 
C goCpls ds om Als [(Pes/ Po) i 1] - to the 
" qorrs? oars yM? PF attack 
, menta 
In order to use these equations, approximate values jy aie 
p., and p, must be estimated inany manner that appeer it as 
feasible, possibly from two-dimensional gas dynam relati¢ 
and low aspect ratio supersonic wing theory. An inj 
cation of the maximum effect of the expansion pre 1 
sures on the lift of a yawed body is obtained by puttin a 
Pep be; = 0. This is shown in Fig. 7 for a typic; Fro 
body. The actual lift (normal force) could be expecte sonic 
to be somewhere between the two curves shown in Fy cone 
6. tion \ 
evel < 
(III) THe CorRELATION OF EXPERIMENTAL Resuy "YS ‘ 
FOR THE SUPERSONIC LirT ON BopIEs OF ReEvo.vurm) “PPT? 
values 
The previous discussions have been concerned @  ogive. 
tirely with approximate methods for predicting the |i indica 
of a body of revolution at hypersonic Mach Numbers tribut 
the order of 10 or 15 and above. Since it is importanf attach 
to have values of body lift throughout the comple{ straig! 
Mach Number range, this still leaves the problem of) begim 
estimating body lift at all lower Mach Numbers. Bef boat 
cause of the lack of theoretical results applicable t) ward’ 
this range of Mach Number (2 < Jy) < 10), a study diame 
has been made of available experimental data, and a! sectio 
attempt has been made to use the indications of thes) */ds : 
results to estimate (interpolate) the body lift in th) 74, 
intermediate range of Mach Number. For this pu site 
pose, use is made of the hypersonic approximation fo aailee 
the upper end of the Mach Number range. This po) 7p... 
cedure is admittedly almost qualitative in some 7 tied 
spects, and it is expected that at least some of the te Bin Fig 
sults given here will have to be modified as more theo} dea 
retical and experimental data become available. citi 
The amount of available systematic supersonic &) 6 at 
perimental data for yawed bodies of revolution appeat need 
to be extremely limited. A survey of all supersomi® ) 4. 
data on yawed bodies of revolution shows that ther) 4, 
exists no complete systematic series of tests (at leas§ that | 
with data in usable form) in which body lift is detety .4.. 
mined as a function of angle of attack, Mach Number) 4., 
and body fineness ratio—particularly for varyilg value 
lengths of cylinder behind the same nose shape. EX? 4; hi 
cept for the complete tests on the A4’ and the Wasselh thece 


fall,’ there is a great scarcity of pressure distributiog j.44:, 
data for yawed bodies. This type of information 4 a 
essential if significant comparisons are to be made witl sugge 
any theoretical results. Even for so simple a shape 2% yy, 


a cone, there are but fragmentary pressure distributiog 
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LIFT 


data for the yawed condition. Moreover, most of the 
atic experimental work in the past has been 


system a 
or less.!° 


p limited to body fineness ratios (l,/d,) of 7 
| The importance of tests on bodies of large fineness ratio 


s in the indications of the discussion below that the 


+ Ors |: on the cylindrical part of the body behind the nose 
t.), the Te not only becomes appreciable, Fig. 7, but also—at high 
enough speeds, Eq. (32)—the lift is directly proportional 
1] to the cylinder length and to the square of the angle of 
- fo attack. W hile much is to be desired by way of experi- 
mental data, by using the hypersonic approximation 
Values ji as a guide together with available experimental results, 
at appear it has been possible to obtain a rather consistent cor- 
dynamic relation for body lift over a wide range of conditions. 
An ind 
HON pres (A) Lift ona Cylinder Following a Cone or Ogive 


VY puttin; 


trom the close agreement in Fig. | between the hyper- 
r § 5 » 


a typied 
expecta! sonic approximation ands the Stone-Kopal values for 
vn in Fe, cone lift, it is evident that the hypersonic-approxima- 
; tion values give a good approximation for cone lift 
even at low supersonic Mach Numbers. It is instruc- 
RESULN tive to compare, at low Mach Numbers, the hypersonic 
OLUTIO approximation for cylinder lift with the experimental 
values for the lift on a cylinder following a cone or 
‘ned @  ogive. It is possible to extract some rather definite 
g the lit indications in this regard from the complete pressure dis- 
mbers¢) tribution data for the German A4 missile at angle of 
nportan§ attack. The A4 body has an ogival nose followed by a 
omplee straight cylindrical section that extends back to the 
yblem «i beginning of the tail surfaces, beyond which there is a 
rs. Bel boat tail. If x is the distance measured from the for- 
cable |) -ward tip of the nose and dg denotes the maximum body 
a study) diameter (dy, = ds, the diameter of the cylindrical 
and a section), the straight cylindrical section extends from 
of thes), x/dy = 3.5 to 6.0, giving a cylinder length of 2.5 cal. 
in the ‘ . a , 
¥ The data in reference 7 give the local normal force 
~ i coefficients Cy’ as a function of position x/d, along the 
we 7 body for a range of Mach Number and angle of attack. 
his pre These local values of Cy’ have been integrated over the 
me TS cylinder section and give the cylinder lift values shown 
be in Fig. 7. These results show that at supersonic speeds 
> theo 


) the cylindrical portion of a body contributes appre- 
ciable lift—which is rather strongly dependent on angle 


Nic eX 
of attack and Mach Number and emphasizes the ex- 
ip peals 
.. tremely poor approximation given by the slender 
ersonit ’ ; Mis i 
© body (zero-order) theory—which predicts zero lift on 
t there h op Je ; a , 
' » the cylinder, and by the linearized (first-order) theories 
t least ' ie oe eee a nee 
vail that predict only limited cylinder lift. Fig. 7 also indi- 
deter: ° 7 if: “6 : 
be cates that for low supersonic speeds the lift on the cylin- 
amber ; mA 
_ | der has already exceeded the hypersonic-approximation 
aryings | ae ban é , 
- value and is increasing rapidly with Mach Number. 
. ; 


‘jon | 
e wi 
upe a 


yu ti0! 


: ) At high enough Mach Numbers, yet to be determined, 
‘assets : 
bution 


these curves must decrease to the hypersonic-approxi- 
mation values. In view of the constant value of the 
normal force coefficient for a cone, for example, this 
Suggests that the normal force coefficient for a cone plus 
cylinder body must go through a maximum with respect 
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Schematic diagram to illustrate the effect of increasing 


Fic. 8. 
distribution of Cy’ along a cylindrical 


Mach Number on the 
section behind a cone. 


to Mach Number. It will be seen later that this is 
exactly what is indicated by the analysis of available 
experimental data. 

The qualitative notions concerning the lift distribu- 
tion on a cylinder following an ogive are indicated by 
At low supersonic Mach Numbers (Mo = 2 to 
’ decreases with 


Fig. 8 
3), the local normal force coefficient Cy 
distance downstream along the axis of the cylinder 
but gives an integrated lift coefficient Cy greater than 
the hypersonic-approximation value. At some higher 
Mach Number, probably in the range 3 < My < 6, 
Cy reaches its highest value. At still higher Mach 
along the cylinder axis 
approaching a constant 


Numbers the variation of Cy’ 
becomes smaller and smaller, 
condition, and Cy approaches the hypersonic-approxi- 
mation value. The Mach Number at which the total 
body lift becomes approximately equal to the hyper- 
sonic-approximation value will be discussed below on 
the basis of the indications of the experimental data. 


(B) The Initial Normal Force Slope for Missiles with 
Cylindrical Afterbodies 
A typical missile body is shown schematically in 


Fig. 9. It consists of a diverging section (nose) of 
length /p, which may be either a cone or an ogive; a 
of length /s; and a con- 


The semivertex 


straight (cylindrical section) 
verging (boat-tail) section of length /¢. 
angle of the nose cone, or of the inscribed cone when the 
nose is ogival, is denoted by @,. The total body length 
is /,, the maximum body diameter is dz, the maximum 
cross-section area is Az = (2/4)d,”, and the diameter 
of the base is d,. All of the body aft of the nose, or 
forebody, may be referred to as the afterbody, and its 
length may be denoted by /4(= ls + lc). The after- 
body may consist of all cylinder, part cylinder, 
part boat tail, or all boat tail. In general, the vari- 


and 
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@, = SEMI- VERTEX ANGLE OF INSCRIBED CONE 
d, = MAXIMUM DIAMETER OF BODY 

A, *-- d= MAXIMUM CROSS SECTION AREA OF BODY 
€ = BOATTAIL ANGLE 

d, = BASE DIAMETER 

4, = TOTAL LENGTH OF BODY 

4, = LENGTH OF NOSE, OR DIVERGING SECTION 

2, = LENGTH OF STRAIGHT OR CYLINDRICAL SECTION 
4, = LENGTH OF CONVERGING SECTION 












































Fic. 9. Diagram of typical missile body—schematic. 
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Fic. 10. Initial normal force slope for missiles with cylindrical 


afterbodies. 


ables upon which the normal force coefficient will 
depend may be indicated by the relation 


be d, 
= func. | 6,, >, Mo, a (36) 


dz dp 


C normal force 
y=: 
QA pz 


For the case in which the entire afterbody is cylindrical 
(no boat tail), it follows that 14 = ls, dg = ds, and the 
dependency of Cy becomes 


: normal force ' ls : 
Cy = = func. 1:0,, —, Mo « (37) 
QoA B ds 


Since most of the available lift data are restricted to 
small angles of attack, such results are employed to 
greatest advantage when they are used to evaluate 
(dCy/da),-». From Eq. (24) and Figs. | and 2, it ap- 
pears that the normal force on the conical nose varies 
approximately as cos? 6,, and this suggests use of the 
parameter (dCy/da),~»/cos’ 0,. However, most of the 
test data fall within the range 5° < 6, < 15°, so the vari- 
ation of cos? 6, from its median value is less than 4 per 
cent and within the accuracy of the experimental data. 
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Furthermore, an investigation of the lift on the cylinde 
immediately following a cone (using yawed cone datg 
and Prandtl-Meyer expansions) indicates that th 
cylinder lift may vary with cone semivertex angle go, 
to counteract the variation of cone lift with 0. Ry 
these reasons 6, is eliminated and the experiment; 
correlation is reduced to the form 


(dCy/da)a “9 = func. (Le ds, M)) (38 


This correlation has been carried out and yields th 
While the experiment; 
points (not indicated) show considerable scatter, ¢}, 
correlation is sufficiently good to define individy 
curves for cone and ogive nose shapes and is consider. 
fairly satisfactory, at least for preliminary PUrposes | 
It is found that when the length of the afterbody q.’ 
ceeds about 3 diameters (calibers), as far as total lift si 


results* shown in Fig. 10. 


concerned it makes little difference whether the nog 
At the lower Mach Number 
the curves go through a maximum with respect to bot}! 
afterbody fineness ratio and Mach Number. As [./( 
becomes large (>°9), the initial normal force slope be 


is an ogive or a cone. 


comes essentially independent of afterbody length 
At the higher Mach Numbers, Wo > 9, the initial Slope 
becomes equal to that given by the hypersonic approxi 
mation for a cone. This follows from the fact that 

cylindrical afterbody contributes nothing to the initial 
normal force slope in hypersonic flow [Eq. (30)]. At 
hypersonic speeds the ogive normal force is approxi-f 
mately equal to the normal force on its inscribed cone 
as has been assumed in Fig. 10. 





(C) Additional Normal Force Resulting from Angle oj 

Attack 

If the normal force curve slope were independent o/ 
a, the normal force could then be obtained immediate) 
from Fig. 10. However, except for small angles 0 
attack and missiles with no afterbody, the normal force | 
is found to depart widely from a linear variation with « 
and to vary in a manner that is more nearly a quadrati 
function of a. This effect of angle of attack on normal | 
force can be studied by introducing an increment a 
normal force coefficient, ACy, defined by 


Cw = (dCx da)a=0 x Qa + ACy (39 


This is illustrated in Fig. 11. Eq. (39) is the relation 


used to correlate angle-of-attack conditions, including 


the effects of Mach Number and afterbody length al 
angle of attack. It appears that, within the accuracy 
of the data, ACy varies as sin? a; therefore, the re- 
sults of the correlation are presented in the form ACy + 


* It will be noted that the abscissa in Fig. 10 is labeled /4/dg 
and therefore is not limited to missiles without boat tail. Its 
pointed out in Section (D) that the effect of boat tail may be ac 
counted for as an increment in initial normal force slope (se¢ 
Fig. 15) which must be added to the values obtained from Fig 
10. When there is no boat tail, 14/dg = 1s/dg and the values 
of Fig. 10 apply directly. 
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sin? a. Before giving the results of this correlation, it is 

worth while to point out how the hypersonic approxi- 
mations may be used to extrapolate from the limited 
range of conditions covered by the experimental re- 
sults to conditions of higher Mach Number and angle 
of attack. 

Ry using the experimental correlation results such as 
those of Figs. 10 and 11, it is possible to obtain experi- 
mentally based estimates of Cy up to Mo = 4, over a 
range of values of l;/ds and a. These are shown by the 
left-hand end of the curves (My, < 4.0) in Figs. 12 and 
13. At the high Mach Number end (My) ~ 15 or 20 
and higher), we have the hypersonic-approximation 
values that are independent of Mach Number. As 
explained in Part (II-D), by extending the hypersonic 
approximation to lower Mach Numbers using p, = 0 
for the surfaces in expansion flow, an indication is ob- 
tained of the limits within which the actual value of 
Cy must lie. These limits are shown as curves 1 and 2 
in Fig. 12, for example. Since there are no experi- 
mental data for fy) > 4.31 and since most of the experi- 
mental data do not extend beyond M/, = 3, these limit 
curves may be used as a guide for extending the low 
Mach Number range into the hypersonic range. This 
is the procedure that has been followed to obtain the 


results shown in Figs. 12 and 13. While this procedure 
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Diagram to illustrate typical variation of normal force 
coefficient with angle of attack—schematic. 
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Fic. 15. Decrease in initial normal force slope due to boat tail. 


admittedly contains certain elements of arbitrariness, 
it is believed to yield fairly realistic results and at least 
gives preliminary working values for lift over a range of 
conditions for which no other information exists. 

By employing results of the type shown in Figs. 12 
and 13, the correlated values of ACy (Fig. 11) from 
test been extended into the hypersonic 
region as shown in Fig. 14. The hypersonic approxi- 


data have 


mation for cones gives Cy as practically a linear func- 
tion of angle of attack for the low-angle-of-attack 
range (a < 26,) which is of most interest. The meager 
experimental cone data indicate that this is also true 
at low supersonic speeds. Thus, the additional nor- 
mal force coefficient AC, for a cone is zero, at least to 
the same degree of accuracy as can be expected from the 
experimental data. Consequently, it follows that the 
ACy values in Fig. 14 refer to the cylinder only and, 
for Mo 


©, are given by Eq. (32). 


(D) Effect of Boat Tail 


A preliminary correlation of available data has shown 
the effect of boat tail on normal force to be limited to 
the initial normal force slope, (dCy/da)a—». No con- 
sistent, pronounced effect on the additional normal 


force, ACy, was found within the angle-of-attack range 





STREAMTUBE 


Diagram of stream tube on surface of body. 


Fic. 16. 


SCIENCES NOVEMBER, 
of the tests (a < 10°). 
at higher angles of attack (a > 10°) would show some 
effect of boat tailon ACy. It appears that the initial 
normal force slope decreases linearly with decreasing 


It is likely, however, that dat, 


base ratio, d,/dz, over the normal range of base ratios— 
0.4 < dy dp A 1.0. 


The effects of boat-tail angle ang 
type (conical or ogival) and the body preceding the boa 
tail appear to be small compared with the effect of the 
The results of this correlation are showy 
by Fig. 15. It will be noted that the boat-tail effect jg, 
maximum at low speeds and decreases to zero at aboyt 
My 9. This is in agreement with the slender-body 
theories'' that predict zero lift when d,/d, 


base ratio. 


0 and 
with the hypersonic approximation that indicates po 
boat-tail effect on the initial slope parameter, |; 
should be mentioned, however, that the hypersonic 
approximation would show an effect of boat tail on AC, 
since the boat-tail portion of the body would be com. 
pletely shielded at angles of attack less than the boat. 
tail angle. 


(E) Center of Pressure 


By using the normal forces that have been estimated 
for any cone-cylinder or ogive-cylinder combination, it 
is possible to estimate the corresponding center of pres- 
sure. Figs. 10 and 14 give the normal force on a cone 
or ogive and the normal force distribution along a 
cylinder following the cone or ogive. The center of 
pressure of the cone and ogive are approximately 0.67], 
and 0.53/p, respectively, aft of the nose. Fig. 15 
gives the increment in normal force due to boat tail. 
It appears satisfactory to assume that this force acts 
at the midpoint of the boat tail. The location of the 
missile center of pressure aft of the nose tip, Xp, 
is obtained by forming the summation of the various 
components according to the formula 


Yep. = > Cues > Cw, 


Preliminary checks of this method have shown satis- 
factory agreement with experimental data. 


(40) 


Appendix 
ANALYSIS OF THE CENTRIFUGAL FORCE EFFECTS 


General Discussion 


The simple analysis of Newtonian flow is incomplete 
inasmuch as the effects of centrifugal forces in the flow 
around the body have been neglected. In the flow 
over plane surfaces at angle of attack and on cones at 
zero angle of attack, since the paths of the air particles 
over such surfaces (that is, the stream lines on the sur- 
faces) are straight lines, no centrifugal forces are pres- 
ent. However, when the are 
curved (as is the case for bodies of revolution at angle 
of attack, for example), centrifugal forces will be pres- 
For these flow problems, the total 


surface stream lines 


ent in the flow. 
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surface pressure coefficient at any point on the body is 
al to the impact pressure coefficient minus the 
centrifugal pressure effect.'2 3 If Cy, = (pi — po)/go 
is the pressure coefficient due to the Newtonian impact 
pressure P; and if p, is the pressure relief due to centrif- 


ugal f 


equ 


orces, then the net pressure coefficient C, is, simply, 
C, = [(bi — be) — pol/qo (41) 


where the net pressure p is specified by p = p; — px. 
Also, when centrifugal forces are present, the limit 
angle 6, denotes the point of zero net pressure co- 
efficient and not the point of zero impact pressure 
coefficient, as in the earlier discussion. 

The Newtonian impact pressures are evaluated ac- 
cording to the methods given in Part (II). The pres- 
sure relief p, resulting from centrifugal forces is evalu- 
ated from the formula 


Cas — (dD, go) = 


pe =*mV, /RAC (42) 
where 

jn = rate of mass flow through a stream tube at 
any point on the surface 

V,, = effective velocity in a stream tube at any 
point on the surface 

R = radius of normal curvature of the stream tube 
at any point on the surface 

AC = width of the stream tube (the height of 


stream tube is the body layer thickness de- 
scribed below) 


In order to evaluate mm and R for a body of revolution 
at angle of attack, the particle paths, or stream tubes, 
onthe body surface must be determined. In the hyper- 
sonic approximation, the shock wave may be imagined 
to wrap itself around the portions of the body which 
are subject, or exposed, to compression flow. On these 
portions of the body, the flow is confined to a thin, 
high-ensity layer that lies on the surface of the body 
and which, therefore, may be referred to as “‘the body 
layer.’ Neglecting friction, the total reaction of the 
body-layer fluid on the surface and of the forces acting 
on the fluid must be normal to the surface. Hence, 
the principal normal (normal radius of curvature) and, 
therefore, the osculating plane at every point on a 
stream line must be normal to the surface. Thus, the 
stream tubes are similar to the geodesic paths obtained 
for the constrained motion of a particle on a curved 
surface.'* At each point along a stream tube, the 
radius of curvature is directed along the inward normal 
to the surface—that is, along the vector n. The rate 
of mass flow m at any point on a body-layer stream tube 
is obtained by finding the sum of all the particles that 
have previously entered the body layer along the 
particular stream tube. 

Let the curve C, Fig. 16, denote a stream tube lying 
on the surface of the body, and let the radius of normal 
curvature of the stream tube at the point P on the sur- 


face be denoted by R. Let the lines of curvature at P 
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have the directions given by the unit vectors ¢ (maxi- 
mum radius of normal curvature) and 6 (minimum ra- 
dius of normal curvature), and let the curve C (that is, 
the stream tube) at this point make the angle y with ¢ 
(Fig. 16). According to Euler’s theorem"™ it follows 
that 


1/R = (cos? y)/Ri + (sin? y)/Re (43) 


where XR, is the radius of principal normal curvature 
in the ¢ direction, and Rs, that in the bdirection. Fora 
body of revolution, which will be the case treated here, 
t is tangent to a meridian and 0 is tangent to a circular 
parallel; the meridians and parallels are the lines of 
curvature on the surface."® R; is the radius of curva- 
ture of the normal (meridian) section obtained by pass 
ing a plane normal to the surface and containing the 
axis of symmetry. RR, is known immediately when the 
profile shape of the body of revolution is specified. R» 
may be determined by means of Meusnier’s theorem 
(see page 505, reference 15), which shows that 


R. = r/cos 6 (44) 


where ¢ is the radius of curvature of a circular parallel 
and @ is the angle between n and r (that is, @ is the direc- 
tion of the tangent along a meridian). From Eqs. 
(43) and (44), it follows that 


1/R = (cos? y)/Ri + (sin? y cos 0)/r (45) 


is the relation for the radius of curvature of the stream 
line at any point on the surface of a body of revolution. 

Consideration must next be given to the velocity 
V,, of the flow over the surfaces of the body which are 
exposed to compression flow. When a free-stream par- 
ticle strikes the body layer at a local point P (see Fig. 
16), it loses its component of velocity normal to the 
surface at this point, while the velocity components in 
the tangent plane remain unchanged. Letting 
denote the angle between Vp and n [n, is identical with 
n in Eq. (3)], it follows from Eq. (3) that 


COS Mm» = Cos a sin 8 — sin a cos @ sin B 


The normal component of velocity Vo cos 7, is lost upon 
impact, and the fluid particle after impact at the point 
P is left with the velocity components along t and 6 un- 
changed. If the angle between V, and ¢ is denoted 
by 7, and the angle between V, and b by m, it follows 
from Eq. (3) that 


cos n, = cos a cos 6 + sin a sin @ sin B (46) 


and 


COS n, = sin a cos B (47) 


After impact, the particle is left with the instantaneous 
velocity components Vo cos 7, along f and Vo cos m 


along b. The resultant velocity vector is of magnitude 


Ves = Vo(cos? Ut + cos? ™). ’ (48) 


It lies in the tb-plane and makes the angle y with the 
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LIFT ON INCLINED 
vector ¢, where 
(Vo cos 72) | 
cos Y = * = = (49) 
Vs [1 + (cos m/cos n,)*] 
From the relation tan y = cos /cos 7, it follows 
that 


tan y = 1/(cot a cos @ sec 6 + sin # tan B) (50) 


Although a single particle after striking the body 
would continue its motion in the tangential plane 
formed by ¢ and 06, the actual flow of a continuous 
medium constrains the particles to follow a stream tube 
on the surface. In order to evaluate the effective body- 
laver stream-tube velocity, V;., a relationship must be 
’ or Vo. Five postulations 


found between V,, and V4, 


have been scrutinized. 


Case z. Vs, => Vo: 
sumption and would overestimate the centrifugal force 


This is the most simple as- 


effects. 

Case 2. V;. = Vs,: This yields the result that each 
particle does accelerate along the body layer, but the 
postulation does not give a velocity gradient in the body 
layer at each point. 


Case 3. Vi. = S'VsdS/fdS: According to this 
method, the velocity of a particle remains constant after 
impact and results in the existence of a velocity gradient 


in the body layer.!* #3, 17-19 


Case 4. V,. = V,;/2: This relation is the result of 
the assumption that the thickness of the body layer is of 
the order of magnitude of that of the boundary layer. 
The velocity distribution normal to the body surface is 
considered as a linear function of the distance from the 


surface throughout the body layer. 


Case 5. Vs. = f(0, B, a) Vs,/2: 
of the pressure distribution as obtained by Sauer” using 
the method of characteristics for the A4 nose at Mach 
Numbers of 3.24 and 8.00 (see Fig. 17), it appears that 
the local pressure coefficient should decrease slower than 
any of the above four cases and that the location of 


Upon examination 


zero pressure coefficient should occur at the position 
where the surface tangent is parallel to the free-stream 
Case 4 yields reasonable pressure 
variation at the start of the body layer, it is used as a 
starting point for a fifth case. The present fifth case 
is then a modification of Case 4, where f{(0, 8, a) must be 
unity at the origin of the body layer, varies as a quad- 
ratic along the body layer, and approaches zero at the 


direction. Since 


position where the surface tangent is parallel to the 
axis of the free stream. As an additional check on the 
validity of Case 5, all five cases are applied to flow over 
spheres (see Fig. 18). The aerodynamic drags are 
shown in Fig. 19 as the limiting values for the experi- 
mental data of Charters.?! the skin friction 
and base drags for spheres at high velocities are neg- 


Since 


ligibly small, Case 5 appears best. 
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Example of Application—Cylinder, Case 5 


The treatment of the hypersonic forces on a cylinder 
will be based on the assumption that the cylinder is 
In the case of an infinite cylinder, the sur- 
This condition makes 


infinite. 
face stream lines are all parallel. 
unnecessary the determination of the individual stream 
lines in the calculation of the centrifugal pressure effects. 
The rate of mass flow in the body layer at any particu- 
lar position on the surface of the cylinder, where the 
centrifugal pressure is desired, is 


m = poVoAcap = poVorls sin a cos B = ppipV,,ts sin 
(51) 


where /; is a specified length of the infinite cylinder. 
From Eq. (45), the radius of curvature R of a stream 


line on a cylinder is (since Ri = © and @ = Q for a 
cylinder) 
R = r/sin? y (52) 


The expression for the centrifugal pressure at a point 


on a cylinder is 


pe = mV,,/lsR sin y (53) 
Combining Eqs. (51), (52), and (53), the centrifugal 
pressure relation may be written 
po Vorls sin a cos B sin? y V 
D6 = f S¢ = 
f rls sin ¥ 
poVoVs, sin a cos B sin y (54) 


Employing Case 5, the effective surface velocity for a 


cylinder is 


V;- = (1/2) sin? B Vo (cos a/cos y) (55) 
since V’;, = Ve (cos a/cos y) as given by Eq. (49). Now 
Eq. (54) may be written 
Pe = poVo? sin acos 8 sin y(1/2) sin* B (cos a/cos y) = 

(1/2)po Vo? sin a cos B sin? B tan y cosa _ (56) 


For a cylinder (6 = 0), Eq. (50) yields 
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tan y = tan acos B (57) 
and Eq. (56) becomes 
Pe = (1/2)p0 Vo? sin? a cos? B sin? B = 
go Sin? a cos? B sin? B (58) 


For a cylinder, the impact pressure coefficient given 
by Eq. (5) is 


Cy, = 2 sin? a sin? B (59) 
and the net pressure coefficient is 
Cp = Cri — (P-/qo) = sin?’ B sin’? a(2 — cos? B) — (60) 
The equation defining 6, is C, = 0, or sin® 6, = 0. 


Thus, for a cylinder, 8, = 0° for all a. 

From Eq. (6) the normal force coefficient corre- 
sponding to the surface of the cylinder exposed to com- 
pression flow is 


; N lL} ae 
Cy = _-— = C, sin B dB = 
Qotrs~ TVs —n/2 
Zils .. (61) 
sin? 
m6 CS 


where rs is the cylinder radius and ds = 27g is the 


cylinder diameter. By comparing this result with 
Eq. (30) for the case in which the centrifugal forces are 


neglected it is found that 


(Cy) with centrifugal force (Cy) Newtonian = 0.900 (62) 


In Fig. 20, the effect of the various postulated V,. on 
the normal force for a cylinder is given. 
V4 
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The Longitudinal Stability of Elastic Swept 
Wings at Supersonic Speed 


C. W. FRICK* anp R. S. CHUBBt 
Ames Aeronautical Laboratory, N.A.C.A. 


SUMMARY 


The longitudinal stability characteristics of elastic swept wings 
experiencing bending and torsional deformations are calculated 
for supersonic speed by the application of linearized lifting- 
surface theory. A parabolic wing deflection curve is assumed, 
and the analysis is simplified by considering the wing panels 
to be simple cantilever beams. The method is thereby limited 
in application to wings of high aspect ratio for which the root 
effects are small. Expressions for the lift, pitching moment, and 
span load distribution are derived in terms of the elastic prop- 
erties of the wing—namely, the design stress, the modulus of 
elasticity, the shearing modulus, and the maximum design load 
The analysis applies to wings with leading edges swept 


factor. 
In all cases, however, the trailing edge 


behind the Mach lines. 
is sonic or supersonic. 

The validity of the assumption of a parabolic deflection curve 
is investigated by comparison with matrix calculations for a 
wing with the moment of inertia varying with the fourth power 
of the spanwise dimension. It is shown that for such wings the 
assumption of a parabolic deflection curve gives valid results 
for the wing stiffnesses and flight conditions in which the designer 
is primarily interested. 

Calculations made for subsonic speed, when compared with 
the results for supersonic speed, indicate that, except for the 
effects of the change in the chordwise position of the section center 
of pressure, the aeroelastic characteristics are primarily a func- 
tion of the flight dynamic pressure and not of Mach Number, 
at least for wings that have subsonic leading edges at all flight 
Mach Numbers. 


INTRODUCTION 


| ha THE PAST SEVERAL YEARS, the N.A.C.A. has 
been engaged in analytical studies of various wing- 
body combinations to determine their suitability for 
efficient supersonic flight. A great deal of this re- 
search has been concerned with wing-body combina- 
tions for which the wing is swept within the Mach cone. 
For such airplane configurations, Jones! has shown that 
lift-drag ratios of the order of 10 to 12 are possible at 
moderate supersonic speeds (./ 1.1 to 1.4). 

In view of the promising flight efficiencies of swept 
wings, the analytical investigation was extended to 
include the stability characteristics of these wings. 
The changes in stability due to the elastic deformation 
of the wing, especially, were thought to be of major 
significance, since it was recognized that the structural 
weight penalty necessary to overcome these adverse 
aeroelastic effects might easily outweigh the aerody- 
The present investigation was 


namic advantages. 


Presented at the Aerodynamics Session, Eighteenth Annual 
Meeting, I.A.S., New York, January 23-26, 1950. 
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therefore undertaken to provide estimates of the aero- 
dynamic effects of the elastic deformation of the wing 
panels. 

The analysis is concerned with aerodynamically efti- 
cient wings—that is, tapered swept wings of high aspect 
ratio with essentially elliptic span load distributions. 
For such wings of usual structural design, the ratio of 
the bending moment to the section moment of inertia is 
essentially constant spanwise so that the deflection 
curve is of parabolic form. The twist of streamwise 
sections therefore varies linearly across the span, and 
the aerodynamic loading due to the twist may be calcu- 
lated by lifting-surface theory. 

The extent of the present paper is limited to the 
treatment of longitudinal stability. Accelerated flight 
conditions only are investigated, since the longitudinal 
stability of the aircraft has greatest significance in 


dynamic flight. 


SYMBOLS 


Cartesian coordinates 


Xi, VM 

x,y = transformed Cartesian coordinates in terms of the 
semispan dimension, s 

fn = x, y coordinates of the apex of any superposed lifting 
sector 

S = distance in the y; direction from the root section to 


the intersection of the flexural axis and the tip 
Mach cone 

- = distance along the flexural axis from the root section 
to the intersection of the flexural axis and the tip 
Mach cone 

yy’ = distance measured from the root section along the 

flexural axis 
spanwise distance in y direction from the root section 


Y = 
to the center of load on the half wing 

Ss = wing area 

nN = taper ratio, ratio of tip chord to root chord 

Cas = average chord 

rs = mean aerodynamic chord (Sc? dy Rg c dy) 

r = local chord parallel to the plane of symmetry 

Co = root chord parallel to the plane of symmetry in 
terms of the span dimension, s 

A.R aspect ratio 

A = angle of sweepback of the flexural axis 

0 = slope of the flexural axis in a vertical plane passing 
through the flexural axis 

n = maximum load factor 

M = bending moment at any point on the flexural axis 

M, = bending moment at the root section of the wing 
beam 

T = torsional moment at any point on the flexural axis 

r torsional moment at the root section of the wing 
beam 

E = modulus of elasticity for the wing-beam material 
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G = shearing modulus for the wing-beam material L = lift 
I = moment of inertia of the wing beam Cc. = lift coefficient (L/qgS) 
5, = moment of inertia of the root section of the wing ce = lift coefficient at maximum load factor 
beam Ms = section pitching moment of a wing section about 4} 
ys = torsional stiffness constant apex of the wing 
a = torsional stiffness constant of the root section of the Cm = pitching-moment coefficient about the apex of ; 
wing beam wing in terms of the mean aerodynamic chord 
c = distance between the flexural axis and the center of the wing area 
pressure of the sectional lift in terms of the local Cra = the rate of change of lift coefficient with the ang 
chord, positive if the flexural axis is behind the attack of the root section 
center of pressure Cng = the rate of change of pitching-moment coefficjg 
Gana maximum design stress with the angle of attack of the root section 
d, = maximum thickness of the wing at the root section 
a = angle of attack of the root section of the wing, rad. ANALYSIS 
Qt = incremental angle of attack at any spanwise station 
of the wing, rad. Lifting-Surface Theory 
Qs = angle of attack of the wing section at any spanwise 
station, rad. For convenience, the aerodynamic effects of th 
ay, = angle of attack of the root section at which maxi- bending deformation of the wing are considered first! 
femme fees Sector is developed, rad. Bending.—The bending deformation of a swept-win; 
B sa ai — 1, where Mis the free-stream Mach Num- peam differs somewhat from that calculated by th 
m = B ies the cotangent of the angle of sweepback of usual simple cantilever-beam theory used for unswe 
the wing leading edge wings because of two differences in the wing structur 
mi = £6 times the cotangent of the angle of sweepback of | characteristics: 
the wing trailing edge (1) The swept-wing beam is restrained along 
t = £6 times the cotangent of the angle of sweepback of a streamwise root section and not along a section per 
ray from the apex of any superposed lifting sector : : : 
r = complete elliptic integral of the first kind with modu- pendicular to the beam centerline. 
ne (V1 ow galt (2) The swept-wing beam does not have a tn 
y = complete elliptic integral of the second kind with flexural axis. The line of section torsional center 
modulus (V1 — 9%) which usually defines the flexural axis has marked curva 
W = airplane weight ture near the root section. 
W/S = wing loading For wings of high aspect ratio, however, the accurac 
' ee ee es re ee ee bending deformations calculated by simple bean 
density and V the velocity of the free stream : ; ; et ‘ 
Ap/q = lifting pressure coefficient theory is not greatly impaired by assuming a straight 
] = load per unit span flexural axis exists coinciding with the straight portia 
of the line of section torsional centers. Further, for th 
purpose of analysis, as is shown in Fig. 1, the root sec 
tion of the wing beam is assumed to be the sectio! 
formed by extension of the wing-beam generators ont 
a plane perpendicular to the flexural axis and passin 
° through the intersection of the flexural axis and th 
streamwise root section. This simplification of thi 
beam analysis has some justification in recent expen 
ments with swept-wing beams performed by Zende! 
Plane of beam and wevtigd of the N.A.C.A. The length echt 
root section! beam s’ is taken as the distance along the flexural axis 
from the root section to the intersection of the flexura 
axis and the tip Mach cone. The semispan s of thi 
wing is taken as extending from the root section to the 
intersection of the flexural axis and the tip Mach cont 
in a direction perpendicular to the plane of symmetry 
The portion of the wing lying within the tip Mach cont 
is ignored, since, as has been shown by Cohen,’ little 
load is carried in this region. 
The coordinate system is selected as shown in Fig. > | 
x The origin of the coordinate system is placed at the ape | 
~~ of the wing, the positive branch of the x; axis lying 
% m a. downstream. 
Tip * The ameliorating influence of the distributed mass of t 
Mach lines wing is not considered, and the change in camber of the airfoil | 
Fic. 1. Geometric characteristics of the wing beam. sections due to the distortion of the wing surface is ignored. } 
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LONGITUDINAL STABILITY 

The mathematical treatment may be made less tedi- 
ous by transforming and nondimensionalizing the co- 
ordinates so that, in the following analysis, 


7 = BY1 S, xX = X/S 
Co = root chord/s 


In general, at both subsonic and supersonic speeds, 
selection of the wing plan form for low drag leads to a 
combination of spanwise loading and spanwise dis- 
tribution of the bending resistance in the wing beam 
such that the wing deflection curve is essentially para- 
bolic; the ratio of the bending moment at any spanwise 
station to the moment of inertia at that station is con- 
stant across the span. The deflection curve deviates 
appreciably from a parabola only if the aeroelastic ef- 
fects experienced by the wing are extremely large. 
Since the designer is primarily interested in wings for 
which the effects of the elastic deformations of ihe 
wing are small, it appears that the assumption of a 
parabolic deflection curve is admissible. This point 
will be developed further later in the analysis. 

Assuming that the deflection curve of the flexural 
axis is parabolic, the slope of the flexural axis, as shown 
in Fig. 3, is 


6 = (M/ED)y’ 


where y’ is measured along the flexural axis. 

The incremental angle of attack of streamwise sec- 
tions of the wing is related to the slope of the flexural 
axis as 


—6@sin A 


ay = 


The slope of the flexural axis in nondimensional trans- 
formed coordinates may be written as 














Fic. 2. Coordinate system for calculation of the characteristics 


of elastic wings. 
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Sketch of the deflection curve and angle of twist of the 
wing beam. 


6 = (M/EI) [s/(8 cos A)]y 


Fic. 3. 


The incremental angle of attack of any streamwise sec- 
tion of the elastic wing is then 


—(M/ETI) (s/8)y tan A 


ar = 
and the total angle of attack of any streamwise section 
is 


a, = a — (M/ED) (s/B)y tan A (t) 


where a is the angle of attack of the root section of the 
wing. 

The distribution of pressure over the elastic wing due 
to twist may be determined by applying known conical- 
flow solutions for supersonic flow. In the linearized 
theory, the principle of superposition of various solu- 
tions may be used to satisfy the particular boundary 
conditions of the problem. For the elastic wing, the 
flow field may be considered to consist of the super- 
position of two distinct flow fields: 

(1) The flow about a flat rigid wing at an angle of 
attack equal to the angle of attack of the root section. 

(2) The flow about a twisted wing for which the 
angle of attack at the root is zero. 

The solution for the first flow field is that given by 
Stewart‘! and Brown’; the second flow field can be ob- 
tained by determining the solution for a differential 
twist da, at one station and integrating this solution 
across the span. 

The solution for the pressure distribution at any 
point, corresponding to a differential twist, must meet 
the following boundary conditions (Fig. 4): 

(1) Outboard of the station of twist, the angle of 
attack must be constant and equal to the differential 
twist. 

(2) Inboard of the station of twist, the angle of 
attack of the surface must be zero. 

(3) Between the swept leading edge and the Mach 
cone, no lifting pressures may exist. 
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The conical-flow solution corresponding to these 
boundary conditions is that for a special lifting sector 
given by Lagerstrom® and is expressed in the notation 


of the present paper as 


Ap 8a m™®” i+? (2 
qd Br (m + 1) m—t 
where ¢ defines a ray from the apex of the sector. Fig. 4 


shows both a sketch of the boundary conditions to be 
met by this solution and a plot of the pressure distribu- 
tion given by Eq. (2). 

The induced pressure resulting from twist due to 
bending of the elastic wing may be found by integrating 
across the span of the wing. This integration corre- 
sponds to the superposition of an infinite number of the 
lifting sectors along the span, each sector having an 
infinitesimal angle of attack da,. 

The pressure due to twist is then given by 


Ap Ss m’* Ms “ i+ 
css —— tan A dy 
qd /r B'r (m+ 1) El 0 m—t 
where 
ie 7-9. wmy- 
 £— £ mx — 7 


The x and y coordinates of the apex of any superposed 
sector are &, n. 

The integration must be carried out from the root 
section of the wing 7 = 0 to the value of » = m corre- 
sponding to the last superposed sector, the Mach cone 
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Fic. 4. The boundary conditions and pressure distribution fo; 
the conical-flow soiution for differential twist. 
tion. The value of m is found by placing ¢ equal to - 


and solving for 7. 
no = [m/(m + 1)] (x + y) | 


The integration yields at any point x, y the pressur 
due to twist 


(**) =— = . tan A = (x + y) X | 
q/r 3B? (m + 1)? EI ; 
x ¥ : | 

mie — ¥ 7 


To this expression must be added the conjugate tem | 
due to the elastic deformation of the opposite wing 


panel. The conjugate term may be obtained by substi- | 
' 


of which encompasses the point x, y under considera- tuting —yfory. Then, 
! 
' 
(“*) 16 m*/? ne Ms behead x+y Po x—-y 
=— — — tan A — x+y x-y 4 
q/r 3B? (m + 1)? EI mx — ¥ mx + y 


It should be noted that the addition of the conjugate 
terms adds some small lifting pressure in the region 
between the wing leading edge and the Mach cone 
where no lifting pressure may exist. These pressures 
may be canceled by the superposition of constant lift 
sectors. Since these extraneous pressures, on the aver- 
age, amount to about 3 per cent of the average pressure 
coefficient over the adjacent wing surface and since 
elimination of those pressures would change the pres- 
sures over the surface only about one-half of 1 per cent, 


4m*a 


Ap 


ee tan A 
q BS Vm? (y/x)? 


 3B2e (m + 1)? 


Ms 
EI 


it seems that, in view of the additional complication 
involved, the cancellation of these pressures is un- | 
warranted. * 

The total lifting pressure for the elastic wing at an | 
angle of attack is then obtained by adding to Eq. (4 
the solution for the flat lifting wing. For the elastic 
wing, then, 

* The authors have recently noted that an exact solution for 
the case of symmetric twist has been obtained by Theodore 
Goodman, of Cornell Aeronautical Laboratory (CAL/CM 524 
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Examination of this equation shows that the relation- isa linear function of the angle of attack, 
ship between M/EI and a must be established before 


. 4° / = (9 é 4 ’ 
the pressure distribution can be calculated. Since for M/I = (2¢maz./dr) (a/ an) 


; Wi arabolic deflection curves the maximum . : ‘ F 
s with parabolic where Omar. is the design stress at maximum load factor, 


wing 
<< occurs < > point of maximum thickness . . ; ; . 

stress occurs at e oe ¥ ; ’ d, is the thickness of the root section, and a, is the 

sssalle > > maximum stress occurring at - : : 
usually the root, the 8 angle of attack at maximum load factor; an expression 

eg ad factor is : : . 
maximum load factor 1 for a, is derived later. 

Omar. = (M/TD),(d,/2) - , eae tt 
we 7 rhe equation for the pressure distribution may then 


and since the bending moment at any point on the span _ be written as 


4m? 32 m’* fui + ah. x + y R x—y 
a = : eras ~ tan A = (x + y) — + (x — y) : (6) 
q > Vm? — (y/x)? = 38?a (m + 1)? E d,a, mx— y mx + y 


The load per unit span can be obtained from an integration of Eq. (6) with respect to x along any streamwise 


(y + mtco)/mt 
l q = § f (Ap q)y = const. dx 
y/m 


station (vy = constant), 


The integration is carried out from the leading edge of the wing, x = y/m to the trailing edge x = (y + m,co)/m, 


and yields 


l 4m?sa fil 32 m'/* ‘ \ Cmer, S* O fol) (7 
me 7 = = , tan J =a (Vv 7) 
q Bz 1 3B°x (m + 1)? E d,a, 


The functions f;(y) and f(y) are given in the appendix since they are somewhat unwieldy. 
The lift coefficient may be obtained by an integration of Eq. (7), spanwise from root to tip. 


8B 
BC, = (2s/S) £ (l/q)dy 
0 


The integration yields 


: 25° | 4m?a 32 m’* Cus. $ @ . 
BC L= . : F, ee ee = tan A : 4 F, (8) 
B Bz 3B°r (m + 1)? Ed, Qn 


The constants F; and F2 are given by equations in the Appendix. 
Eq. (8) may be used to determine the angle of attack at maximum load factor a, which is needed in the fore- 


going equations: 


An 


B2> a 32 m’* ik Os 
= = C.. + , , tan A ; F, (9) 
4m?*F, (2s? 3B°r (m + 1)? ze «© 


The pitching-moment characteristics of the elastic wing may be determined by an integration of the pressure dis- 
tribution given by Eq. (6). 
For any spanwise station, the section pitching moment about the apex of the wing is 


(y + m,co)/mt 
m,/q = —s? (AD/@) yo conct, % GX 


” y/m 


This integration yields 


m 4m*a 32 m'?* Cute. S @ 
— = — s? 3(v) — tan A —— (y) 10) 
q | Br fy) 3B°x (m + 1)? Ed, a fc | ( 


The functions f2(y) and f;(y) are given in the Appendix. 
The total pitching-moment coefficient about the apex of the wing in terms of the mean aerodynamic chord is 


B 
BC,, = (2s/Sé) f (m,/q) dy 
0 


found by integration across the span, 
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: 2 s* | 4m*a _ 32 m’* Ciey Sb o. 
BC» one a | . Fs ee 9 tan A = FP, l] 
Se B> 3B°r (m + 1)? KE dot 
or 
me 2 53 | 4m?F3 32 m’* a oe 
Ss, = —=- es - tan A = 
- Sé Bz 3B°r (m + 1)? Ed, " 


The constants F; and F;, functions of the aspect ratio, taper, and sweepback, are given in the Appendix. 


Torsion. 
effects of wing twist due to torsion. 
tained are, in reality, those for wings of infinite tor- 
For wings having large angles of sweep- 


‘The previous analysis has ignored the 
The solutions ob- 


sional stiffness. 
back such as are necessary for efficient supersonic flight, 
the aerodynamic twist due to torsion has been calculated 
to be about 15 to 20 per cent of the twist due to bend- 
ing (for thin wings). Insuchcases, the effect of the tor- 
sional deformation on the spanwise loading may be 
Eq. (7) 
may be utilized in the calculation of the torsional mo- 


neglected in calculating the torsional moment. 


ment in this instance. 

An expression for the torsional moment at the root 
section of the wing beam (perpendicular to the elastic 
axis) may be obtained by assuming that the distance 
from the center of pressure to the flexural axis for any 
section of the wing is a constant percentage of the local 


Then, 
° B 
sf cos A 
70 


where c is the local streamwise chord given by the equa- 


chord. 
Z5 _ 
q 


c dy 
q 


8 


tion 
c = sto [1 — (1 — A) (y/B)] 


where \ denotes the taper ratio of the wing and ¢ is 
the distance from center of pressure to the flexural 
axis in terms of the streamwise chord. The function 
describing the spanwise loading //q is given by Eq. (7). 

The equation for the torsional moment at the root 
may be written as 


B] y 
s*&cg cos A / E — (1 — A): | dy 
0 @g B 
8] 1 BY] 
s*tco cos A | f dy — —- (1 — A) [ y ay| 
0 q B #¢ ¢ 


As will be shown later, it is convenient to derive the 
ratio of the torsional moment at the root to the bending 
The bending moment at the root 


moment at the root. 
is given as 


The equation for the angle of attack at maximum load 


pz { S 


7 4m?F, l2s2 


32 m’* 


ie * —— = 
7 383a (m + 


Qn 


M, s? ~B l 
p - y dy 
q BeosA Jo gq : 
and 
T,/M, = §co[(8/Y) — (1 — d)] cos? A 
where 


i “B / wi 
Y= / y dy, / dy 
ws 2 qd i wet ¢ 
and corresponds to the spanwise center of pressure {or 
the load on the half wing. The value of Y may be 
determined by a graphical or analytical integration ¢ | 
Eq. (7). 
When the assumption is made that the twist due ty | 
torsion varies linearly across the span (or that the rati 
7T/GJ is constant across the span), the incremental angle 


of attack of any section of the wing due to torsional de. | 
flection may be written as 


cos A T'S’ Ts y 
8 GJ 


( 


and by adding this expression to the angle of twist due 
to bending [Eq. (1)], the total angle of twist of any 


Ms y Li Ts 
— —— — | tan A — | — (14a 
EI B J,G/ M, 
ay Eek 1, : 
cs tan A — ( ) | (14b 
An B Jas M, 


Expressions for the 
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d, E 
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Combined Bending and Torsion. 


" Ee 


ie) | 


aerodynamic properties of swept wings experiencing | 


both bending and torsional deformation may be ob- | 


tained from Eqs. (6) to (12) if tan A is replaced by 


[tan A — (1,E/J,G) (T,/M,)| 


factor for combined bending and torsion is, then, 


Ly 7 
| tan ee ( :) | ryt 
J,G/ M, f 
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(1) where 
T,/M, = tceol(8/¥) — (1 — d)] cos? A 
In applying the foregoing analysis to a specific wing, it is convenient to use the equations to obtain the ratio of 
Cn, OF Ci, for the elastic wing to the value for the rigid wing. Multiplying this ratio by the value of C,, ot Cm, 
for the rigid wing, as determined by the complete theory wherein the region within the Mach cone of the tip, ete., 
‘s considered, will then give more accurate parameters for the elastic wing. Then, 
] Am?s . 32 m’* Omar. 8? fol y) | (2) | 
= ily) — : ; 5 tan A — (16) 
qa B= ~~ 3B°’x (m+ 1)? E d, a J.G/ M, 
an 8 2TVm omer. S I,E\ T, | Fe s 
Ceo. . —_ 3 Br(m + 1)? E day | tan ane (Fc) z| F, aa) 
Grigid 
(13 rn 8 TVm omar. S$ .E\ T, | F, 
C, |  3px(m +1)? E dian tan _— a 4 F; us) 
rigid 
In using the preceding equations, it is necessary to solve for a». This, in turn, involves finding the ratio 7,/ \/,, 
which is determined by the parameter Y (usually has a value of about 0.40). 
ire for 
ay be A solution of the combined bending and torsional It is known from recent work by Brown and Adams? 
on deformation effects can be obtained by assuming a_ that a form of strip theory similar to that given by Eq. 
value of Y, solving for a,, and checking the value of (19) gives the correct load distribution for wings with 
ue to | Y from a moment and area integration of a plot of asymmetric linear twist. In the present case, sym- 
ratio | Eq. (16) to see whether a second approximation is re- metric linear twist is involved, but, since the mutual 
angl quired to determine a, more accurately. interference of the two wing panels is small, it appears 
al de The previous equations apply primarily to flat lifting 
wings or to twisted and cambered wings for which the —_ en 
: loading due to twist and camber is the same essentially . | | 
as the loading due to change in angle of attack. For a ~# a = 
wings with somewhat arbitrary camber and/or twist, — == 2 — a 
these equations apply to all accelerated flight condi- — | 
tions. A solution for the aeroelastic characteristics O 20 40 60 80 100 
in steady level flight for such wings must involve a Percent semispan 
consideration of the effects of the loading due to the (a) Twist distribution due to deformation 
- due | known arbitrary twist. 24 = so 
any Strip Theory Strip theory 
It is of interest to compare the results of the fore- . 
14a) | going analysis, which utilizes lifting surface theory for — Lifting surface theory 
the aerodynamic load calculations, with analytical re- 8 16 
sults that more carefully consider the structural char- = 
: acteristics of the wing. Such a comparison involves the s 
4p) | development of a strip theory of good accuracy which 3 12 
may be used in matrix calculations. One form of strip d 
the theory may be derived by using the basic loading for aS ia 
“ing | a flat swept wing at an angle of attack. The load 8 
= at any spanwise station produced by an angle-of-attack 
change for that station is assumed equivalent to the O4 
load produced at that station by the same angle-of- 
attack change for the entire wing. The expression for 
the load may be written as 
‘ 0 20 40 60 80 100 
l/qa = (4m*s/BY)fi(y) [1 + (az/a)] (19) Percent semispan 
Comparisons of the load so computed with lifting-sur- (b) Load distributions 
15) face theory show that this form of strip theory is not of Fic. 5. Comparison of the load distributions as given by 
strip theory and lifting-surface theory for a wing with parabolic 





i sufficient accuracy, deflection curve. 
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that the strip theory resulting from Brown and Adams’ -6 relatet 
work may prove satisfactory for aeroelastic calcula- a beat 
il “ae: ; a . 
tions.* The expression for the load at any spanwise re P 
ot. bo 1c st72 ae —_ — 
station 1s given as Ln 
1 4m?*s > ar (a) 20 40 60 80 700 The § 
-—— feet + — > = Percent semispan ; swept 
qa B= _ Sed... m= oe T} 
—- : (a) Twist distribution due to deformation Phe 
b er Ll — 9" , 
2. In pal 
(2( ) Rigid wii the el: 
A comparison of the loading given by Eq. (20) for a — 
, . ‘ ; ; : yal 
wing with a parabolic deflection curve with the loading Lifting surface vathal 
‘ a P functs 
for the same wing derived from lifting-surface theory is & oa i —_ 
: ‘sme. fh r P a 5 i) Or, son 1S 
given in Fig. 5 fora Mach Number of 1.414. The span- ~e ? ‘4 “ 
= . . ‘ - ‘ ‘ valcu 
wise twist distribution corresponding to the parabolic x ; 
: : a as (1) 
deflection curve is shown in part (a) of the figure. The $ 
, : : oa ~~ wing 
loading for a flat wing is also shown. The same com- Ss 8 

‘ ° : x: , : . owe! 

parison is made in Fig. 6 for a cubic deflection curve. $ I : 
. ‘ é (2) 

In both cases, the strip theory gives a close approxima- 2 

‘ ar . dai : d assull 
tion ‘to lifting-surface theory. Since the deflection 8 3) 
. : : ; (: 

curve of most wings can be closely approximated by baa 
A a d igi : : to be: 
suitable superposition of segments of parabolic and 04 Fo: 

. . . ° . . , 40 
cubic arcs, it appears that this form of strip theory : il th 

© — a 
should be of sufficient accuracy. nate 
‘ : +e ° , , ‘ sign | 
If this simplified form of aerodynamic loading is = 
as m - : : Sih rs) 20 40 60 80 100 siona! 
utilized, Eqs. (15) to (18), which are based on lifting- 

‘ ; Percent semispan show 
surface theory, may be rewritten as sa nee 
- (6) Load distributions carve 

* The work of Goodman of Cornell cited in < “vious footnote . ae F ae ee . 
aencap tet te : pn diaageaneae omg eet a i Fic. 6. Comparison of the load distributions as given by strij up 
can be used to develop a similar strip theory that the authors have theory and lifting-surface theory for a wing with cubic deflection amou 
found possesses equalaccuracy. — curve. Wi 
form 
are 0 
> effect 
l 4m?sfi(y) a Ley T, sai 
= — i-——= Pits = mm _ m* : tan A — ‘ : (21 as th 
qa Bz Bd,Ea, 2 = J,.G/ M, there 
= ye Ll = 
assur’ 
- 
B?z oe x. 8m? $ Omar - = Lact i, . ~ 
% = — te. + —— : 2— m? _ m?* , tan A — : - iF, (22 
4m?F, | 2s? Bd, E -=— + J,G/ M, 
L = oe 1 — m* 
-~ 
C ae z — 
Vactastic 2 Ss Oo mar. 9 z ° I,E 7 r F 5 99 
; =l1-— : a= Mm m* ; tan A — : : (23 
... B d, Fay, -z=- = J,G/ M,\ F; 
“rigid | le | = ms 
C 2 
W tetanic 2 S Omar. 9 = - LE T; 7 F, 9 , 
: =1- 2—m*_ m? . tan A — : : (24 y 
"9 ae B d, An me z- v I JG M, F; 
rigid 1 — m* 1 — m* 
The constants F; and Fare given in the Appendix. 
/ 
: / 
Matrix Calculations these calculations, it was assumed that the section ’ 
/ 


F : ae ; . moment of inertia of the wing beam varied with the 
In order to investigate the validity of assuming that. f : i : 
’ : ‘ fourth power of the spanwise dimension so that 
the wing deflection curve is parabolic, calculations were 
eqe . . ° = om — ’ ;) 4 

made utilizing strip theory with the matrix methods of I i, (1 d) (nn/s)] 


a recent paper by Diederich,* of the N.A.C.A. In and that the torsional stiffness constant was directly 
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related to the bending moment of inertia of the wing 


beam as 


J = (3/2)1 


The geometric and structural characteristics of the 
swept wing of Fig. 7 were used in the analysis. 

The results of the calculations are shown in Fig. 8. 
In part (a) of Fig. 8, the ratio of the lift-curve slope of 
the elastic wing to that of a rigid wing of the same plan 
form is shown as a function of the dynamic pressure. 
In part (b) of Fig. 8, the shift in the neutral point as a 
function of the dynamic pressure is given. A compari- 
son is made between the aeroelastic characteristics as 
calculated by: 

(1) Strip theory with matrix calculations for the 
wing with the moment of inertia varying as the fourth 
power of the spanwise dimension. 

(2) Lifting-surface theory with the deflection curve 
assumed to be of parabolic shape. 

(3) Strip theory with the deflection curve assumed 
to be of parabolic shape. 

For the comparison of this figure, the calculations for 
all three cases were based on the same maximum de- 
sign stress at the root section, and the effects of tor- 
sional deformation were not included. The results 
show that the assumption of a parabolic deflection 
curve gives reasonable results for all dynamic pressures 
up to a value for which the loss in lift-curve slope 
amounts to 40 per cent. 

When the effects of both bending and torsional de- 
formation are considered, the results shown in Fig. 9 
are obtained. Matrix calculations that treat torsional 
effects more accurately give essentially the same result 
as the analysis of the present paper. It is evident, 
therefore, that, for the types of wings considered, the 
assumption of a parabolic deflection curve gives a rea- 








4 % A.R.= 3.28 ri \ 
E=/05 x 10® 4,78 = 0.052 
G=385 x 10° A= 0.25 

IJ = 0.67 f=.2 


Sketch of the wing used in the calculations. 
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(6) Neutral point shift 


Fic. 8. Variation of lift-curve slope and neutral point shift 
with dynamic pressure at \J = 1.414 for the elastic wing of Fig. 
7, experiencing pure bending deformations. 
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(b) Neutral point shift 


Fic. 9. Variation of lift-curve slope and neutral point shift 
with dynamic pressure at JJ = 1.414 for the elastic wing of Fig. 
7, experiencing combined bending and torsional deformations. 


sonable estimate of the longitudinal stability charac- 


teristics. 


Mach Number Effects 


It is of interest to compare the results obtained for 
supersonic speed for the wing of Fig. 7 with calculations 
for the same wing in incompressible flow utilizing the 
aerodynamic loadings given by the Weissinger method 


as developed by Stevens.’ A comparison is shown in 
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Appendix 
MATHEMATICAL DERIVATION OF LOADING FUNCTIONS AND PLAN-FORM CONSTANTS 
The functions f(y), fo(y), f(y), and fs(y) and the constants F,, F2, F3, Fs, Fs, and Fs, which appear in Eqs. (7) to 
(23) of the text, are given in this Appendix. The evaluation of the constants from the analytical expression may 
prove to be tedious, and it may be easier to evaluate the corresponding integral equations graphically. 
The functions f;(y) and f2(y) were developed from the following integral: 
(y + mtco)/mt Ap | 
l/q=s dx 
y/m q fy const 
from which 
; (y + mtco)/me dx 
hi(y) = 7 
y/m V m* — (y/x)? 
which yields 
fily) = (1/m?*) V [(m?/m?) — 1] 9? + (2m?/m,)coy + mc? 
and 
' 
(vy + mtco)/mt x 4 y (v + mtico)/mt die ae y ' wher 
fro(y) = (x + y) dx + (x — y) dx (m, - 
J y/m mae ¥ / y/m mx -- ¥ 
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which yields 


2m(y + mo) + my(5m + 3) 
f(y) = 4m*m 


3y*(m + 1)? 2m(y + mo) — my(1 — m) 2m(y + mo) — miy(om + 3) 
| the : [eost : . — : + : x 


Sm m(m-+ 1)y 4m?m 


a 
Ib (y + meco + my) (my + mmo — my) | + 


3y?(m + 1)? 2m(y + meo) + my(1 — m) 
= cosh~! = aa 
m(m + l)y 


Ee + 1) V2(1 — m) ke > i? | 3 | 
i _ : , cosh~! 
4m”? Sm’? m+ 1 


B 
The constants F, and F are evaluated as follows: Fy = fi fily) dy, which yields, for m, ¥ 1, 
0 


[~ (y + mio — my) (my + mmo + my) | = | 


8m’? 


; ap — mf ,— mmc” — mz,f? , m*f — ap _,m 
Fi = — V mf? + 2pm*f — af? + > 5, */2 (COS — cos 
2m ma ; 2a ™~* mm, f mM, 
and, for m, = 1, 
_ Bl — m?) — mo ——= - = - 
Fi = 7 = V mco? + 28m2%co — B2(1 — m2) + 
2m*(1 — m*) 


mC” Co" m*co — B(1 — m?) 
+ cos”! — cos~'m 


2(1 — m?) 2(1 — m7?) mCo 


*B 


where a = m,?> — m?, f = mo, and Fy, = J fo(y) dy, which yields, for m, ¥ 1, 
e { 


) 


e Bes : f(de + 4mb) 
F, = ( ) [(mf?) ’? — (Bdf + mf? — B*b)’] + 1° x 


12m?m 7b, 8m?m 7b 


j28b —df , if 5 + m “df? L m(m + al _ af — 2pb ' d 1 
VB 2— §? - : cos : — Cos 
| 1b pay - my ’ Ab 85” mf(m + 1) mim + 1) Jf 
B3(m + 1)? Bd + 2mf 38h + 2mf 3—m f(m 1)? 
| : [eos l — + cosh! : “— — cosh! a T x 
Sm Bm(m + 1) Bm(m + 1) m-+ | Sm-* 
jf?(3d? + 4mb) | , af — 2b . d | 3df + 2Bb oy er ; ; 3m ‘“df*\ 
‘OS — cos — ’ 2 r __ pe : 
{ Sb mf(m + 1) mm + 1)] 4b? my" + Bas B*) + 4b? f 
B3(5m + 1) V2(1 — m) | (mf? + Bhf — Bg) ’* — (mf?)” f(A4mg — hj) 
=m], (4) se" — cor) [sme = 1) 
12m “’ 4m?m ,? 3g Sm*m (°¢g 
j262 —hf ;,— 4 m'*f2h n f?m ,?2(1 + “| ’ hf — 28g ; h ) 4 
V - Bhf — Bg ; : = cos > — cos a 
{ tg my" + PRY £ tg Sg mf(m + 1) m(m + 1) Jf 
fim + 1)? ] (PP Bh? + 4mg) hf — 26g h 28g + 3hf 
= +: : cos : — cos! —- —- X 
Sm? ( di mif(m + 1) m(m + 1) 1g? 
[= : 3m hf?) 
V mf? + Bhf — B’g + 19? 
a 
where b = (m, — m) (m, + 1), d = 2m + mm, — m, e = 2m + 5mm, + 3m, f = me, g = (m+m,) X 
(m,—1),h = 2m — mm, + my, and j = 3m, + 5mm, — 2m. Form, = 1, 
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2m l6m*(m — 


| 6(m — 1) 6(m — 1)9 


j| 480" — 1) + ep(3m — 1) Co(38m — 1) 
) [Vv mco? + Boo(3m — 1) + 282(m — og 8 vx 
(m — 


S(m — 1) 


(v ) | co2(m + 1) “all ,48(m — 1) + co(3m — 1) _3m — | a r+ (* (m + 1)2 
/ -2 cos — COs 
wed 1677 2(1 — m)” Colm + 1) m+ | Sm » 


B(3m — 1) + 2mco Bim + 1) + 2mco 3—-—m Co°(m + 1)2 
cosh~! -+- 


F (@ + *) f [mco? + Beo(3m — 1) + 282%(m — 1)]*” (mco*) ”* \ E Co(3m — 1) (7m + ra 
= ; x 
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x< 


+ cosh! — cosh 


B(m + 1) B(m + 1) m+ | 
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) 0 — ) _— 
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> 4 
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Co 


327 2m(1 — m) ‘o(m + 1) m+ | 


j2[Sm?co" — 4Bmce(m + 1) + 3887(m + 1)7]) 16m “°c 


L5mcy *(m + 1)8 
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Ee +1) V2(1 — m) S(mceo?) * ) | 
; — 2 
12m lom(m + 1)] 


The functions f;(y) and f;(y) were developed from the following integral: 


(LV mco + B(m + 1) | i 


L5(m 


; 
; f 
msg ROS ee 
q 7 y/m q y = const 
from which 
Pin F, 
; Vv + mMico)/mt x dx 
daly) = 
Jf y/m Vm? — (y x)? 
- mM, 
Then, 
y+ OF / y? m(y + mo) 
(vy) = < / , a 242 + — -osh—! - 
Fs(y) a a m*(y + mo)? — m4 —— onl 
2m*m ,” 2m' my 
Also, 
™(y + mico)/mt x + y ~(y + mtco)/mt t—y 
filly) = (x + y) — xdx+ (x — y) — x dx 
y/m mis— ¥ y/m mx + y 
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[Vv 
fily) = LV (vy + mo + my) (my + mmo — my) | xX 
8(yv + mo)? 4 2mm y(y + mieo) (7m + 5) + m,y?(8m? + 22m + 15) r pa — 3m> — 9m — ”) X 
24mm ,? 24mm ,? 16m 
(2m — mm, + m,)y + 2mm co 3—m (2m + mm, — m,)y + 2mm co 
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| ~B 
"| x The constant F; is evaluated as Ff; = / fs(y)dy, which yields, for m, ¥ 1, 
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) . ; « . - . / or 
X mf? — (m2f? + 2Bm?f — Ba)” B* mB+f) flap — 3m*f) V mf? + 2Bm'*f — ap? 
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A System for the Excitation of Pure Natural 
Modes of Complex Structure 


ROBERT C. LEWIS* ann DONALD L. WRISLEY?# 
Massachusetts Institute of Technology 


SUMMARY 
Analogy is drawn between analytical and experimental ap- 
proaches to the evaluation of structural dynamic behavior, 
pointing out the relationship between method complexity and 


adequacy of results 

The theory, construction, and operation of a system using 24 
shakers controlled from a console in the systematic excitation of 
pure natural modes are described. Practical operation is made 
feasible by the coordinated display of the force and resultant ve- 
locity at each shaker station, together with a system for orderly 
discriminatory adjustment by the operator, which is analogous 
to a mathematical iteration. Sample laboratory and full-scale 


data are shown. 


INTRODUCTION 


— OF THE DYNAMIC BEHAVIOR of complex 
structure has become increasingly important in 
the design of modern aircraft whose dynamic loadings 
influence both aerodynamic and strength considera- 
tions.' The adequacy of these calculations is dependent 
on the thoroughness of the mathematical treatment and 
on the correctness of the engineering evaluation of 
boundary conditions required by physical circum- 
stances. 

Development of mathematical techniques and _ in- 
genuity in their application are continually being em- 
ployed to extend analytical accomplishment. Since it is 
by experimental verification of the analytical results 
that the adequacy of the engineering evaluations are 
proved or found wanting, a similar degree of develop- 
ment and ingenuity is required in the experimental 
evaluation of these same complex structures. 

Experimental determination of the natural mode 
responses of the complex aircraft structure is required 
by various contracting agencies.” Unfortunately, the 
results of many shake tests are in a form difficult to 
compare with the calculated data, and a large part of 
their value is thereby lost. 

The lack of clarity in comparing calculated and ex- 
perimental data on natural mode responses does little 
to aid in checking the validity of the theoretical as- 
sumptions; it may help to confirm the analyst to more 
exclusive dependency on mathematics as the ultimate 
definition of structural behavior, or it may further con- 


Presented at the Aeroelasticity Session, Eighteenth Annual 


Meeting, I.A.S., New York, January 23-26, 1950. 
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vince the experimentalist of the perversity of nature 
oversimplified by the theorist. Neither of these re- 
sults is desirable. 

It would seem worth while to examine the possible 
reasons for the confused picture of the natural mode 
responses obtained from some shake-test data: 

(1) The shake test on a new airplane is usually made 
on the prototype prior to flight. Schedules may be 
lagging and pressure high. Even in the presence of a 
contract requirement, the shake test smacks of flogging 
a dead horse to the project nearing the already long- 
delayed fruition of flight. 

(2) Budget and time may not allow check of data 
for adequacy during the tests. They almost certainly 
prohibit the rerunning of any data found inadequate 
later because of the extensive set-up time required. 

(3) The experimental procedures used may in some 
cases be naive by comparison with the mathematical 
methods whose results they are expected to check. 

There is certainly more than a superficial analogy 
between the mathematical analysis and the experi- 
mental analysis, and the results in either case will be 
no more elegant than-the methods devised to treat the 
problem. The inability in some cases to obtain clear 
definition of the mode shapes experimentally in com- 
plex structure, by the use of techniques adequate for 
one or two degrees of freedom, should not be surprising. 

The instruments available to the experimentalist 
have developed in a manner similar to the methods of 
the analyst, and there is a far greater profusion of good 
instruments here in the United States than elsewhere. 
It would seem worth while to utilize this material ad- 
vantage by the continued exercise of ingenuity in the 
development of experimental techniques capable of 
evaluating the physical complexities described by the 
language of the theorist. 

An attempt along this line is embodied in the equip- 
ment and method of use developed under Government 
contract{ and described in this paper. 

THE GENERAL PROBLEM 

Aircraft structures are characterized by distributed 

flexibility and mass and a general complexity of struc- 


tural and spatial arrangement that inany 
These modes may have frequencies 


leads to 


natural modes. 
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close together and may be complicated in mode shape, 
with all parts of the ship more or less active. 

This general complexity and activity is, at one and the 
same time, (a) the major reason for making experi- 
mental check on a somewhat simplified analysis and 
(b) the major difficulty in obtaining clear experimental 
definitions of the actual responses. 

In general, the excitation of the structure is applied 
at arbitrary points, and the resultant motion of the 
structure at the forcing frequency may be thought of as 
the summation of the forced responses of each of the 
structure’s many natural modes. Since some of the 
modes have natural frequencies close together, each of 
several modes may be excited to substantial amplitude 
by an arbitrary excitation, and various parts of the 
structure, each perhaps active in a different mode, will 
exhibit substantial phase differences. If under these 
conditions the alternating forces are removed, each of 
the several modes decays at its own frequency, with 
the production of beats. 

The objective sought for the experimental method is 
to produce, in the complex structure, by means of ad- 
justable forces, oscillations that consist essentially of 
one natural mode of motion. For the purposes of ex- 
plaining the basic concepts of the development, it is 
convenient to assume that such a pure natural mode 
is by some means sustained. The properties of natural 
modes are such that all points in the structure will 
reach their maximum displacements at the same in- 
stant. At the instant when the structure reaches the 
maximum of its complex deflection form, we can 
visualize this elastic deformation of the structure as 
the result of a specific distributed inertial loading, and 
we can further see that the distributed inertial loading 
produces shears and moments on the structure which 
are exactly balanced by the shears and moments in the 
elastically deformed structural members. 

In the absence of damping forces, such a natural 
mode would continue in the steady-state oscillation we 
have imagined without the application of external 
forces. If, on the other hand, damping shears and 
moments did exist in the structure, we could still pre- 
serve the steady-state natural mode if we exactly canceled 
these shears and moments by externally applied oscillatory 
forces. 

This is the basic concept around which the method 
was formed. This assumption of the accomplished 
existence of a pure natural mode allows certain infer- 
ences of importance in devising the experimental 
method : 

(1) Since the natural mode motions are all in phase, 
the damping forces and, therefore, the driving forces 
will all be in phase when the final desired condition of a 
pure mode is obtained. This means that phase manip- 
ulation of the driving forces will be not only unnecessary 
but undesirable. 

(2) If our externally applied forces are not properly 
adjusted to cancel exactly the internal damping shears 
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: 
and moments, the resultant motions of the structy, 
will depart from the “in-phase’’ condition. This jg 

valuable criterion for mode purity but not a suitaty 
criterion for corrective adjustment, at least for J 
present. ' 

(3) If the motion existing under the action of ¢) 
externally applied oscillating forces is a satisfacto, 
approximation to a pure natural mode, a smooth deca 
of the mode should be obtained on removal of i, 
excitation; freedom from beats will indicate at least, 
practical approximation to the required pure mode g 
citation and will indicate conclusively that the » 
perimentally determined frequency, mode shape, apj 
damping coefficient describe a true natural mode, 

(4) 
distributed points on the complex structure in ord 
that a unique adjustment of the magnitude and sen 
of the shaker forces will be obtained for each natury) 


A large number of shakers will be required y 


mode shape.* 

The major fundamental problem in applying thi 
simplified concept concerns the mechanics by which 
external forces can be adjusted to produce shears an( 
moments to cancel exactly those produced by damping 
within the structure. We are brought face to face with | 
the need for knowing more of the specific mechanisns | 
of structural damping and its distribution. As a rea. | 
sonable but simple assumption, we can take the concept 
of structural damping literally and assume that each 
elastic shear and moment in the structure has asso- 
ciated with it a proportional damping shear and mo- 
ment in quadrature phase relationship. 


Since shears and moments from the inertial loading | 


exactly cancel the elastic shears and moments, it can | 


be seen that, because of the assumption of damping | 
distribution, the externally applied forces required to | 
cancel the damping shears and moments must be dis- 
tributed over the structure proportional to and in 
quadrature with the inertial loading. That is, each one 
of the numerous oscillating forces acting on a complex 
structure at a natural frequency must be adjusted in 
magnitude to be proportional to the product of the 
mass of the structure on which it acts and the ampli- 
tude of that mass in the mode being excited. 

Thus, the result of assuming a distribution of damp- 
ing proportional to stiffness provides a practically use- 
ful criterion on which to base discriminatory action. 
The question remains as to how nearly a pure mode 
might be excited by use of the criterion, even in prac- 
tical cases where the simple assumption of damping 
distribution is violated immoderately. 

In this connection it is timely to point out that nor- 
mal mode concepts apply strictly to systems with zero 
damping, or a “small amount of damping,’’ and are 
stretched with difficulty to cover the case of distributed 


* This concept was formulated in memorandum Aer-RD-223- 
LEB, January 28, 1946, in which Lt. L. E. Beckley, U.S.N.R., 
outlined the essential advantages and requirements of such a 


system. 
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SYSTEM FOR EXCITATION 
systems with discontinuous energy sinks or even lumped 
systems with nonuniform damping elements. 
“Unfortunately, the actual structures tested by the 
experimentalist do have concentrated energy sinks, and 
reasonably clear data are still desired for comparison 
with the calculations unaware of such circumstances. 
For example, conventional shaking of a wing at the tips 
may evoke enthusiastic response from the empennage 
at one of the “‘wing’’ natural frequencies. The experi- 
mental data will be confused by large phase shifts be- 
tween empennage and wing motions, lack of clearly 
defined nodes, and beating on decay of the “natural’’ 
mode so excited. This situation may well be caused 
in large part by the phase shifts, due to transmission 
through the flexible wing and fuselage to the tail, of 
the forces introduced at the wing tips. The energy is 
fed in at the wing tips and is dissipated throughout the 
structure and, perhaps, in large sinks at the wing- 
fuselage joints. The dragging along of other modes 
consequent to this procedure might be expected. At 
the same time, much might be done to eliminate this 
damping- or phase-coupling by simply giving the 
empennage and other major parts of the structure 
firsthand knowledge of what is going on by the use of 
other shakers at these locations, all tied to the same 
source. 

A good understanding of the concepts presented here 
in descriptive terms is essential to advanced experi- 
mental work. A complete mathematical treatment of 
the theory of small oscillations of conservative systems 
and of systems with small damping may be found in 
von Karman and Biot,* chapters 5 and 6, and in Ray- 
leigh,* Section 4. A summary of this theory is pre- 
sented in Appendix I of this paper. Kennedy and 
Pancu’ have made a masterful presentation of con- 
cepts pertaining to natural modes as a set of gener- 
alized coordinates and have presented a method com- 
bining experimental and analytical procedures in the 
separation and definition of structural natural modes. 


THE SPECIFIC PROBLEMS 


Evolution of the general concepts outlined into an 
integrated set of equipment for testing full-scale air- 
craft required four major steps. 

(1) Analysis of the theoretical system requirements 
for exciting natural modes, relating the structural 
properties, the shaker properties, the mode amplitudes, 
and the power requirements. 

(2) Construction of experimental equipment, de- 
signed with consideration of the theoretical analysis. 

(3) Experimental development. (a) Investigation 
of suitable criteria for mode purity by means of 
plex model with an adjustable distribution of 
sinks. (b) Determination of suitable means 
criminatory action in adjustment of multiple 
forces to achieve mode purity as indicated by the 
(c) Evolution of a working system, 


a com- 
energy 
of dis- 
shaker 


criteria discovered. 
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based on practical requirements, on which to base the 
design of complete full-scale equipment. 

(4) Construction and test of the full-scale equip- 


ment. 
ANALYSIS OF THEORETICAL SYSTEM 


Power Expressions 


Assume now a structure whose mass distribution is 
represented by the function m = m(x, y) oscillating in 
the ith normal mode. The deflection of any point in 
the structure can be represented by the product of a 
real, normalized function f; = f;(x, y) and a complex 
coefficient do (= do sin wt). The analysis can be ex- 
tended to include the z coordinate also. 

According to the assumption of 
damping forces proportional to restoring forces (and, 


distribution of 


therefore, inertia forces), the absolute value of the 
damping force at any point (x1;) is 
dfa\ = \gidfe) = |gaowfi(xiyi)m(xiy1) (1) 


where df, and dfg are damping and restoring forces 
and g, is the damping coefficient for the 7th mode. 
The power absorbed by the force df, is 


Gf. = [dfa: dow f (x11) | 2 


and the total power absorbed by the structure is 


l 
P, = > gow, f(x, y)m(x, y) dow f(x, y)dx dy 
_ x y 
(3) 
gow; way = 
P, = = m(x, y)fir(x, y)dx dy (4) 
“ FJ tJ 9 


In evaluating the power supplied by the shaker forces, 
it is necessary to recognize that the physical limits of 
the actual shaker equipment must not be exceeded. 
It is desirable, then, to define the following quantities: 

i, = maximum force rating of an individual shaker 
M = the value of the function m = m/(x, y), where 

the product m(x, y) fi(x, y) is a maximum 

for the 7th mode 
the value of the normalized deflection function 

fi = f(x, y) corresponding to M 

force at any point (x, y) 


tu _ 


i= 

It was stated previously that the magnitude of the 
applied shaker force should be proportional to the 
product of mass and amplitude at the point of applica- 
tion 


1 = cmp; (5) 


and the maximum force applied to the structure will be 
equal to the maximum shaker force 
1; = cMfu (6) 


from which the proportionality factor, c, can be evalu- 
ated for the ith mode under the condition that the 
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Fig. 1 Pure Bending 
m, = total mass of the nt) section 
x, — distance from the center of gravity to the elastic axis 
= radius of gyration of section n about its center of gravity 
fin ~ Value of the normalized deflection function f,,(y) at station n, 


fiqj ~ Value of the normalized deflection function at the station where the 


product - fin iS @ maximum, 
4, ~ circular frequency of the bending mode in radians/sec. 
8, ~ damping coefficient for this bending mode. 
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Fig. 2 Pure Torsion 
m, = total mass of section n 
x distance from the center of gravity to the elastic axis at station n 


/, ~ radius of gyration about the center of gravity at station n. 


fin the value of the normalized torsional deflection function at station n 


fen the value of the normalized torsional deflection function at the station 


fr 
where the product —2 +f i 
produ - (x, n?) fey 1s @ maximum. 


®, = circular frequency of the torsion mode 


cA damping coefficient for this mode. 
. m, 

M value of _ corresponding to fey 
1 maximum shaker force rating 





shaker called upon for the largest force is adjusted to 
its maximum rated output. 


c= (1, ma tar) (7) 
the force at any point 1s now 
t= (i,/M)(fi/fu)m (S) 


The power delivered to the structure by the force 7 is 


| (‘:) ff ; 
=} — m-aqw,f; 
2\M tu ; 


and the total power delivered by the distributed forces 


dP, = (9) 


1S 
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AoW; 


t p ‘ 9 
2 *) far / J M(x, y) fP(x, y) dx dy (1 


Amplitude at Resonance 


Py = 


Equating power delivered to power absorbed, thy 
magnitude of the complex coefficient a, which evaluates 
the deflection, is obtained. 


l (> l 
a= -,> = (11 
, wg; \M tu 


The amplitude that the structure attains (subjec 
to shaker stroke limitations) is seen to be dependent oy 


frequency, w;; damping coefficient, g;; and the quan. | 


tity (4; M)(1/fw), which includes mass distribution 
mode shape, and the force limitation of the shaker. 


Power Requirements 


The power required to excite a structure to an am | 


plitude dp in a normal mode is, from Eq. (4), 


a “s ad 
g 9°; »e 
- m(x, y) fir(x, y) dx dy 
= x vy 


and, substituting the value for a from Eq. (11), 


' i\2 ae 
P= ( +) —" | / m(x, v) f(x, y) dx dy 
2u ifr M B; uo" x y ; 


(12 


ri * 


While normal modes may be, and usually are, a com- 
bination of bending, torsion, and possibly chordwise 


motions, they may often be considered predominately | 


“‘pure’’ bending or torsion or may be resolved to pure 


bending and pure torsion components. It is useful 


— 


therefore, to write Eq. (12) in terms of pure bending and | 


pure torsion. 


Since in an actual test program a finite number of | 


shakers is available, the expressions for pure bending | 


and pure torsion will be written in terms of summa- 
tions. 

The structure can be approximated by a dynamic 
system where the mass distribution is replaced by 
lumped masses at appropriate distances from the 
elastic axis and the axis of roll to approximate the 
actual inertia properties. 

For the case of pure bending (i.e., no twist about 
the elastic axis, see Fig. 1), we can represent the dis- 
placement of the elastic axis by 


Sp = fan(y) (13) 


where f,(y) is the normalized deflection curve and dp 1s 
a complex coefficient. 

The power required to excite the structure to an 
amplitude do in the pure bending mode is 
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do = (1/wr7gn) (is/M) (1/foas) (15) 


For pure torsion, there is no displacement of the 
elastic axis (see Fig. 2), and we can express the twist at 
any spanwise station by the normalized function f,(y) 
and complex coefficient ag: 


(16) 


a= ao fily ) 


and the displacement of any chordwise point due to 
twist, assuming no chordwise bending, is 


x aof (y) (17) 


s, = 


f 


Now, the power required to excite the structure to 
an amplitude apo in the pure torsion mode is 


I (; y 1 > - 
i 2w 2; M fis? (Eu + pw)? =1 
My (Fn? + prn?| fin? (18) 
and 


I 1, | 
i és _ - (19) 
wee \M/ feu (Ear + par) 
Assuming no chordwise bending, any normal mode 
motion of the wing can be resolved into two compo- 


nents—pure translation of the chord line and twist about 


the elastic axis. 
ponent may be found by using the above expressions, 


The power required for each com- 
and the total power is their sum. 


CONSTRUCTION OF EXPERIMENTAL EQUIPMENT 


It was desired to have the experimental program serve 
in some measure as a proving ground for the equip- 
ment later to be procured for delivery to the contract- 
ing agency, and it was also desired to have as much as 
possible of the equipment either commercially avail- 
able or conducive to the use of standard switch gear. 

Electrodynamic shakers were the obviously prac- 
tical means of delivering forces to the structure, variable 
in both frequency and magnitude and fixed in phase 
but reversible in sense. 

Selection of the required force rating of the shakers, 
as well as the power-handling capacity of the variable 
frequency source, the distribution system, and the 
control system, was based on application of the ampli- 
tude and power relationships derived in the previous 
section. Typical mode shapes for aircraft structures 
of known mass distributions were analyzed, and it was 
concluded that the final system should consist of a 
variable frequency alternator having a capacity of 2.5 
kva. at 100 cycles per sec. powering 24 shakers, each 
capable of delivering an alternating force of 25 Ibs. 
through a stroke of +!/.in. 


Alternator 


An alternator was required to furnish good wave 
lorm at as low as SO per cent power factor over a fre- 
quency range of from 4 cycles per sec. or lower to 100 
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cycles per sec. or higher. Although the standard 115- 
volt output at 60 cycles per sec. favored the use of stand- 
ard switch gear, it seemed desirable to have a special 
alternator built rather than use a standard 60 cycles 
per sec. alternator. The specific advantages of the 
special alternator obtained are (1) reduction in har- 
monic content of the current output; (2) provision of 
forced cooling for low speed, high current operation; 
(3) specific design for rotational speeds required to pro- 
duce 100 cycles per sec. output. 

A consideration of the overall system requirements, 
discussed later, determined the required alternator 


rating. 


Power Distribution 

The power distribution system for the experimental 
work was essentially the same as that described later 
in connection with the final equipment. Essential ele- 
ments in the control circuit were (1) a large Variac 
on the alternator output to control the supply to a 
feeder line and (2) smaller Variacs to control each shaker 
individually. 

The current to each shaker from its Variac (and, 
hence, the force delivered) was reversible in direction 
by means of a standard switch. 


Shaker Design Considerations 


Elementary texts in either physics or electrical 
engineering show that, when a current is caused to flow 
in a conductor placed in a magnetic field, a force is 
exerted on the conductor perpendicular to the con- 
ductor and to the lines of magnetic flux. This relation- 
ship is expressed mathematically when reduced to 
pound-inch-second units by 


t = 0.885 X 10-"BLIN (20) 
where 
it = lbs. force 
B = flux density = ¢/A, lines per sq.in. 
L = conductor length in flux field, in. 
I = current, amp. 
N = number of turns of length ZL in flux field 


If we divide through by J, we have the relationship, 


i/IT = 0.885 XK 10-7 BLN (lbs. per amp.) (21) 


For any particular shaker, B is a function of the field 
structure design and L and N are functions of the arma- 
ture coil design. 

Therefore, the force/current ratio is a constant for 
any given shaker (assuming constant field strength), 
which indicates its performance as an electromechani- 
cal transducer, in the conversion of current to force. 

If we had only to produce large forces, it would be 
comparatively easy to do so by using an armature coil 
of many turns of small-diameter wire. However, it 
must be remembered that we wish to produce a force 


to act on a moving or resonating structure. In other 
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words, it is not sufficient to produce a high ‘blocked 
armature’ force, but we must so design the system that 
the force can be delivered for any expected velocity of 
the moving system—we must be able to deliver power 
to the mechanical system. 

The second basic electromagnetic equation shows 
why it is difficult to maintain the current (and, hence, 
the force) when the armature is moving. This equa- 
tion states that, when a conductor moves in a mag- 
netic field, a voltage is generated. This voltage is such 
as to cause a current to flow in a closed conductor which 
will produce a force opposing the motion. For this 
reason, the generated voltage is known as a back volt- 
age, and its value is given by the equation 


Ep = 19 BLNe (22) 
or, dividing through by e, 
E,/e = 10°* BLN (volts/in. per sec.) — (23) 


where /, is the back voltage (volts) and e is the velocity 
of the moving conductors relative to the field (in. per 
sec.). 

If we now solve for BLN from Eq. (21) and substi- 
tute in Eq. (23), we find 


Ep/e = 0.1125 1/] (24) 


’ That is, the higher a force/current ratio we design 
for, the higher a back voltage will be obtained for a 
given velocity. 

By referring to Fig. 3a, it is seen that, when an open 
circuited coil is moved in a field, the open circuit volt- 
age is equal to the back voltage. 

On referring to Fig. 3b, it is seen that, in order to 
drive a current through the coil in a direction producing 
a force in phase with the velocity, it is necessary to 
have a supply voltage equal to the voltage drop through 
the coil electrical impedance plus the back voltage. 
That is, 


Es = IRe +- Ep 


(25) 


Thus, it is seen that provision of a high force /current 
ratio is not the complete answer to a suitable shaker 
design, because the higher this ratio, the more voltage is 
required to get current through the coil. 


SCIENCES—NOVEMBER, 1950 

The proper approach to the design problem jg 
through consideration of the mechanical mobility of the 
systems that will be driven and the matching of this 
mechanical load by a suitable reflected electrical jp. 
pedance to a practicable power supply. 

The impedance-mobility relationship is found py 
dividing Eq. (22) by the relationship J = 7/(i/T) after 
suitable transposition, giving 


Es I = 0.1125(4 T)e 1 (96) 
or 


Rp = 0.1125(4/1)% 


(26a 


where F,// is the reflected impedance (resistive at 
resonance) and e/7 is the ratio of the velocity of the 
driven system to the force applied or the mechanical 
mobility of the system being shaken. 

When the two latter quantities are evaluated, then 
the required 7/J is established. Following Eq. (21), it 
was remarked that this factor is controlled by B, or 
field coil design, and the quantities L and N that belong 
to the armature coil design. Actually, the two designs 
are closely interrelated. 

The practical limitation of current for any given 
armature coil is dependent on the /°R power losses 
that can be dissipated without exceeding commercially 
acceptable temperature rises. A _ similar limitation 
applies to the design of the field coil. The magnetic 
flux obtainable with a given structure design is also 
dependent on the material properties affecting magnetic 
reluctance. 

A further design consideration was necessary in at- 
tempting to provide the maximum force output possible 
with a shaker still light enough to be considered port- 
able or movable by a single man. A shaker designed to 
meet these requirements was available, and, after due 
consideration of the mobility relationships, an arma- 
ture coil having an 7/J of approximately 3.5 was de- 
signed to provide a suitable match for the power supply. 


Shaker System Characteristics 


It is of paramount interest to know the characteristics 
of the overall power system in order to understand the 
source of various limitations. 

For a given setting of the field excitation, the alter- 
nator terminal voltage at no load varies linearly with 
the alternator frequency, and the voltage at any par- 
ticular frequency can be varied by a change in field 
excitation. It is desirable to have as high a voltage as 
possible at low frequencies, as will be seen later, and 
yet the output voltage at any frequency and light load 
must not be sufficient to cause saturation and distortion 
of the Variac controls. 

It was determined experimentally that a suitable 
no-load terminal voltage curve produces not over 120 
volts at 60 cycles per sec., and the alternator is oper- 
ated at a fixed excitation. 
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SYSTEM FOR EBACITATION 

The approximate relationships of the voltages and 
currents in the system may be better visualized by 
reference to Fig. 4, representing a single shaker con- 
nected to an alternator by means of a single Variac. 
(The second Variac, shown in the diagram of Fig. 10, 


may be considered to be set for a transformation ratio 


of unity. ) 

In the section on ‘Shaker Design Considerations,’ it 
was pointed out that one of the main factors govern- 
ing the system design was the characteristic of the 
structure being shaken. It was further pointed out, 
in Eq. (26), that the electrical ratio of K,/I was di- 
rectly related to the mechanical ratio e/7 through the 


force-current characteristic of the shaker. That is, 


E,/I = 0.1125 (i/D)? e/i 


If a force i is applied to a certain structure at its 
resonant frequency, and the structure responds with a 
velocity e, then, since the characteristics of the electro- 
mechanical transducer or shaker are known, the elec- 
trical quantities / and / can be calculated as follows: 


I = s/@/f) (27) 


E, = 0.1125(4/D)e (28) 

The voltage applied to the shaker terminals, /s;, must 
be equal to the back voltage /, plus the voltage drop 
through the armature coil, Re. 


Es = Ey + IRe 


(29) 


The shaker voltage is supplied through the auto- 
transformer or Variac the alternator terminal 
voltage /;, which is related to Es by the Variac setting 
N. as the alternator current is re- 


from 


or turns ratio N, 
lated to the shaker current. 


Ey = E(Ni/ Ne) = (Ni/Ne) (Eg + 1Re) (80) 
I, = I(N2/N;) (31) 
Dividing Eq. (30) by Eq. (31) gives 
nn tx) ~ Qu) Gite) 
or 
Rr = Rs(Ni/Ne)? = (Mi/Ne2)?(Reg + Re) (33) 


That is, the terminating resistance of the alternator is 
equal to the effective shaker resistance times the square 
of the turns ratio set on the Variac. Also, it is seen 
that the effective shaker resistance, Rs, is equal to the 
reflected resistance of the mechanical load, R,, plus 
the coil resistance, Ro. 

R, is therefore a function of the variables N,/N2 
and R,; Rr, can be held constant for different values 
of Rz by means of a compensating setting of the Variac 
Ni/No. In addition, the terminal voltage E, differs 
from the no-load voltage Ey, because of the internal 
resistance of the alternator, R,, complicated by satura- 
tion and armature reaction effects. 


OF PURE NATURAL MODES 
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It can now be seen that, for a given frequency of the 
alternator producing a certain no-load voltage Ey,, 
there will be a certain value of R;, which cannot be re- 
duced without causing the current rating of the alter- 
nator to be exceeded. 

i. = Exe (34) 
(Racy. + Rrmin.) 

It can also be seen that, if Rz is low (a highly damped 
structure), the setting N2 can be decreased [Eq. (33)], 
causing R, to remain above the minimum value. 

These relationships can be visualized as a function 
of varying alternator frequency by reference to Fig. 5. 
It will be noted that, for a certain assumed current, 
there is a certain frequency, fi, at which all the avail- 
able voltage, yz, is dissipated by losses, leaving none 
available to overcome Ey. 

This simply means that a current equal to (No/Ni)Jq 
can be pushed through the shaker coil at the frequency 


fi only when it does not move (blocked armature) and 


that, if the armature is allowed to act on the resonant 
structure, the current deliverable must be less. 

It is obvious that the most critical design condition 
is at low frequency, where there is little voltage avail- 
able from the alternator to overcome the voltage drop 
at high current values. 

An understanding of the relationships between cur- 
rent, back voltage, power to the mechanical system, 
and losses is better realized by examination of Fig. 6, 
which shows these relationships at 4 cycles per sec. 


alternator frequency. The line /:, on this figure 
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i ; a particular value of /, at which Py is a maximum. 
voltage drop due to [Rg losses in the armature coil and |, sane ae 5 
; : From the P,, curve of Fig. 6, it is seen that a maximum J 
its connecting cable, as reflected to the alternator by 
the Variac turns ratio Ni/N. = 1. For any given 
current, then, the ordinate between these lines is equal 


power transfer to the mechanical system is obtained 
when /, = 2.1 amp., at which time /,(.Ni/N») 1s 
; . ‘. 3 " .32 volts. 
to E,(.N,/ Ne), as can be seen from Eq. (30). 
Now, the power delivered to the mechanical struc The lower plot of Fig. 6 shows the approximate per : 
I I 4 PI I f 
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SYSTEM FOR EXCITATION 





Experimental structure with variable damping, set up 
for test. 


Fic, 7. 


sec. It will be noted that, with the turns ratio setting 
of N,/N2 = 1 as shown, the maximum power cannot 
be delivered because the shaker stroke limit of 0.5 in. 
gives a maximum back voltage of only about 3.4 volts 
rather than the 4.32 volts required at maximum power. 


Similar system limits of power delivery and EF, vs. / 
relationships can be drawn for other turns ratios within 


the possible range of about 0.8 to 50. Because of 


OF PURE NATURAL MODES 


EXPERIMENTAL DEVELOPMENT 


A test structure was devised to allow experimental 
evolution of a system whereby natural modes might 
be excited by application of orderly methods based on 
the concepts described earlier. In particular, it was 
desired to know: 

(1) How good an approximation to a natural mode 
could be produced in a uniform structure when forces 
were adjusted proportional to the product of mass and 
amplitude at the station of force application. 

(2) How much better this excitation was than some 
more arbitrary excitation. 

(3) How satisfactory an excitation of pure modes 
might be obtained by the described method of force ad- 
justment, even in cases of nonuniform structure with 
irregular damping characteristics. 

The test structure as set up for testing is diagrammed 
schematically at the head of Fig. 8. The structure 
consisted of nine lumped masses interconnected by two 
continuous straps or flexures bolted to the top of the 
masses at each side and two identical flexures bolted 
to the bottom of the masses on each side. By this ar- 
rangement, the flexures constrain the masses to trans- 
late relative to each other along axes perpendicular to 
the flexures. 

Motion of one mass relative to the next mass was re- 
sisted not only by the flexure forces but also by damping 
forces that were produced electrodynamically. Each 
mass, except the center one, consisted of a field struc- 











the alternator characteristic, however, only slightly 
higher limits than those indicated for N,/N2 = 1 are 
obtainable. 
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FIFTH MODE ASYMMETRIC STRUCTURE 
f+ 45.5 cps. 
9= 0.044 
“A= 0.930 
DECAY 


“METHOD” EXCITATION 


STEADY STATE 
“METHOD” EXCITATION 


STEADY STATE 
STA. | EXCITATION 


STEADY STATE 
STA. 9 EXCITATION 


Fic. 9. 
the following conditions: 
at Station 1; bottom 


top 


ture and field coil identical with those used for the ex- 
perimental shakers. In the air gap of each mass other 
than the center one was a single shorted turn or tube 


as armature, which was rigidly attached to the next in- 


board mass. Thus, motion of one mass relative to 
another caused damping forces between the two 
masses. 

The mass distribution of the test structure was 


changed by (a) adding lead blocks to the original 
mass and (b) removing the field coils. In the latter 
case, of course, the damping between the light mass 
and the next inboard mass was eliminated. 

‘ The damping distribution was changed by (a) re- 
moving the damper field excitation, (b) introducing a 
mixture of iron powder and oil into the flux gap, and 
(c) removing selected field coils, as in (b) above. 

It was also possible to change the damping by a 
change in thickness of the damper tube or shorted turn. 
This variation was not found to be necessary. 

Tests were made on two variations of the test struc- 
ture: (a) a uniform structure, as diagrammed at the 
head of Fig. 8; (b) an asymmetrical structure, as dia- 
grammed at the head of Fig. 9. 

In the former case, masses, stiffnesses, and dampers 
were as near uniform as practicable, and, in the latter 
case, masses and dampers were made to differ in mag- 
nitude and placement in an asymmetrical fashion, with 
the stiffnesses remaining uniform between masses. 

The test structure was supported from an overhead 
frame at a low frequency by means of springs attached 





Data obtained on the fifth natural mode of the asymmetric structure for various excitations. 
decay of structural mode from ‘‘method”’ excitation; 
decay of same mode from excitation at Stations 1 and 9. 
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Oscillograph records show ‘ 


middle—decay of same mode from excitatior 


near the nodes in the first free-free mode of the struc. 
ture. A vibration pickup was attached to each mass 
element of the structure for measuring the steady-state 
amplitude at each mass station. A shaker was attached 
to each mass and was mounted below the test struc. 
ture. Each shaker was controlled from a_ console, 

whose arrangement and auxiliary systems were evolved 

during the experimental program to allow systematic 

convergence excitation. This | 
system, with refinements, is to be described later in con- 

nection with the final equipment, and further reference 

to it here will be omitted. 


on a natural mode 


The test procedure, described more fully in reference | 
6, was to excite the test structure, first, by the method 
of pure mode synthesis devised and, then, in one or { 
more arbitrary ways, but using good technique. In 
each case the amplitude of each mass was measured 
by means of its pickup, and phase relative to some 
reference was established. A record was also taken ol 
the mode decay when the excitation was removed. 


Experimental Results 


Results of these tests are exemplified by Fig. 8, which 
shows comparative data for the various excitations 
applied to the sixth mode of the uniform structure, 
and Fig. 9, which shows similar data for the fifth mode 
of the asymmetric structure. 


Evidence of the superiority of the multiple shaker | 
adjustment is obvious in the purity of the decay record | 
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SYSTEM FOR EXCITATION 
and by the close degree of similarity between the steady- 
state and decay amplitude ratios obtained by the mul- 
tiple excitation. 

On the other hand, the amplitude ratios obtained from 
the arbitrary excitations bear little similarity to those 
in the pure mode, and the large phase shifts and poor 
decays give some indication of the obscurity in which 
data of this type would leave the analyst. 

Referring to Fig. 8, some similarity between the 
mode shapes obtained in the various ways is observed 
on comparison of the absolute amplitudes of the 
masses. These values alone are of little value, however, 
especially on a continuous structure, unless the posi- 
tion of the nodes, indicating a reversal of phase, can be 
ascertained. The large phase shifts indicate the ab- 
sence of true nodes on the structure and, hence, makes 
their location doubtful and only approximated pain- 
fully by careful amplitude surveys or phase analyses. 
The mode is obviously there but makes concise definition 
difficult, as shown in the poor decays. 

Specifically, it can be concluded that the system of 
multiple shakers, adjusted according to a scheme based 
on the assumption of damping distributed in propor- 
tion to stiffness, offers a great clarification in direct 
experimental evaluation of mode properties not only 
for a structure that has the assumed damping distribu- 
tion but also in at least some cases where the actual 
radically from that 


damping distribution departs 


assumption. 


Approach to the Experimental Problem 


The theoretical analysis indicated that, with an as- 
sumed simple distribution of damping and with multiple 
shakers distributed to act on each element of the equiva- 
lent lumped system, the natural modes could be ex- 
cited by the adjustment of each shaker force to be pro- 
portional to the product of the assigned mass and its 


amplitude in the mode under excitation. That is, 


cmp; 


Since the mass assigned to any shaker is a constant for 


1 = 


a particular setup, it can be introduced as an electri- 
cal constant, attenuating the signal that must represent 
The 


shaker force is proportional to its armature current, 


the magnitude of the adjusted shaker force. 


which is regulated by the voltage across its control 
Variac. This voltage is attenuated by a potentiometer 
calibrated directly in mass ratio, and the resultant 
signal for each station, proportional to i/m, is placed 
on a switch whose selector goes to the horizontal axis 
of a cathode-ray oscillograph. 

The amplitude of the mass in the mode f; is measured 
by an accelerometer pickup attached directly to the 
shaker armature, which moves the mass. Its output 
is integrated to give proportionality to velocity, and the 
resultant signal for each station is applied to a switch 
Whose selector goes to a normalizer and then to the 
vertical axis of the scope. 


OF 


NATURAL MODES 


“I 
wo 


PURE 


Proof that ratioing Variac voltage to pickup output 
is equivalent to ratioing force to amplitude is as follows. 
Since mass is a constant for each station, it will not be 
considered in the proof. 


Let 

Es, = Variac voltage at any station 

yn, = pickup output voltage at any station 

I, = current in voice coil of shaker at station 
(amp. ) 

i/I = calibration of shaker = constant (Ibs. per 

amp.) 
1, = shakerforce = (z/J) I, (Ibs.) 
S, = amplitude of motion at station n (in.) 


Then, following Eq. (29), 


Esp = Ez, fe I,,R 


and 
Ey, = 0.1125 (i/De, = 0.1125 (i/D2rfs, (39)* 
for a given mode 0.1125(¢//)2xf = constant A. Then, 
Es, = As, + [R/(t/I)|t, (40) 
E,», = calibration constant: acceleration 
= (const.).(27f)*s, = (const.)-s, 
E,, = Bs, (41) 
then, 
Es, As, | [R/(i/I)]1, 
E,, Bs, Bs, 
Es,/ Ep, = (A/B) + [(R/B)/(4/D] (tn/sp) (42) 
If the adjustment is made such that Es,/F,, = ¢ 
forn = 1,2,3,..., thenz,,/s, will have a constant ratio 
1n/Sn = (4/I)/(R/B)[c — (A/B)] (43) 


Procedure in Exciting a Natural Mode 


Operation of the system in exciting a natural mode 
starts out with the introduction of alternating forces to 
the structure by means of only one of the many shakers 
driven from the same system. The frequency of ex- 
citation is varied until the structural response indi- 
The 
electrical signals corresponding to the amplitude and 


cates proximity to a natural mode frequency. 


force at this station are selected on the cathode-ray 
oscilloscope. As the resonant frequency is approached, 
the velocity of the structure approaches an in-phase 
condition with the force; hence, the observed figure 
is an ellipse of variable fineness ratio. The gain con- 
trols on the scope are adjusted to give the major axis 
of this ellipse a 45° slope. By selecting another sta- 
tion, its force may be adjusted proportional to its am- 
plitude merely by adjusting the Variac to obtain a 
45° slope of the observed ellipse. 

Energy introduced into the system by the second force 
will cause the amplitude of response to increase, thereby 


* There are no Fqs. (36)—(38). 
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changing the ratio of force to amplitude at the first 
station. The purpose of the normalizer is to maintain 
the signal to the scope constant (producing a 45° ellipse) 
for the changed ratio of amplitude to force at the first 
(or master) station. This is done by dividing the in- 
creased pickup output of the first station by a quantity 
sufficient to re-establish the original ratio of signals to 
the oscilloscope and dividing any other selected pickup 
output by the same quantity. 
normalization, automatically performed as described 
later. 
into adjustment to succeedingly better approximations 


This is essentially a 
By this means, each successive force is brought 


to the finally required constant ratio of force to ampli- 
tude, and readjustment of only the first few forces is 
required to obtain the desired condition : 


(1/m), (1/m)s (1/m)s 


C 


Sy So Se 


Before turning on another shaker at a selected sta- 
tion, it is necessary to make sure that its force, when 
applied, will act in the proper sense, putting energy 
into the system rather than drawing it out. Applica- 
tion of the original force at the master station evoked a 
velocity response of the structure that was in phase at 
that station, and the phase of every other point in the 
structure (in the natural mode in question) is either in 
the 
applied is measured by the Variac voltage and the phase 


or out of phase with this motion. Since force 
of all Variacs is the same, the phase of the structural 
response at any station can be found by comparing the 
velocity signal from the pickup with the Variac voltage 
before the shaker circuit is closed. The wiring can then 
be arranged so that an in-phase indication will produce 
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a lissajous figure to the right, calling for a phase-switg, 
position to the right. Similarly, a lissajous figure With 
axis inclined to the left requires a phase-switch setting 
to the left. 


can be turned on and adjusted to level as previous) | 


After this adjustment is made, the Shaker 





A 
described. b. 
It is an experimental fact that, although phase Shifts 
between the responses at various stations may be large ' 
under the excitation of the first applied force even whey 
adjusted to a resonant frequency, as more and more | 
forces are applied and adjusted the more uniform jy | 
phase are the various points on the structure. This js 
to be expected. Since we end up with adjustments 
urging the structure at each of 24 points to move in the 
natural mode form, we have an excitation: practical) 
orthogonal to all other modes. 
This means that, with the shakers so adjusted as to 
phase and force amplitude, the structural response wil FIG 
with - 
Note 


be in the one mode regardless of the frequency of excita. 


tion. 
obtat 
the s 
and 


serve 





Fic. 11. Control console for systematic excitation of pure 
natural modes. Contains all power supplies, instruments, and 
control circuits for 24 shakers and associated pickups. Data 
recorded by oscillograph housed in bay at right. 


Shaker with accelerometer attached. 
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SYSTEM 





Fic. 13. Chassis containing all master and channel selectors, 
with 24 pickup coupling circuits and sensitivity adjustments 


Note: Dash numbers refer to drawing No. 17CO447 


tion. The frequency can, of course, be adjusted to 
obtain a maximum velocity amplitude as indicated on 
the scope or to obtain an in-phase adjustment of force 
and velocity as indicated by degeneration of the ob- 


served ellipse to a straight line. 


THE FULL-SCALE EQUIPMENT 
A schematic diagram of the overall system is shown 
in Fig. 10. The (exclusive of 
shakers) was housed in a control console requiring 
but a single connection to a 220 volt, three-phase power 


complete system 


line. 

Thirty shakers were provided, any 24 of which could 
be controlled from the console. All cables have stow- 
age space in the console. The essential features of the 
variable frequency power supply and distribution sys- 
‘tem have been discussed earlier, and it is now in order 
to describe some of the more important features of the 
instrumentation and control circuits provided. 


Pickups 

A pickup is attached directly to each shaker arma- 
ture. Thus, with a single attachment to the structure, 
the force applied and the amplitude resulting at the 
point are obtained. 

The pickup signal is used in two ways: (1) to provide 
an amplitude signal for reference in the adjustment of 
shaker force and (2) to provide simultaneous signals 
for recording the resultant structural motion on a multi- 
channel oscillograph. 

The pickup used was of a new type, allowing the 
operation of the oscillograph elements without ampli- 
fication, which affected a considerable degree of sim- 
plification in the instrumentation and its maintenance. 
Coupling circuits for each of the 24 channels are car- 
ned on a single rack mounting panel. Provision is 
made for convenient calibration of the pickups either 
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on an absolute or comparative basis. The latter is 


mode-shape measurements. 


generally sufficient for 
Sensitivity adjustments are provided for each channel, 
and means are also provided for rapid overall electrical 
calibrations of the several attenuator-oscillograph ele- 
ment combinations. 

The pickups used are of 
275 cycles per sec. natural frequency. 
sensitivity exceeds 0.5 volt per g and, when coupled to 
the oscillograph, give signals flat in amplitude response 
to the required 100 cycles per sec. with a sensitivity of 
25 in. of record deflection per g, thus making adequately 
clear the need for attenuation at reasonable accelera- 
tion levels and the ability to produce good record de- 
flections at low structural accelerations. 


l5g capacity and about 
Their voltage 


Normalizer 


A simplified diagram of the normalizer is shown in 
Fig. 15. The amplitude and force signals from any 
particular station, selected as master, are fed to the 
normalizer. 

The pickup output from this station is applied to the 
master potentiometer, from whence the attenuated 
signal goes through a current amplifier and is applied 
to the heater element of an isolated construction 
vacuum thermocouple. At the same time, the output 
of the miass-ratioing divided 
by mass) is applied to a similar thermocouple. The 


potentiometer (force 





Left panel of control console, showing typical con 
ind 


Fic. 14 
trol stations, frequency control, scope indicator, master 
channel selectors, and normalizet 
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Fic. 15. Simplified schematic diagram of normalizer. 
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-3b _ hidden 


sensitivity 


adjustment 


1700430 Amplifier 


2 -21 


Fic. 16. Normalizer chassis, showing servo amplifier and 
motor, dual potentiometer, thermocouples, and autématic over- 
load protectors. Note: Dash numbers refer to drawing No. 
17D0431. 


outputs of the two thermocouples are connected in 
series opposition and are applied to a d.c. servo-ampli- 
fier and associated balancing motor. If the current 
through each thermocouple heater is the same, it is seen 
that their outputs cancel and there is no input to the 
servo-amplifier. If the current through the pickup 
thermocouple exceeds that through. the force signal 
thermocouple, a polarized signal is applied to the servo- 
amplifier and causes the balancing motor to rotate in a 
certain direction. The master potentiometer is geared 
to the motor in such a manner that its output is further 
attenuated until the outputs of the two thermocouples 
again cancel. Conversely, if the output of the force 
signal thermocouple exceeds that of the pickup thermo- 
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couple, a signal of the opposite polarity is applied t 
the amplifier motor combination and the direc 


tion of 
rotation is reversed. 


The master potentiometer js now 


| 


driven so that it attenuates the pickup signal to a lesser § 


extent until the system again comes into balance. The ' 


normalizing potentiometer is identical to the Master 
potentiometer and is geared to the motor so that its 
motion is the same as that of the master potention. 
eter. Thus, the output of any pickup applied to th 
normalizer potentiometer through the pickup selecto; 
switch is divided in the same ratio and is thereby ef. 
fectively normalized. 


The normalizer assembly utilizes a standard an¢ 


commercially available servomotor and amplifier, 




















11, 12 (Pitching engine) 


Shaker Channel No Assigned Weight Mass Patio 
Pa es 75 Ibs. 1.14 
Bs Gi Ble oe 152 2.8 
7,38 1996 19.5 
o, ©, 6, % 6 1.0 
&, 17 450 6.8 
13, 14 1707 S35 
uw, 7 _— (15) 
Note: Value of Iass ratio in ( ) for shakers 11, 12, 13 & 14 actually used. 


Fic. 17. Shaker locations and mass ratios for SB2C-5 sym- 


metric modes, 





Fic. 18. 


Shakers on mounting stands beneath wing. 
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Provision for Setup 

The first actual tryout of the equipment was made, 
after delivery, on an SB2C-5 airplane. A diagram of 
the shaker locations and mass ratio assignments is 
shown in Fig. 17. To facilitate locating the large 
number of shakers used on such a setup, a demountable 
stand was provided. The top section consisted of a 
plywood pyramid with a metal mounting table adjust- 
able in height from approximately 25 to 43 in. on a 
coarse lead screw provided with a clamp. Either one 
or two subsections, each assembled from flat plywood 
sections with trunk latches, could be used, thus pro- 


viding a maximum height of shaker attachment of 83 in. 


Attachment to the structure was made by means of 
metal plates provided with tapped bushings and either 
cemented to the structure over a thin intermediary 
layer of sponge rubber or screwed to existing convenient 


‘ 


locations. 


Connection between the shaker-mounted accelerom- 
eter and the structure was made by means of the rod 
and flexure arrangement. The rod and flexure links 
provide a simple means of compensating for misalign- 
ment and provide a safety link to prevent inadvertent 
overload of the shaker armature by excessive motion of 
the aircraft, which is generally spring-supported. 


Initial Full-Scale Tests 


Records of three modes set up during the first test 
runs are shown in Fig. 21. These are acceleration rec- 
ords, and it will be noted that the amplitude ratios are 
clearly definable in either the steady-state or the decay 
portions of the records. The modes meet the criterion 
of purity, as evidenced by the lack of beats in the de- 
cays. 

Some hash is in evidence on the acceleration traces, 
particularly due to excitation of the transient caused by 
abrupt removal of the shaker excitations. Some of this 
hash is due to pickup of miscellaneous structural buzz 
and rattle, and some is due to excitation of the pickup 
natural frequency. It should be pointed out in this 
connection that use of a phase criterion in setting up the 
modes requires uniformity of phase response, as well 
as a minimized phase shift. For this reason, the pickups 
are not damped to the conventional 65 per cent of criti- 
cal. In any case, the hash does not seem to be of serious 
consequence in view of the clarity of the recorded 
modes. 


Comments on Operation of Equipment 


It is worth remarking that the modes recorded in 
Fig. 21 were set up, or synthesized, experimentally 
without any direct reference to the aircraft structure. 
Rather, the shaker force adjustments were made at the 
control console solely by reference to the cathode-ray 
oscilloscope for information as to required adjustments 
in frequency, phase of shakers, or force. 
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Fic. 19. Attachment of shakers to structure 


en, i 








demountable shaker stand 


Fic. 20. Adjustable and 

It is further worth remembering that this orderly 
experimental method of natural mode synthesis bases 
its action in force adjustment on a specific and simple 
assumption of damping distribution. Its virtue, there- 
fore, must not be expected to lie in either its exactitude 
or its infallibility but rather in its ability to yield any 
useful result in a complicated experimental situation. 
This can only be fully appreciated by one who is faced 
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MG. 2 
PPh Sona for modes at 8.4, 14.7, and 26.5 cycles per sec. 
100. 
with the myriad uncoordinated adjustments possible in 
a multiple shaker setup. 

The criteria of any natural mode, whether excited 
by a single- or a multiple-shaker system, are in its 
ability to exist in the steady state without phase shifts 
or to decay without beats. In many practical cases, 
moderate phase shift or beating is not intolerable. 
Should the orderly system devised fail in this practical 
test, it will do so because the adjustments made are 
called for by an assumption that is violated. It then 
becomes a matter of individual insight and common 

sense on the part of the operator in modifying these 
adjustments better to meet the criteria. 
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Recorded structural modes of the SB2C-5 airplane excited with equipment. 
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3 


Freedom from beats is evident in the de 


Attenuation on the various traces of acceleration response varies from | to 


It is of particular importance to note that the operat- 
ing adjustments are to be based on phase comparisons. 
If and when a single mode is excited to the major ex- 
clusion of all others, the motion of all points will be in 
(or 180° out of) phase, regardless of the proximity to 
the natural mode frequency. It is therefore evident 
that a valid experimental evaluation of the in-phase 
condition requires a high degree of uniformity in the 
phase response of the pickup circuits, which must be 
verified. The in-phase condition of force and velocity 
is conveniently used as a criterion for the frequency 
setting. 


It must be realized that the absolute phase | 
shift of the circuits may be appreciable even though the 
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SYSTEM FOR EAZACITATION 
ive phase shift between circuits is negligible, and, 


relat , 
fore, such frequency adjustments are not accu- 


there 
rately the mode natural frequency. This in itself is not 
too important, for, if the excitation is orthogonal, the 
mode response will still be in phase, but not at maxi- 
mum amplitude, and the true frequency will become 
evident on the decay. Complication with regard to 
the frequency adjustment arises in actual cases where 
two modes are close together in frequency. 

In following the procedure set up in the experimental 
method, it is desired to make convergence on one or the 
other of these modes. Bringing in successive shakers 
tends to produce large phase shifts in the responses 
because of elimination of one of the near-resonant 
modes. If, under these conditions, readjustment of the 
frequency is made to close the ellipse formed by force 
and velocity signals, the change in frequency may be 
sufficient to cause a node line to shift from one side of a 
shaker to another because of the effect of changed 
frequency on the two as yet unseparated modes. 
Under these conditions the incorrectly phased shaker 
will draw power from the system. Either the adjust- 
ments become divergent or the frequency 
changed and the other mode chased inconclusively. 


is again 


This tendency must be guarded against by awareness 
of the operator: (1) that both frequency and force dis- 
tribution are changing the indicated phases; (2) that 
the indicated phases are not necessarily absolute and, 
therefore, that the true mode frequency is not necessar- 
ily obtained with the force-velocity ellipse closed; and 
(3) that, even without the true mode frequency, con- 
vergence on one mode will be indicated by the same 
phase shift at all stations. 


CONCLUSIONS 


A coordinated set of equipment has been provided 
for the experimental evaluation of the pure modes 
of complex structure. The equipment has been made 
as simple in operation as possible under the complex 
requirements of its use. 

It has been shown that the multiple shaker system 
may be adjusted to excite pure natural modes in ac- 
cordance with a simple procedure, utilizing an automatic 
normalizer and signals readily interpreted as to phase, 
frequency, and amplitude on a cathode-ray oscillo- 
scope. 

If and when this method, based on a simplified as- 
sumption of damping distribution, fails to produce a 
pure mode, it will be evident on the directly recorded 
data. In such cases, the adjustments may be modi- 
fied with full ability to evaluate the effects on the 
structure. 

Much careful attention has been given to providing a 
complete set of equipment in which the power circuits 
and instrument circuits are basically simple and utilize 
commercially available components. Coordination of 
the extensive equipment required into a single console 
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has helped to provide quick and reliable means of overall 
check-out on the circuits prior to or during tests. 

It is hoped that the emphasis on application of ex- 
perimental techniques during the test period wili be 
amply repaid by the quality of the data obtained and the 
savings in data analysis and reduction time. 


APPENDIX——NORMAL MODE THEORY 


A complete mathematical treatment of the theory 
of small oscillations of conservative systems and of 
systems with small damping may be found in reference 
3, chapters 5 and 6, and reference 4, Section 4. The 
results of these analyses have been summarized below. 
Reference 5 was an invaluable aid in making this sum- 


mary. 


Conservative Systems 


Free Oscillation.—Lagrange's equation of motion 


for the free vibrations of a conservative system is 
(d/dt) (OT'/0g;) + (0U/eq;) = Q; = 0 


where 7 and U are kinetic and potential energies, re- 
spectively, g; are the generalized coordinates (7 = 1, 
2,...mn), and Q; are the corresponding generalized ex- 
ternal forces (Q; = 0 for free oscillation). 

The most general form of oscillation of such a system, 
with static coupling only, is represented by the super- 
position of 7 pure harmonics, 


n 


> A; sin (w,t + y,) 


di = 


If the amplitudes of the principal oscillations are 
introduced by putting A; = ¢,°"c,, it is possible to 
evaluate 7 and U in terms of the principal oscillations. 
Further, using the normalizing condition 


n 


> mo,”?? = M 


i 1 
and calling 
& = c, sin (w,t + y,) 


the displacement of the rth normal mode, then, 


T = (1/2)M > £,? 
r= i 


n 
U = (1/2)M >> ow, 
r 1 
and the equations of motion become simply 


£, + w,*§, = 0 


The kinetic and potential energies of the original 
system are seen to be the sums of the kinetic and poten- 
tial energies of m single degree of freedom systems of 


9 


mass ./ with a spring constant A, = Ma, 
Forced Oscillations.—The equations of motion for 
the same system under the action of periodic external 


forces are 
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d (=) ‘ Ou 
dt \Og; Od: 


where F; is the amplitude of the 7th component of the 
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= F, sin (wt + y) 


generalized force. The component of F; corresponding 


to the displacement of the rth normal mode is 


n 


L Fe” 


i= 1 


o. = 


The equations of motion in terms of the normal co- 
ordinates £, and the normal force components ®, are 
» 

é, ee w,*&, = 


The steady-state solution is 


(,/M) sin (wt + y) 


®, ] , 
sin (wt + yp) 


"* K, 1 — (w/w,)? 


If w? — w,*, resonance occurs, and the amplitudes 
of all normal modes are small compared to the rth mode. 
The amplitude ratios of the forced vibration are ap- 
proximately equal to the amplitude ratios of the rth 
normal mode. 


Damping (Nonconservative Systems) 


If, now, forces that are linear 
functions of the velocity components q; (i = 1, 2,..., ) 


Viscous Damping. 


are introduced to the system described above, La- 
grange’s equation becomes 


d (") 4 oD rs Ou 0 
dt \Q¢, a 


where 7 and U are kinetic and potential energies as 
before. D is the dissipation function, and Q, are the 
generalized external forces corresponding to the gener- 
alized coordinate q;. 

In terms of the normal coordinates of the undamped 
system, the equations of motion for free vibrations 
are 


Mi,+ DY Brot: + Kt, = 0 
s=] 


were 8,, are constant damping coefficients. Further- 
more, if the dissipation is small (i.e., B,, < M and £,, 
< K,), the rth solution of Eq. (10) is largely made 
up of £,, which is the solution of 


Mt, + Brt, + Kt, = 0 
and small £,’s whose amplitude and phase are given by 
f, = &[8,5 w,/(Mw,* + K;)] 
and w, is the complex root of 
M w,? + Bw, + K, = 0 


The imaginary part of this complex root gives the 
decay frequency that is nearly the undamped natural 
frequency. The real part is a measure of the rate of 
decay. 


AERONAUTICAL 


SCIENCES—NOVEMBER, 1950 


| 

In the case of forced oscillations of the slightly i 

damped structure (Q; = F; sin wt), the equations of t 

motion in terms of the normal coordinates of the gop. 
servative system are 


n 
Mé, + Z, Brsks + K,é, = ®, sin wt 
s=1 
Again, for small values of 8,,, the cross terms 8,, are 
virtually negligible, and the solution is predominantly 


®, l 


& = 
K 
[ 


&; : sin wt + @ 
6 = tan 


“OT 
W, K, 
' {(Brr @/K,)/[L — (w/c)? | 


This is recognized as the solution for the forced 0. 
cillation of a single degree of freedom system with vis. 
cous damping. 

Special Case. 
damping that is proportional to, and in quadrature 


Of particular interest is the type of 


with, displacement. The equation of motion in terms 


of the normal coordinates of the conservative system 
are 


n 
Mt, + Ki(& +7 SO gpk) = &, sin of 
AY = 1 
Again, for small dissipation (i.e., g,, < 1), the rth 
solution is predominantly the solution of 
Me, + K,(1 + jg,)é = ®, sin wt 


For the free oscillation (®, = 0), the decay frequency 
V K,/M and the rate of 
decay is approximately equal to that of a single degree 


will be approximately w, = 


of freedom system with a viscous damping coefficient 
Ber Bevis a 4 2) 2Zrr- ; 
The solution of the equation for forced vibration is 


®, l ; 
SS sin (wt + 6) 


A. w\2 2 
(OT. 
W, 


@ = tan 


Str 


r 


{ grr/[(w/w,)? — 1]} 
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The Flow over a Wedge Profile at 
Mach Number | 


G. GUDERLEY* ann H. YOSHIHARAT 
Wright-Patterson Air Force Base 


SUMMARY 


The flow over a wedge-shaped profile at zero angle of attack 
ind with a free-stream Mach Number of 1 is computed using 
the hodograph method as simplified by the transonic law of simi- 


larity. Drag coefficients are computed for several profiles 


(1) INTRODUCTION 


bbe FLOW PATTERN Over a two-dimensional profile 
with the free-stream Mach Number of | has been 
discussed previously.'_ The example given there showed 
that steady flow patterns at Mach Number | do exist. 
However, in the example computed there, the con- 
tour of the body is not suitable to carry out further 
e.g., to determine the influence of an 

For this purpose the study of a wedge 
Furthermore, such a prob- 


investigations 
angle of attack. 
profile is more promising. 
lem represents an example where the boundaries in the 
physical plane are predetermined so that it does have 
an interest of its own. 

The formulation of the boundary value problem in 
the hodograph plane does not offer essential difficulties. 
The singularity that arises at the point of the hodograph 
which corresponds to the free-stream velocity has 
already been investigated.' Also, a family of particular 
solutions suitable to fulfill the boundary conditions by a 
superposition is easily found. Serious difficulties, how- 
ever, are found when one tries to carry out this super- 
position such that the boundary conditions are actually 
fulfilled. A direct attack leads to an infinite system 
of equations. One might attempt to fulfill this system 
directly by taking a finite number of terms in the sys- 
tem; however, this procedure did not prove too satis- 
factory. 

The present paper shows how to overcome this diffi- 
culty. First, the boundary value problem is changed 
in such a way that the solution remains the same while 
the supersonic part of the boundary is more conveniently 
located. This change does not yet eliminate the 
necessity of solving an infinite system of equations, 
but it reduces considerably the amount of work re- 
quired to establish this system. Furthermore, in the 
present case of zero angle of attack, the problem can be 
formulated in terms of an integral equation whose kernel 
in a close approximation corresponds to that appearing 
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in the theory of a thin airfoil in an incompressible 
flow. Thus, the solution finally requires only a suitable 
Fourier analysis of the boundary conditions. 
(2) Basic RELATIONS 

Let x and y be the Cartesian coordinates in the physi- 
cal plane; w, the absolute magnitude of the velocity 
vector; 96, the inclination of the velocity vector with 
the x-axis; p, the pressure; p, the density; a, the ve- 
locity of sound; y, the stream function; y, Legendre’s 
transformed potential; w*, the critical velocity at which 
a; and k the ratio of the specific heats (k = 


con- 


w= 
stant = 1.4). 

If one introduces the quantity » by means of the 
relation 


y*) w* } 


[(w — w (1) 


n = (k + 1) 
then the differential equation for the stream function 
simplified by the transonic law of similarity assumes the 


form 


Von = nWoo = 0 
Similarly, the transformed potential can be expressed as 


Yan — 109 = O (3) 


The values of the x and y coordinates in the physical 
plane corresponding to a certain point (ny, @) in the 
hodograph plane, within the above approximation, are 
given by the expressions 
[(k + 1) 


(1 agw* )Qa f 


x(n, 0) = ‘/ayw* |e, t (4) 


y(n, 0) = 
Here, dy» is the velocity of sound for the stagnation con- 


dition. 


In general, the boundary value problem in the hodo- 
graph plane will be given in terms of the stream func- 
tion. This formulation has the advantage that along 
the boundaries one has conditions that correspond to a 
standard form of boundary value problems—that is, 
to a case in which the value of the function to be deter- 
mined is prescribed at the boundaries. However, be- 


cause of the relation 
p* oo (5) 


y = 


which connects the transformed potential with the 
stream function, the boundary value problem can be 
easily expressed in terms of the transformed potential. 
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STREAMLINES 


Fic. la (top). Flow in the physical plane. Fic. 1b (bottom). 


Flow in the hodograph plane. 


This form is more suitable in the transformation back 
to the physical plane. 

The pressure coefficient C,, is expressed in terms of 7 
by means of the equation 

C, = (p — p*)/(p*-w*?/2) = —2(k + 1)~ (6) 
This form is consistent with the above approximations. 

The characteristics in the hodograph plane of the 
Eq. (2) or (3) are given by 
(7) 


@ x (2 3)n° ‘= const. 


where the plus and the minus signs refer to the two pos- 


sible families. The corresponding directions of the 
Mach waves in the physical plane are given by 
dy/dx = =(k+ 1) '%_-" (S) 


Here, the signs are to be used consistently with those 
of Eq. (7). 
of » represent supersonic velocities, while negative 


Note also from Eq. (1) that positive values 


values correspond to subsonic velocities, so that only 
the former give real characteristics as is obviously re- 
quired. 


(3) THe BOUNDARY VALUE PROBLEM FOR THE WEDGE 


The hodograph representation of the boundary value 
problem for the flow over a wedge with the free-stream 
Mach Number of 1 is given in Fig. 1b; the flow in the 
physical plane is given in Fig. la (In the following 
description only the flow in the upper half-plane will be 
described, since the flow in the lower half is the same.) 

All stream lines start at negative infinity of the 
physical plane with the free-stream sonic velocity. 
Accordingly, they all begin at point A, which represents 
a singular point. 
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The upper half of the physical plane is limited by the 
zero stream line. Let us follow the zero stream line it 
the hodograph plane. It starts at infinity with th 
sonic velocity (point A). As one proceeds along th 
zero stream line toward the profile, the velocity veetg, 
will always have a horizontal direction, and its magni 
tude will decrease. Accordingly, the map of this por 


tion will lie along the line AB. At the nose of th 





wedge one has a stagnation point; the nose of the pro 
file will therefore correspond to point B. Here ,§ 
branch point of the stream lines arises, and the zer 
stream line then follows the side of the wedge which ha; § 
a constant slope (line BC). The side of the wedge js 
now assumed to end with a convex corner at point ( 
At the convex corner one obtains the sonic speed; th 
The floy 


over this corner will result in a Meyer expansion. The 


sonic line will therefore start from point C. 


boundary stream line will therefore map to a char 
(Considerations 0 
The 
amount of expansion occurring in the flow over th 


acteristic of the hodograph plane. 
this kind are carried out in detail in reference 2.) 


corner depends upon the change of direction taking 
place at this point. 

In general, only a certain portion of this expansion is 
of importance for the subsonic part of the flow pattem 
This can be seen as follows: In the previous example ofa 
flow with a free-stream Mach Number 1,! two classes 
of Mach waves that start at the body could be differ. 
entiated. The Mach waves of the first class end at the 
sonic line, while the Mach waves of the second class 
never reach the sonic line. If the contour of the body is 
changed such that the resulting perturbation propagates 
along one of the Mach waves of the first class, then 
this perturbation will reach the sonic line and thus in 
fluence the subsonic part of the field. A perturbation 
that travels along one of the Mach waves of the second 
class has no effect upon the subsonic part of the flow 
There exists one Mach wave that is the last to reach the 
sonic line and thus separates the above two classes; it 
forms the limit of the supersonic region necessary t 
We shall therefore de- 
note this wave as the ‘‘limiting Mach wave.” In the 
hodograph plane, Mach waves that start from the 


determine the subsonic field. 


body are represented by epicycloids, such as AD in Fig. 
lb. The Mach wave that reaches the sonic line at in 
finity 
AD since point A is the point of the sonic line which 
Accordingly, 


that is, the ‘limiting Mach wave’’—is given by 
maps to infinity in the physical plane. 
only the part of the supersonic boundary CD is of im- 
portance for the determination of the subsonic part of 
the field, although the expansion may go beyond point 
D. One might try to follow the continuation of the 
field downstream of the “limiting Mach wave’’ in the 
hodograph plane. Actually, this will not be necessary, 
since the continuation of the field is purely a super- 


sonic matter, and the computation may be carried J 


graph. 


out in the physical plane without the aid of the hodo- | 
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FLOW OVER A WEDGE 


The singularity at point A has been investigated in 
detail in reference 1. Accordingly, we obtain the fol 
lowing boundary value problem: We seek a_ solution 
of Eq. (2) in the shaded part of Fig. 1b which has a 
given singularity at point A and assumes the value 
= 0 along the boundary ABCD. Along AD no 
The type of boundary 


y 
conditions can be prescribed. 
value problem considered here is closely connected 
to the one dealt with by Tricomi.* Therefore, one 
sees that the solution of the problem exists and _ is 
unique. 

In Fig. 2a the hodograph plane is shown in terms of 
the variables » and 6. In the transonic law of similar 
itv, a limiting process occurs in which the thickness 
ratio of the body tends to zero. In Fig. 2a the changes 
of the location of the boundaries that occur in the 
course of such a limiting process have been counteracted 
by an enlargement of the scale. As a consequence of 
this scale transformation, the stagnation point B will, 
in the limiting case, lie at negative infinity. (For a 
more detailed representation of this limiting process 
the reader is referred to reference | or 2.) 


(4) AN EQUIVALENT BOUNDARY VALUE PROBLEM 


It will be seen later that it is relatively simple to ful 
fill the boundary conditions along the lines AB and 
BC. 
if one tries to fulfill the condition along CD. 
useful that the condition along CD can be transformed 
constant 


However, one encounters considerable difficulties 
Here, it is 
into conditions prescribed along the line 6 = 
which runs through point C—+.e., along CE of Fig. 
2b—and along a line » = constant 
We shall first represent the procedure in a somewhat 


generalized way. 


i.e., along ED. 


Let the prescribed condition along the characteristic 
CD of Fig. 2b be expressed by y = (0). Now if one 
imposes the arbitrary condition along the line CE y = 
io(n), then the flow field in the region DECD will be 
uniquely determined. It can be computed by means 
of the method of characteristics. The function y and 
its normal derivative y, along the line DE can especially 
They may be denoted by /3(@) and /,4(@), 
respectively. 


be found. 


Conversely then, the flow field in the region DECD is 
uniquely determined if the functions h(n), /3(0), and 
h.(6) are given; it can be found again, for example, by 
means of the method of characteristics. In particular, 
it will give the function /,(@) along the character- 
istic CD. Accordingly, if one has a function ¥(n, @) 
which fulfills the differential equation, Eq. (2), and 
which (1) assumes along the subsonic boundary the 
values required in the original problem, (2) assumes 
along CE the values ho(n), (3) assumes along DE the 
value h3(@), and (4) whose normal derivative assumes 
along DE the values h,(0), then this solution will assume 
the values of /,(8) along the characteristic CD. 
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Fic. 2b 


problem 
Equivalent boundary value problem 


Fic. 2a (top). Original boundary value 


(bottom). 


The solution of the boundary value problem modified 
in the above manner be the 
as the solution of the original boundary value prob- 


will therefore same 


lem. 


In the present case, the prescribed boundary condi- 
tion along the characteristic CD is y = /,(@) = 0, and 
the function f(y) is chosen as zero; therefore, the 
function W will vanish identically in the region DECD. 
Consequently, the prescribed conditions for y and y, 
along DE will be y = h;(0) = 0 and y, = h,(@) = O, 
respectively. 

(5) THE SINGULARITY AT Pornt A 

The desired singularity at point A can be expressed 
by means of the particular solutions obtained by the 
hypothesis (cf. reference 1) 


g= in "f(g, n) (9) 


where 


¢ = (9/4) (07/n’) (10 


and 7 is a constant. 


Inserting this expression into the differential Eq. (3 
will yield the following ordinary differential equation for 


IG, #): 
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| 
9 n(n — 1) f = 0 


(11) 

Eq. (11) is a hypergeometric equation; the properties 
of its solution are well known. The considerations for 
determining the value of the parameter n in the case of 
the flow with a free-stream Mach Number | are carried 
out in reference 1, and in brief they are as follows: 

(1) The ‘limiting Mach wave’ AD must not lie at 
infinity in the physical plane. Since the gradient of 
the Legendre potential is proportional to the space vec- 
tor in the physical plane [cf. Eq. (4)], one obtains the 
requirement that along the Mach wave AD the func- 
tions f(¢, m) and f’(¢, 2) must remain finite. For neg- 
ative values of » (which are required to produce the 
singularity at A), this restricts the choice of f(¢, ”) to 
one of the linearly independent solutions which is regu- 
lar at the point ¢ = 1. This point corresponds, accord- 
ing to Eq. (7), to the “limiting Mach wave.”’ 

(2) Along AB one has the condition y = 0. 
¢o is proportional to y according to Eq. (5), one has the 
condition that the transformed potential is symmetric 
to the line AB of the hodograph. A solution that ful- 
fills the above condition will, for an arbitrary value of 
n, in general be represented by a function along the 
line 6 = 0 which consists of a symmetric portion and an 
Setting the coefficient of the 


Since 


antisymmetric portion. 
antisymmetric portion equal to zero will yield a condi- 
tion for n which is fulfilled by an infinite number of dis- 
crete values of n. 

(3) Among the values of ” found above, only those 
can be used for which the value of y tends to infinity at 
point A and for which the representation in the physical 
plane is single-valued. 

From these considerations the suitable value of 1 is 
found to be nm = —1. Accordingly, the desired singu- 
larity is given by 


fF, — 1) (12) 


(6) OTHER PARTICULAR SOLUTIONS 
Another family of particular solutions of Eq. (3) may 
be found by the hypothesis 


g = g(n, m) cos (m7O/8) (13) 


where m is a positive integer and @ is an arbitrary angle. 
Inserting this into the differential Eq. (3) yields the or- 
dinary differential equation for g(n, m)—namely, 


gz” + (mm/0)°*ng = O (14) 
The solution of this equation is 
g = V7Z,,,[(2/3)(mr/Oo)n *| (15) 


where Z,,, represents any linear combination of Bessel 


functions of the order '/;.t For negative values of 


tT Cf. reference 4. 
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g(n, m) 
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n, this solution is represented by 


, 2 .\/mr 
g = Vind, (3 a Jim fp 7 <0 (I 


Bessel functions of an imaginary argument have th} 


asymptotic behavior of exponential functions 
one linearly independent solution tends to zero for lary 
values of the argument while the other tends to in 
finity. 

Since the stagnation point at the nose of the wedg 
corresponds to negative infinity in the 4, 6-plan 
only the solution that tends to zero as the argumer 
tends to negative infinity is permitted, inasmuch , 
the nose of the wedge must lie in the finite portion of th; 


physical plane. The desired particular solution js 


therefore, 
os a) [2 mr. 
g=1°V \nlily, |; tn) *), a< 0 (Nf 
5 A 
where //, ,‘" is the Hankel function of the first type 
the order '/;. The analytical continuation of this 


solution to positive values of 7 is given by 
2 
/5 Vn x 
V3 
mr 


E (| , ' 4g (; mmr , ) ~ 
at : F ant ’ h) 
‘ls g n \s 9, 0 n \ 


where J .../, are Bessel functions of the order + '/;. 


The derivatives of these functions are given by 


dg mr. yo (: ort 50 ) <0 (I 
= 1 . 2 v : as 
dn Ay . 3 9p . " 
and 
dg 2° mx 
=— n 
dn V K A 


vs 2M® 5; 
n ) — Jz, ( n )I n>O— (20 
3 % 


represents Hankel functions of the first 
are Bessel functions oi 


2mr 
| 
3 % 


where J/:,," 


type of the order 7/3; and J .:2,, 
the order +?/;. 

For large values of the argument, the following 
asymptotic representations are obtained (cf. reference 


4): 
34, V/9 n 4 
e 
m tg 
(=) ao (Co 
cos 
mM 1 


3 
0, one has 


( 


— (2/3) (mm/0)'n 2? n < 0) | 


2 Tv 
=> "), 1 >0 


(21 


g(n, mM) ~ 


For 7 = 


ar 2 (“") sa ‘ 
~~ t—-ifsary °°” 


And, finally, for the derivatives, one obtains 
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J om ~ (2/3) (ma /Oo)\n\3 72 
dt + (7) gr. ce 


an (22) 


. 4 y) 
dg en) er (; mr s), 1 *) _ 
; ( i: °° “ua” 


dn 
J 2 mn\' 
dg _ ah = ( ) , n=O 
dn (—2/3)! \ % 


From these expressions the following asymptotic 
formula, which is basic for the method of solution used 
subsequently, is found: 


2 MT «; T ; 
tan | - n° — \)" n>0O (23) 


> 
” 0 


dg/dn mt 


g A 


(7) REFORMULATION OF THE BOUNDARY VALUE 
PROBLEM IN TERMS OF THE LEGENDRE POTENTIAL 


In Sections (3) and (4) the boundary value problem 
was formulated in terms of the stream function. In 
this form the boundary value problem conforms to 
the type considered by Tricomi.* From the stand- 
point of the actual computations, however, the use of 
the transformed potential is preferable, because of the 
fact that the coordinates in the physical plane can be 
obtained from the transformed potential merely by 
partial differentiations whereas an integration is neces- 
sary in the case of the stream function. The corre- 
spondence between the two may be obtained from 
Eq. (5). One merely replaces in the boundary condi- 
tions the function y by (1/p*)g,; then 7 is replaced by 
(1/p*)¢,9. Since the latter expression is prescribed only 
along a line » = constant—.e., along the line DE of 
Fig. 2b—one can integrate it with respect to @ and ob- 


tain 
¢, = S G0 + const. 


Consequently, ¢, is prescribed along DE except for a 
constant. Since the transformed potential is symme- 
tric with respect to the horizontal axis, the same con- 
stant would occur in the lower half of the hodograph. 
A solution of the hodograph equation that has g, = 0 
everywhere and for which ¢, = constant along DE is 
given by 
yg = const. 7 

Therefore, the boundary conditions for ¢ are deter- 
mined from those for y with the exception of a term re- 
sulting from a solution given by the last equation. 
The presence of such a solution would change the x- 
coordinates by a constant; in other words, it represents 
a shift of the system of coordinates in the physical plane. 
Naturally, such a shift is unessential. 


Accordingly, the boundary value problem in terms of 
¢ has the following prescribed conditions: g = 0 
along ABCE and ED (Fig. 2b) and ¢, = 0 along 


ED. 
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(8) PROCEDURE FOR THE FULFILLMENT OF THE 
BOUNDARY CONDITIONS 


The computations start with an expression that 
yields the desired singularity at the point that corrre- 
sponds to the free-stream conditions. This expression 
does not fulfill the boundary conditions. However, 
since the condition g, = 0 along the line BCE can be 
considered as a condition of symmetry, it may be ful- 
filled rather easily by an image system. Then, since 
¢» is both continuous and antisymmetric along AB, 
the condition g, = 0 there will be automatically satis- 
fied. 

In Eq. (15), if one sets 4) equal to half of the angle 
at the tip of the wedge, then yg, formed with these par- 
ticular solutions will be zero for 6 = 0) and @ = 0—that 
is, along the surface of the body and along the line AB. 
These solutions can be arbitrarily superimposed with- 
out influencing the boundary conditions along BCE and 
AB and will, therefore, be used to fulfill the boundary 
conditions along ED. 

In the present case, the fulfillment of these condi- 
tions will be carried out in two different steps. (1) 
One forms an expression that fulfills the condition 
for ¢g» without any regard to the conditions prescribed 
for ¢,. (2) One adds a solution that does not con- 
tribute to the value of g, along ED but fulfills the pre- 
scribed conditions for ¢,. 

Accordingly, the transformed potential will be ob- 
tained as the sum of three portions—that is, 

g=¢gt at ¢s (24) 

Here ¢; is the portion that originates from the image 
system—that is, it arises by a superposition of an in- 
finite number of functions given by Eq. (12). The 
singularities are located along the 6-axis at the interval 
of 26. 

¢2 arises by a superposition of the particular solutions 
given by Eq. (13) such that the condition g, = 0 along 
ED is fulfilled. (¢g. is not uniquely determined by this 
requirement as will be seen later. ) 

¢3 arises from a superposition of functions given by 
Eq. (13) such that along ED (or, in other words, 7 = 
no and 6/2 < 6< 4) one has ¢;, = 0 and $1, 4 $2, aa 


¢3, = 0. If ge is chosen, then ¢3 will be determined 
uniquely. 


(9) COMPUTATION OF ¢; 


If one differentiates Eq. (12) with respect to 6, then 
one obtains 


9 


go = —3ini OV elf’, sign 9-0 


where the symbol sign 7-6 denotes that the preceded 
quantity is to be multiplied by the sign of the product 
n:°6. Accordingly, ¢» is an odd function with respect to 
the n-axis. If one now has a series of these singularities 
located at an interval of 20) along the 6-axis, then 


along a line @ = constant which lies in the middle be- 
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tween two singularities—e.g., at 6 = 6 )—the con- 
tribution of g, from the singularities in the upper half 
will be canceled by that from the lower half-plane. In 
other words, the boundary conditions will be fulfilled. 
Consequently, ¢; is given by 


ga =—|gi-' dS fltm —1) (25) 
where 


Gm = (9/4) [(8 — 2mO%)?/n*] 


The partial derivatives of Eq. (25) with respect to 
n and @ are 


go =i|n-? Dd (Binf’ +f), signn (26) 


ip =e » 2 Vv Cm tA sign v- 2m6o)n 


m 


(27) 
respectively. 
Let us check the convergence of these expressions. 
For large values of the argument, the function f(¢,,, —1) 
is represented by 


* sas Cin ss Com 


where c,; and c. are constants (see reference 1). Ac- 
cordingly, the particular solutions (12) for large values 
of (, are given by 


9 2/9 » $ 
1, ~ o1 E ae 2m) | + Cn E ie — 2m) 
” oO 


where Pip, is the mth term of the series in Eq. (25). 
Hence, one sees that for large values of ¢,,—1.e., for large 
values of m—the mth term of the series in Eqs. (26) and 
(27) will be less than, or equal to, c’m~** or c"m~"*, 
respectively, where c’ and c” are constants independent 
of m. Consequently, the convergence of the partial 
derivatives given by Eqs. (26) and (27) will follow. 


(10) COMPUTATION OF ¢» 


Let the value of along the line ED be denoted by 
m. From the equation for the characteristics, Eq. (7), 
one finds 

m = [(3/4)00.]” (28) 

If one now inserts this expression into the asymptotic 
formulas given by Eqs. (21) and (22), one finds that 
for » = m the absolute values of the functions g(y, 77) 
and g’(n, m) will be larger than a finite constant for all 
values of m. Therefore, the family of solutions given 
by Eq. (15) will furnish a complete set of functions sin 
m0/0) and cos ni7@/ 0) in the interval 0 < @< @. 

If one considers solely the fulfillment of the condi 
tions for gs, then one has at his disposal an entire set 
of orthogonal functions defined in the interval 0 < 


6 < 6. However, the condition for gs along ED e 
pressed by the requirement that $1, + ox, = Ois giver 
in the interval 4/2 < @< @ only. Obviously, one cay 
choose the continuation of ¢), arbitrarily in the intery, 


from 0 = 0 to @ = 6/2; then the determination of g, yy 





be unique. The arbitrary element introduced in this [ 


way will vanish if one determines ¢3. For our coy 


putations, we choose $2, in the interval 0 < 6 < 4: 


symmetric with respect to the line @ = 6/2 with th 
values in the interval 0/2 < 0< @. 

The actual computation is carried out in the following 
manner: One first finds the values of ¢), for » = », 
and 6)/2 < 6< @, and then continues the curve into th, 
interval 0 < @< @/2 as described above. This function 
is then represented as the Fourier sine series 

(2m + 1)26 


») —S/s x 
»o 
Ply = —3 ( 0) y Bom+1 sin (29 
} m 0 Ay 
where #B,,,,; are constants determined by the Fourier 
analysis. Next, ge is written in a corresponding form 
g(n, 2m + 1) (2m + 1)76 
g2 = >» Aom +1 9 cos 
m 0 2(n, 2m + 1) A 
Here, the coefficients A»,, ,, must be determined from 
the condition that ¢1, + ¢,, = 0. If one now differen- 
tiates Eq. (30) with respect to @ and equates the te 
sulting expression to Eq. (29), one finds the coefficients 


Aom ;,aS 


¢ ) * Ay Bom +1 
FO 3, = —3 A 
4 x 2m-+ 1 


Accordingly, the partial derivatives g2, and ¢»_ will be 
b 4 


3 ; ; 3(n, 2m + 1) 
P25 = 3 (* 4) > Bom +1 = 
m 0 


2(mo, 2m + 1) 
(2m + 1)70 


sin (3 
A 


3 7 ae 
b = —3 0) _ x 
. ” (; T ey 2m + l 


g’(n, 2m + 1) (2m + 1)76 
cos (02 


2Z(m, 2m + 1) A 


Since, according to the asymptotic representations 
Eqs. (21) and (22), the order of magnitude of the func- 
tions g(n, 77) (as well as for the functions g’) is always the 
same for all values of » > 0 with the exception of the 


4 


universal factor n ‘, the convergence of the function 
¢2, (as well as for ¢2,) is assured if the corresponding 
series converge at n = m. Since the latter series are 
Fourier expansions of closed functions, the expression 
for ¢g, and g: , Eqs. (31) and (32), will converge for all 


positive values of n. For negative values of n, the con- 


vergence is assured because of the presence of the ex- 


ponential function [ef. Eqs. (21) and (22)]. The fac- 
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‘ occurring in this consideration seems to indi- 


tor 7 . 
i.e., for 7 = 0. 


cate a difficulty at the sonic velocity 
Here, the functions g and g’ can be determined directly 
and are given in Eqs. (21) and (22). One finds that the 
mth term in the expression for ge, at 7 = 0 is compar- 

- 1 6 
able to the mth term at 7 = mo except for a factor mm’, 
whereas at the mth term for ga factor m must be 
added to the corresponding term at 7 = %. Accord- 
ingly, the convergence of $2, for » = O is somewhat 
better than for » = m and that for ¢g», is decreased as 
compared to the expression at 7 = m. However, since 
even the Fourier representation of the function gg at 
ny would converge, the latter factor will be harm- 


n= 
less. 


(11) CALCULATION OF ¢3 


The superposition of the functions ¢; and ¢g» yields an 
expression that has the desired singularity and which 
satisfies the condition that g, is zero along the boundary 
ABED. In this regard the boundary conditions are 
fulfilled. However, the condition ¢, = 0 required along 
DE has not yet been considered. Therefore, a par- 
ticular solution ¢; will be added to the previous ones 
which has ¢3, = 0 along the boundaries ABE and ED and 
assumes along ED such a value of $3, that gi, + 2, os 
g, = 0. The determination of such a function will, in 
general, lead to an infinite system of equations. We 
shall present here a method based on the asymptotic 
representations of the particular solutions which pro- 
vides a simpler method of attack. The function ¢3 
is represented by a superposition of particular solutions 


given by Eq. (13)—that is, 


mr 
‘OS (23:3) 
1 Mm £(no, m) 0 


A bm g(n, m) 
_ Cc 


ao 

T m 
The conditions along BCE will 
Along ED, one ob- 


where b,, are constants. 
then be automatically fulfilled. 
tains the conditions 


where F(@) has been determined previously so that it 
can be considered as a known function. If one inserts 
Eq. (33) into the above boundary conditions, one ob- 


tains 


> by sin (m7O/0)) = 0 
1 


m 


9% bm fg'(n, m) mrO 
< » cos 
Tm=1mM 2( n, ™) n No A 


(8/2 < 6< H%) 


_= F(@) 
(34) 
If, in the second equation of Eqs. (34), the asymptotic 


expression for g’/g for 7 = m is inserted from Eq. (23), 
then it is simplified to 
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mn 


0) ~ (-1)" bp cos = 


(00/2 < 6 pa 60) 


F(@) 
(35) 


The error, due to the use of the asymptotic formula, is 
3.9 per cent for m = 1 and decreases as m becomes 
larger. The first condition of Eqs. (34) can be fulfilled 
in the following way: The function ¢s3, is zero in the 
interval 0)/2 < @< 6, while in the interval 0 < @< 6/2 
it will, in general, be different from zero. Let the value 
of ¢3, in the latter interval be denoted by p(@). Then 
the coefficients 5,, are expressed by means of a Fourier 


analysis of the function $3, i.e., one obtains 


») 60/2 
- . mru 

= p(u) sin du 
A 0 A 


If the coefficients 


Dus (36) 
where wu is a variable of integration. 
are determined according to the last equation, the first 
condition in Eqs. (34) will be fulfilled for any choice of 
the function p(w). If one inserts the expression for },, 
as given by the above equation into Eq. (35), then one 


obtains 
2 (3 ; ws 
- ( a) > (- 1)” / p(u) X 
A 4 m ] « 
mmru m0 ” A . 
sin du cos = F(@), = <0< 6 (37) 
0 0 - 
This is an integral equation for the function p(w). The 


previous problem of finding the coefficients 5,, has been 
transformed into that of finding the function p(x). 
It should be noted that in the above equations the 
range of 6 is between 0/2 < 6 < 6, while u varies from 
0 to #/2. In order to have the same interval of the var- 
iables, one introduces the following transformation: 


Eq. (36) with the above transformation becomes 


aia Ea . 0/2 
al , % » » p(u) X 
- ” 1 
mTrUu mTrv 


du cos = 
0 0 


sin Fi(V) (38) 


where 0 < u< @/2and0< v < 6/2. 





AY 


~~ R= 0 ; 
° WZ = FY) 
> —~ = 
£2 < A? | 
—00 ~<— 4 B= O 8 x 


Fic. 3. 














Boundary value problem of the incompressible analogy. 
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In the incompressible flow theory, there exists a prob- 
lem that yields exactly the same equation as the present 
problem—namely, Eq. (38). A discussion of this in- 
compressible example will be instructive for the subse- 
quent analysis. 

Consider a boundary as given in Fig. 3 where the 
horizontal axis is denoted by x and the vertical axis by 
y. The differential equation to be fulfilled in the 
shaded region of Fig. 3 may be that for the incompres- 


sible fluid—namely, the Laplace equation 
¢rx + Gy = 0 


The boundary conditions prescribed are that yg, = 0 
along ABC and along DE, and ¢, is a given function 
F\(y) along CD. If one replaces gy, by a new function 
¢@ and maps the region ABCDE onto a half-plane, the 
above boundary value problem can be interpreted as 
that for finding the lift distribution over a thin airfoil. 

A harmonic function that fulfills the boundary con- 
ditions along AB and DE is given by 


BA" 2 
oe ( Go } 
4 ma = 
where the 3,,’s are constants. 
zero along BC and will assume some value p(y) (un- 


known as yet) along CD, the coefficients b,, above can 
then be determined by means of a Fourier analysis; 


the result is 
9 8/2 ‘ / 
2 _ mry 
= — / p(y’) sin dy’ 
6 0 A 


where y’ is a variable of integration. 
If one now expresses the condition given for ¢, along 
CD, one obtains 


ye mry 
cos 
1m 0 


max /Oo 
é 


g(x, ¥) = 


Since the function g, is 


bm 


mry’ 


3 (° ‘ ) - = by /2 
= 0 
2\4 m= 1 0 A 


p(y’) sin 


’ mry 
dy’ cos = 


0 


Fi(y),O< y< 6/2 (39) 
This is exactly the same equation as the one found 
previously for the main problem—that is, Eq. (38). 
[To be exact, the variables u and v of Eq. (38) must be 
replaced by y’ and y, respectively. | 


The Solution of the Incompressible Analogy 


As mentioned above, it is clear that the solution of the 
incompressible problem can be found by a suitable map- 
By means of this method one will finally 
We shall 


give here a method of derivation of the result which is 


ping process. 
obtain p(y’) expressed in terms of F;(y). 


not rigorous from a strictly mathematical point of 
view; however, it gives the desired expression for p(y’) 
directly in a convenient way. A more rigorous deriva- 
tion of the result can be obtained either by following 
the mapping process or by using other properties of the 
incompressible flow theory. 
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Eq. (39) will first be rewritten with the order of the 
integral and the summation signs interchanged—thg 


3 (; ) : 7 my’ 
_ A b(y’) sin : 
2\4 : 0 PU p> 1 Oo 


mry |, 
cos dy’ = Fi(y) (3% 
0 


One now introduces 





in this equation the closed ¢. 
pression given by 


emits conaabeie 


my’ mry (1/2) sin (ay’/) 


> sin cos = 


“= % 0 o cos (ry/A) — cos (ry'/, 


(4() 


(The use of this correspondence in the integral equatio, 
can be justified by carrying out the mapping proces 
mentioned above.) Eq. (39a) with the use of Eq. (40 


becomes 


ne -/ - b(y’) XX 
‘G ) I PC 


sin (ry’ /o)dy’ , 4s 
— : = Fi(y),0O< y< &/2 (4 
cos (ry/A) — cos (ary"/6) 

The singular kernel of this equation offers no essential 
difficulties, being similar in nature to that in the down. 
wash integral of incompressible theory. 


The integral Eq. (41) will now be re-expressed in terms 

of the new variables 7 and s defined by the equations 
, 
Ty 

cos = 
0 


l 
5 (lL + cost), ply’) = pilt) 
7 (42 


cos— =-—(1+ oss), Fi(y) = Fi(s) 
A z 


The integral equation in terms of these new variables is) 


17/3 a " sin ¢ dt 
= A pi(t) = 
wr \4 0 cos s — cos ft 


(O< s<m) (4 


F,(s) 


In order to solve this integral equation, one first intro- 
duces the correspondence in Eq. (40) (expressed in terms } 
of s and ¢ instead of y and y’) into Eq. (43); the results 


2/3,.\" 
— ( m) / pi(t) >> sin mt cos ms dt = Fils 
Tv 4 x 0 m 1 


(O<s<-m) (4)! 


In this form the soluticn is rather obvious. First, one | 
expresses F,(s) in terms of a Fourier cosine series in s— | 


that is, 


F,( $s) = 


7. dm COS MS 


m 0 


where 


T 
dm = (2 nf F.(s’) cos ms’ ds’ 
0 


Here, s’ is a variable of integration. 
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It was mentioned in Section (7) that the expression 
_can be changed by a constant without an essential 
If the given 


a 
. 


effect upon the flow in the physical plane. 
values of ¢, are changed by a constant, then F2(s) will 
also be changed by a constant amount. Accordingly, 
the constant in the expression for F:(s) may be dis- 
regarded. 

One now makes the hypothesis 


p(t) = do Cy sin mt (46) 
m 1 

Inserting Eqs. (45) and (46) in Eq. (44) and using the 

orthogonality relationship for the sines, one obtains a 

relation between the unknown coefficients C,, and the 

known coefficients d,. has C,, dé. 

With this result the function p;(¢) can be expressed in 


One simply 
the form 
» x "x . 
nd=- dL / F,(s’) cos ms’ ds’ sin mt (47) 
T m i «© 
or expressed in the closed form 


/ F3(s") 


where the correspondence in Eq. (40) has been used. 


ds’ 


/ 
cos s — cost 


sin ¢ 
pill) = (48) 


us 


Eq. (47) represents the solution of Eq. (39). In these 
equations no serious difficulties of convergence will 
Eq. (47) has been found by means of the in- 
However, it must be emphasized 


occur. 
compressible flow. 
that the validity of Eq. (47) is naturally independent of 
the manner in which it has been obtained; it can be used 
without reference to a specific example. One might 
question the uniqueness of the solution. In the case 
of the incompressible analogy and in the case of our ex- 
ample, the uniqueness is obvious from general unique- 
ness theorems concerning the flow problems (in the 
case of the transonic flow from Tricomi’s uniqueness 
theorem’). 


Application of the Incompressible Analogy to the Present 
Problem 


The results of the previous paragraph can now be 
applied to the main problem. It was seen previously 
that the equation in this case [Eq. (38)] was identical 
to the above incompressible analogy if the variables y 
and y’ are replaced by v and u, respectively. Accord- 
ingly, Eqs. (47) and (48) may be applied directly to the 
present problem, provided in the transformation Eq. 
(42) the variables v and w are inserted in the place of y 
and y’, respectively. 


3) (36), together with 


either of the above results, comprise the final results 


Consequently, Eqs. (33) and 


lor gs. In the actual computation, however, the func 
tion p(t) in the form of Eq. (47) was used. The con 
vergence of the function g; and its partial derivatives 
with respect to n and @ can be investigated by following 
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| / 
9 
a 
Fic. 4. gi, as a function of @ for n = 7 
the influence of the different transformations. It can 


be shown that ¢; and its derivatives will converge. 


(12) DETAILS OF THE COMPUTATIONS 


The starting point of the computations is the evalua- 
that is, for 7 


tion of ¢, along the line DE ny and 


A 2 A 0 ‘< 0. 
tion of the function f(¢,,, —1) and its derivatives are 


Here, the formulas for the determina- 
given in reference |. In some cases more accurate val- 
ues were determined by means of series representa- 
tions of these functions. The resulting values of ¢), are 
shown in Fig. 4. 

In order to determine go, this curve was continued 
symmetrically with respect to the line 6 = @)/2 (dashed 
1). 
; are found by means of a harmonic an- 


line of Fig. From this curve, the Fourier coeffi- 
cients Bo, 
alyzer; they are listed in Table 1. 

In order to determine ¢;, one first forms the expres- 
sion F(@) 


In F(@) the transformation from the variable 6, which is 


— (¢1, + ¢») by means of Eqs. (26) and (32). 


essentially the variable v, is now made to the variable 
s. The relation between v and s (as well as between 
uw and ¢) is given by Eq. (42), in which y and y’ are re- 
placed by v and u, and is shown in Fig. 5. A Fourier 
analysis is now carried out for the function F2(s) (cf. 
Fig. 6), and the coefficients d,, of Eq. (45) are com- 
puted. From this result the expression /;(/) as given 
by Eq. (46) may be found, since C,, = d,,. Because of 
the singularity that occurs in the relation between 
A = 6/2, the function p(@) will 
have a vertical tangent at 0 
coefficients b,, for Eq. (33) by a Fourier analysis of the 


6 and s for v 


7 = 
7 = 


& 2. One now finds the 


ABLE | 
Values of Be, ,, for Eqs. (81) and (32 
m B.,, m B 
0 —0 OS875 5 0.0011 
l 0.0124 6 0. 00059 
= —() 0042 7 0 0050 
: 0 O0OR3 & —0 OQOO490 


} 0.0012 
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Fic. 5. Transformation of the variables v and u to s and f. 





3 
" 


Fy(s) forn = nw andO<cs <r 


Fic. 6. 


function that is equal to p(@) for 0 < @< /2 and equal 
to zero in the interval 0)/2 < 06< @. The coefficients 
are given in Table 2. The vertical tangent mentioned 
above for the function p(@) at 8 = 6/2 will cause this 
series to converge more slowly. 

As a check for the solution 3, the fulfillment of the 
original boundary condition that y, = 0 along the 
characteristic CD will be verified. Fig. 7 shows the 
values of Pig» Ply + Hrs and ¢1, + $2, + ¢3, as functions 
of @ along the characteristic CD. It can be seen that 
this boundary condition is satisfactorily fulfilled. 

The next step of the computations represents the 


transformation of the results back to the physical 


THE AERONAUTICAL 


SCIENCES—NOVEMBER, 1950 


TABLE 2 
Values of b» for Eq. (33) 


3.\° 3 VM 
bm 0 m } 
m . ( . ) m b ( w) 





1 0.01796 11 0.00069 
2 0.01346 12 — 0.00251 
3 —0.00354 13 — (0.00072 
4 —(0.00858 14 0.00200 
5 0.00061 15 0.00069 
6 0.00569 16 —0.00164 
7 0.00034 17 —(0.00052 
8 —(0.00412 18 0.00144 
g —0.00041 19 0.00050 
10 0.00320 20 0.00118 


plane. It is carried out by means of Eqs. (4). \ 
shall restrict ourselves to the determination of the 
ordinates of the profile and the corresponding pressy 
distribution. The x-coordinates as a function of 9 a 
found according to Eqs. (4) with 6 = 4. The pressy 
coefficient is then found as a function of n by means; 
Eq. (6). 
in Fig. 8. 


The resulting pressure distribution is show 


In the original simplifications of the basic equation 
by means of the transonic law of similarity, it was a 
sumed that the entire velocity field was in the neigh 
borhood of the critical velocity. The stagnation regia: 
at the nose of the profile will, therefore, be out of th 
range of validity of the present theory, just as in th 
case of the Prandtl-Glauert approximations. However 
Fig. 8 shows that this region is indeed small for con 
ventional thickness ratios of the profile. 
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Distribution of ¢, along the characteristic CD. 
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Pressure distribution over the nose. 
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to right.) 


; Construction of the flow from the sonic line to the 
limiting Mach wave.’ (The free-stream direction is from left 
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Fic. 10. Construction of the rear part of minimum drag. (The 
free-stream direction is from left to right.) 
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Fic. 11. Pressure distribution for the profiles of Fig. 12 with 
different positions of the maximum thickness. 





PROFILE B WITH REAR PART OF MINIMUM DRAG 
Fic. 12. Profiles for the pressure distributions of Fig. 11. 
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Fic. 13. Comparison of the pressure distributions over wedge 


and cusped nose profiles of 10 per cent thickness. 
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Fic. 14. Construction of the straight-sided rear part. (The 


free-stream direction is from left to right.) 


The position of the sonic line is then computed, and 
the flow downstream to the “‘limiting Mach wave’’ [cf. 
Section (3)] is constructed by means of the method of 
characteristics. The result is shown in Fig. 9. 

The flow downstream of the “limiting Mach wave’”’ is 
arbitrary and can be constructed by the method of 
characteristics if the continuation of the zero stream 
line is prescribed. First, we shall assume that the rear 
portion of the body is such that it gives a minimum 
contribution to the drag at a given length of the pro- 
file. 


condition is fulfilled if one first allows a Meyer expan- 


In reference 1, it was shown that this optimum 


sion and then forms the rear part of the profile such 
that the compression waves coming from the sonic line 


are not reflected at the body. The length of the rear 


NOVEMBER, 1950 

portion of the profile is determined by the amount gj 
the Meyer expansion allowed at the corner. The gop. 
tour of the body downstream of the corner is the, 
found by integrating the slope of the stream lines fo, 
y = 0. The construction of the flow for a profile with 
a rear part of minimum drag is shown in Fig. 10. Jy 
Fig. 11 is shown the pressure distributions for profile 
obtained by using different amounts of expansions 4 
the corner; here, the different profiles are reduced to the 
same thickness ratio by means of the transonic law 9 
similarity. 
files in which the position of maximum thickness varies 


By this process one obtains a family of pro. 


The profiles corresponding to the pressure distribution 
of Fig. 11 are shown in Fig. 12. 

In Fig. 13, a comparison is shown between the pres. 
sure distribution over a wedge-nosed profile with that 
over a profile with a cusped nose as determined jn 
reference 1; the rear parts in each case are those for 
minimum drag conditions. 

The second type of rear portion is chosen so as to ob- 
tain the standard diamond profile with the maximum 
thickness located at the mid-chord. The construction 
of the flow in this case is shown in Fig. 14, and the 
corresponding pressure distribution is given in Fig. 1], 
It has been shown! that at a corner such as that occur. 
ring at the mid-chord of the diamond profile, a shock 
wave must arise; the shock in this case is formed by the 
compression waves coming from the sonic line which 
coalesce after being reflected from the side of the pro- 
file. In the present case, this shock was so weak that 
the mesh interval used for the construction was too 
coarse to show its presence. 
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Variation of Cp with the position of maximum thickness 
(thickness ratio 10 per cent). 
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Finally, in Fig. 15 the drag coefficients for the wedge- 
shaped nose with rear parts of minimum drag are given 
as a function of the position of maximum thickness. 
Also shown are the values for the diamond shape pro- 
file and for the bodies of reference 1. It should be 
noticed that the drags for the profiles with the cusped 
nose are lower than those with a wedge nose. Accord- 
ing to physical plausibility, the cusp at the nose must 
be considered as unfavorable. Naturally, the con- 
tribution to the drag that is caused by the tip of the 
profile in such a case is rather small; therefore, if the 
tip portion of the cusped profile is omitted, the body 
as a whole for a given length of the profile will appear 
to be blunter than a body with the same thickness 
ratio but without a cusp at the nose. In this respect 
the above results are astonishing. On the other hand, 
one sees from Fig. 13 that along the subsonic side of 
the profile the expansion is much more rapid in the 
case of the profile with the cusped nose than in the 
case of the wedge-nose profile. In the former case one 
has even a suction on the portion of the profile up- 
stream of the maximum thickness. 


(13) FINAL REMARKS 


In the present paper, the flow over a wedge-shaped 
profile at Mach Number | has been computed. The 
direct procedure for satisfying the boundary conditions 
for this flow in the hodograph plane leads to an infinite 
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system of linear equations. By reformulating the 
supersonic portion of the boundary value problem, it 
was shown that one could avoid solving this infinite 
system directly and reduce the problem to that of find- 
ing suitable Fourier expansions of the boundary con- 


ditions. 


The final results indicate that one may expect a rea- 
sonable value (as compared to the values that are usu- 
ally predicted) of the drag coefficient at Mach Number 
1. The present results are based upon the assumption 
that the flow is potential, so that the resulting drag 
coefficients must be considered only as an optimum re- 
sult. 
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SUMMARY s = distance around the perimeter of the box f 


Curves are obtained which show the total increase in weight of 
an optimum-designed sheet-stiffener compression panel when 
the skin thickness is increased by an arbitrary amount. These 
example curves are based on data from tests conducted by the 
N.A.C.A. on 75S-T Y-stiffened panels. However, the method 
outlined may also be applied to the plotting of other sheet-stiffener 
compression panel data. 

The N.A.C.A. data, when put in the form outlined here, indi- 
cate the minimum amount by which the total weight of sheet 
and stiffeners would be increased in order to achieve a given in 
crease in sheet thickness. If optimum design principles are fol- 
lowed, this increase in total weight is considerably less than the 
increase that would result if the amount of stiffener material 
were left unchanged. 

The method outlined here for presenting data is of particular 
use in preliminary design work where the effects of aeroelasticity 
requirements must be evaluated in terms of increased wing 
weight. 


INTRODUCTION 


O™ OF THE PROBLEMS encountered in the structural 
design of the wing for high-speed airplanes is that 
of providing sufficient structure to prevent aeroelastic 
effects, such as loss in aileron effectiveness, from ex- 
ceeding acceptable limits. Most of the aeroelastic 
effects experienced by the unswept wing are due to a 
reduction in effective torsional stiffness which accom- 
panies an increase in flight speed. For a wing having 
the flexural axis located aft of the line of aerodynamic 
centers, there is not a great deal that can be done to pre- 
vent this apparent loss in stiffness. What the designer 
usually does is provide sufficient structural stiffness so 
that at the speed in question the wing retains sufficient 
effective stiffness to keep the aeroelastic effects within 
specified limits. 

One method of increasing torsional stiffness, without 
redesigning the entire wing structure, is to increase the 
skin thickness of the structural box. Assuming that 
the shear flow is constant around the structural box, the 
expression for the section torsional stiffness, GJ (in.Ibs. 
per rad. per in.), is 


GJ = 4A/G/fds/t, (1) 
where 
A, = enclosed area of the structural box 
G = modulus of elasticity in shear 
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( 


t, = skin thickness 


Writing the denominator in series form, 


4A 2G 
Gis : 
S} 4 So 4 53 4 ‘ Sn (2 
he "oe es 


It can be seen from Eq. (2) that, since the terms iy 
the denominator are divided by skin thickness, the 
most efficient} way to increase the torsional stiffness 
at any spanwise station is to make the skin thickness 
uniform around the box. After this has been done, the 
section torsional stiffness may be determined from the 


expression 
GJ = 4A,°Gt,/s 3 
where 
t, = skin thickness (assumed to be uniform around 


the structural box) 


Ss = perimeter of the structural box 


Eq. (3) shows that, for a uniform skin thickness around 
the box, the torsional stiffness varies directly with the 
skin thickness; hence, any further increase in stiffness 
is most efficiently obtained by a uniform increase in 
skin thickness. 

At this point it is desirable to know how the weight 
of a sheet-stiffener compression panel is affected by an 
increase in skin thickness, the loading remaining un- 
changed. It is incorrect to assume that an increase in 
torsional stiffness would require an increase in skil 
thickness without any decrease in the amount of stiff- 
ener material. Actually, as the skin thickness is in- 
creased, the amount of stiffener material may be de- 
creased, since the axial loading to be carried remains 
constant. This, of course, can be done only in the ear- 
lier stages of design. 

By using data from tests made by the N.A.C.A., it 
is possible to obtain curves showing the total increase in 
weight caused by increasing the skin thickness by an 
arbitrary amount. 


METHOD OF OBTAINING CURVES 


The data on which the curves of this report are based 
were obtained from Fig. 20 of N.A.C.A. Technical 


+ ‘‘Efficient”’ is used here in the sense of obtaining the max!- 
mum increase in stiffness for the amount of weight added. 
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Note No. 1389, which contains data on tests of sheet- 


stiffener compression panels. Fig. 20 shows a plot of 


g/t vs. q/t, for lines of constant g/Lo, where 
q = average load per inch width (lbs. per in.) = 
total load divided by total width 
i = average thickness of the panel (in.) = total 


area divided by total width 
t, = skin thickness (in.) 
ly = effective pin-end column length (in.) 


Since the curves were composed of a series of arcs, 
envelopes were drawn to represent a continuous vari- 
ation. These envelope curves are shown in Fig. 1. 
The term q// represents the average compressive 
The term q/Lo is the 
structural index,’ which is a measure of the loading in- 


stress developed in the panel. 
tensity. The parameter g/t,, the average load per inch 
width divided by the skin thickness, is a convenient 
means for introducing the skin thickness as a parameter. 

A straight line through the origin (Fig. 1) represents 
a line of constant value of the ratio /,/?, the slope being 
equal to /,/7. By cross-plotting the curves of Fig. 1 as 
g/t vs. q/Lo for lines of constant ¢,/7, we can show how 
changing the ratio of skin thickness to average thick- 
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ness affects the allowable stress for any value of stry. 


turalindex. This cross-plot is shown in Fig. 2. 


It is apparent from Fig. 2 that, as the value of 1.) 


increased, the allowable stress is decreased. 
eurve, for ¢,// = 1.0, represents the limiting cag 
wherein all of the material is in the skin. This cyry, 


was calculated from Eq. (15b) in reference 2. This ¢ 


pression for the allowable stress F, of a flat plate use 


as an Euler column is 


F. = 0.97(Er)'" (q/Le)” (4 
where 
FE = modulus of elasticity (elastic range) 
t = E,/E, where /, = tangent modulus 
The position of this curve indicates that, as /, L ay 


proaches unity, the allowable stress drops off rapidly 


At higher values of q/Zo, it is likely that the rate of 


decrease of F, with increases in ¢,/f will be less pro. 


nounced. At the time of writing, however, there wer 
no experimental data available which covered values o/ 
q/Lo above a value of 800; therefore, this could not ly 


shown by the curves. 
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It can be seen that the curve representing optimum The average thickness may be written 
(minimum-weight) design may be attained only by us- : } : 
t= C.Lo + (q F,.) (6) 


ing an extremely low value of ¢,/7.. This indicates that 
such a value of allowable stress is not likely to be at- 
tained in compression panels designed for use in wing 
structures. For this reason, it is desirable to qualify 
allowable-stress curves, of the type shown in Fig. 2, 
with values of ¢,/7. The ratio /,/? is a useful parameter 
for estimating the relative amount of torsional stiffness 
which would be contributed by the panel when used in a 
wing structure. Its use introduces a further conven- 
ience because, according to the design charts of refer- 
ence 1, for a specified value of /,// there exists a variety 
of combinations of stiffener size and spacing which 
will yield very nearly the same allowable stress as the 
maximum for that value of ¢,/7 shown in Fig. 2. 

As shown by Shanley,’ the average thickness may be 
determined from curves of the type shown in Fig. 2 
and approximated by a straight line as follows: If the 
value of g/Ly at each point on the curves is divided by 
the corresponding value of qg/?7, we can obtain curves 
These curves can be 
For the set 


of ?/Lo vs. g/Lo as shown in Fig. 3. 
closely approximated by straight lines. 
of curves in Fig. 3, the straight lines are not distinguish- 
able from the curves except at extremely low values of 


q/Lo. The expression for the line is given by 
/Ly = Co + [(q/Lo)/Fe,] (5) 
where 
C, = intercept value at g/Lyo = 0 
F,. = effective stress = 1/slope of the line 


For any value of t,/7, the total average thickness / re- 
quired at any value of g/L» may be found from this 
expression. For a value of f,/7 = 0.6, the expression 
for the average thickness is 


? = 0.00155Ly + (¢/83,000) 


For example, if a panel 20 in. in length and having pin- 
end conditions is subjected to a loading of 300 lbs. per 
in., the total average thickness required is 0.0246 in., of 
which the skin thickness is 0.0148 in. It should be re- 
membered that these numbers refer to panels made of 
75S-T Alclad skin 75S-T extruded Y-section 
stiffeners. 

Fig. 2 shows how a change in the ratio /,// affects the 
However, it does not show 


and 


allowable average stress. 
how the average thickness is affected by changes in 
skin thickness, the loading remaining constant. This 
can be done as follows: 

Using Fig. 1 and choosing some value of f,/] to be 
used as a reference value (designated by f,,//o), we can 
obtain curves of ¢,/t,, vs. //fo for each line of constant 
g/L», where 

t, = skin thickness of a panel having a value of 

t,/? equal to the reference value 
iy) = average thickness of a panel having a value of 

t,/? equal to the reference value 
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The term ¢,/t,, represents the ratio of an arbitrary skin 
thickness to the skin thickness of the panel having ¢,/7 
equal to the reference value. The term //fy is the cor- 
responding ratio of the average thicknesses. 

The value of t,/? = 0.3 will be chosen as the reference 
value because it represents a limiting case for the curves 
of Fig. 1 and all values of /,/t,, and ?/% will be greater 
than unity. 

For each line of constant g/J,, the intersection with 
the line of ¢,/2 = 0.3 gives values of q/f) and q t,, for 
that particular value of g/L». Then, for each line of 
constant g/L, the values of g/f) and q/t,, are divided by 
the ordinates and abscissas defining that line. The re- 
sulting values are plotted in Fig. 4 as 7/f, vs. ¢, t. for 
lines of constant g/L». 

For example, Fig. 4 shows that, for a panel having an 
initial value of t,/7 = 0.3 and a structural index of 300 
Ibs. per sq.in., the skin thickness may be doubled at the 
cost of only 5 per cent increase in total panel weight if 
optimum design methods are employed. The dashed 
line shows the variation of panel weight based on 
the assumption that any increase in skin weight 
will increase the total panel weight by the same 
amount. 


CONCLUSIONS 


The curves of allowable stress F, vs. ¢/Lo represent- 
ing ‘‘optimum”’ design for given types of sheet-stiffener 
panels show that this minimum-weight design can be 
attained only at extremely low ratios of ¢,/?. Since 
wing panels are not designed on the basis of allowable 
compressive stress alone, curves of this type should 
show how much of the panel material is made up of 
skin and can be counted on to resist torsion. One 
method of indicating the relative amount of material 
in the skin is illustrated by the curves of Fig. 2, which 
are drawn for constant values of ¢, /?. 
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The allowable stress curves of Fig. 2 can be crog. 
plotted as shown in Fig. 3 and closely approximated py 
straight lines, resulting in a simplified analytical ¢. 
pression for the relationships involved. 
approximations of this type have been used conyey 
iently in wing-weight-estimation procedures.’ 

Fig. 4 illustrates the manner in which panel weight 
must change with changes in skin thickness and shoys 
the overconservatism of the assumption represented by 
the dashed line. For an actual design, where the values 
of both g and g/L are known, the curves of Fig. 1 coy) 
be used to calculate directly any change in pang 
weight for a given change in skin thickness. 

As the skin becomes an increasingly greater fractioy 
of the total panel thickness, the construction ap- 
proaches the ‘“‘panel’”’ type? in which a lighter construc. 
tion may be obtained by putting all of the mazterial in 
the skin and running the separators spanwise rather than 
chordwise. In this type of construction, which uses 
all the panel material to resist torsion, the compressive 
strength of the panel will depend on the width rather 
than on the length. If the data had covered a higher 
range of values of g/L) and skin thickness, the curves oj 
sheet-stiffener construction could be extended to show 
a comparison of these two types of construction and 
would then indicate the values of ts/ts, and q/Ly at which 
it would be advantageous to change from one type of 
construction to the other. 
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RIEF REPORTS of investigations in the aeronautical 

sciences and discussions of papers published in the JOURNAL 

are presented in this special department. The publication wi!l be 

completed approximately S weeks after receipt of the material. 

The Editorial Committee does not hold itself responsible for the 

opinions expressed by the correspondents. Contributions should 
not exceed 800 words in length. 


Correction for ‘‘The Effect of Variation of Mass 
on the Dynamic Stability of Jet-Propelled 
Missiles’’* 


M. V. Barton 
Professor of Aeronautical Engineering and Research Engineer, 
Defense Research Laboratory, Thz University of Texas 


August 1, 1950 


R. A. H. Fax, of the Cornell Aeronautical Laboratory, 

Inc., has pointed out that there is an error in the expres- 
sions for the jet damping force and moment of the air-jet-type 
missile. This error comes about by neglecting the sudden change 
of momentum of the incoming air at the entrance of the duct. 
The air at the entrance of the duct has a velocity component of 
Ql, at right angles to the undisturbed stream. The sudden 
change in momentum thus produces a force at the lip of c,Q/2 and 
a moment of c,Q/,?. When these quantities are added to the 
expressions for the forces and moments representing the integrated 
coriolis forces [Eqs. (10) of the paper] the jet damping force be- 


comes 
2[c(lo — hh) + cglolQ 
and the damping moment is 


[cqlo? + ¢(1p2 — 1,2) Q 


The net effect of this change can be effected by setting /» 
in the final equations. 

The general conclusions in the paper are not generally modi- 
fied by this change except that the possibility of negative damp- 
ing in the case that J» > Jp is precluded. 


* JOURNAL OF THE AERONAUTICAL SCIENCES, Vol. 17, No. 4, 197-203, 


April, 1950. 


PP 


Roll to Sideslip or Yaw Ratios 


Dunstan Graham and Allan McCaskill 
Flight Research, Cornell Aeronautical Laboratory, Inc., Buffalo, N.Y 


August 23, 1950 


R 


directional oscillation of airplanes as a possible pilots’ criterion 
for the upper limit of acceptable effective dihedral. Relative 
implitudes of roll and sideslip in unforced oscillations should not, 
however, be confused with the bank angle to sideslip ratio in cross 
controlled rectilinear flight which may be shown to depend en- 


EFERENCE | SUGGESTS the ratio of amplitudes of the rolling 
motion to the yawing or sideslipping motion in the lateral- 


tirely on other aerodynamic characteristics. 
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Since the solutions to the equations of lateral, stick-fixed mo- 
tion of an airplane are of the form? 


B, p, or r + Coed + Ck! (cos yt — C) (1) 


the ratios defining the relative magnitudes of rolling and yawing 
or rolling and sideslipping are functions of the three values of C;. 
These constants of integration are expressions containing the 
various airplane stability derivatives, mass factors, disturbing 
functions, and the roots of the characteristic equation raised to 
various powers and multiplied together in an involved manner 
Fortunately, however, the use of sweeping, simplifying assump- 
tions permits the development of approximate formulas for the 
roll-to-sideslip ratio and the roll-to-yaw ratio whose accuracy 
does not seem to have been unduly affected by the simplifica- 


= Ce 


tion 


Thus, if the lateral motions are assumed to be described by the 
emaciated differential equations 


d 
- +p - cS rin @ FE 
d(t/r) . 
j (2) 
Sse a wn, ri rd(t/r) = YON 
d(t/r) 


(these equations assume an undamped directional oscillation with 
the rolling motion coupled only through the dihedral effect and 
no lateral motion of the c.g.), the solution for the oscillatory com- 
ponents of p(t) and r(t) following a step yawing moment disturb 


ance gives 


PL) | os ¢ ¢ wl 3 
r(t)| y B n?V 1 + (1,/n)? 
where 

¢ = rolling oscillation amplitude 

vy = yawing oscillation amplitude 

8 = sideslipping oscillation amplitude 

m = airplane relative density parameter = m/pSb 

ly = nondimensional rolling acceleration due to sideslip (di 
hedral effect) = (1/2)(b kz)*Cig 

ly = nondimensional rolling acceleration due to rolling 
(damping in roll) = (1/4)(b/kz)?Ci, 

n, = nondimensional yawing acceleration due to sideslip 
(directional stability) = (1/2)(b kz)*Cng 

n = nondimensional frequency of oscillation (radians per 
unit ¢/r) = imaginary part of complex root of char 
acteristic equation 

k, = radius of gyration in roll (ft.) 

kz = radius of gyration in yaw (ft.) 

Ci, = rolling moment coefficient due to sideslip (per rad.) 

f 

Ci, = rolling moment coefficient per unit wing tip helix 
angle, pb/2U 

Cng = yawing moment coefficient due to sideslip (per rad 

T = unit of aerodynamic time = m/psU 

This result may also be derived by neglecting the effects of 


product of inertia and inclination of the flight path, assuming 
zero damping of the oscillation and the stability derivatives /, = 


Nr n, = Oin the complete expression for ¢/8. 
Similar assumptions applied to the complete expression for 
¢/y give as an approximation 
¢ by l 
y | my V1 + (p/n)? 


The same results are obtained no matter what the form of the 


(4) 


yawing or rolling moment disturbance may be. 


Under certain circumstances [as it is in Eq. (2)] 9 V un, is 
a valid approximation to the imaginary part of the complex root 
of the characteristic equation. The expressions for ¢/8 and for 


¢/y are then equivalent. At high altitude and high lift coeffi- 
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cients, however, n + V uns, and ¢/8 and ¢/y may be markedly 
different. 

The approximate expressions given by Eqs. (3) and (4) check 
numerically with the complete expressions (within 10-12 per cent) 
even for some extreme flight conditions and radical airplane 
configurations. Where the inclination of the airplane’s principal 
longitudinal axis is significant, the root » must be computed 
taking the changes in inertia and the product of inertia into 
account. Otherwise, the approximations are still good. For the 
combination of high dihedral effect and large positive inclination 
of the principal axis, the approximation fails completely, and the 
same is true for conditions of extremely low directional stability 
Some error exists in predictions (by the approximate formulas) 
of the airplane’s response to disturbances from its controls be- 
cause of the facts that the ailerons do not produce a pure rolling 
moment and the rudder does not produce a pure yawing moment. 
This is particularly true of some ailerons. In one typical case, 
however, other errors inherent in the approximation slightly more 
than canceled the error due to neglecting aileron adverse yaw, 
and the result was again in good agreement with the more exact 
calculation. 

Effects of design changes or trends on the roll-to-yaw and roll- 
to-sideslip ratios may perhaps be best discerned by considering 
the effects on the stability derivatives contained in the approxi- 
mate expressions. For instance, since /, is proportional to the 
dihedral effect and inversely proportional to the rolling moment 
of inertia, an increase in the dihedral effect will result in an almost 
exactly proportionate increase in g/8 and ¢/y. Decreasing the 
rolling inertia will tend to have the same effect, modified, however, 
by the corresponding change in /, 

If the roll-to-sideslip or -yaw ratios should achieve reputability 
as a Significant parameter, the approximate expressions may be 
useful for checking calculations of these quantities and as a 
guide to designers in the labyrinthine choice of fin area and di- 
hedral. 
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On the Plastic Buckling of Plates 


P. P. Bijlaard 
Cornell University, Ithaca, N.Y. 
August 8, 1950 


IT AN INTERESTING Readers’ Forum item (JOURNAL, June, 
1950, p. 378) Dr. Cicala calculates approximately the reduc- 
tion, caused by eccentricity in the buckling stress of a compressed 
plate with one unloaded side simply supported and the other free. 
Assuming the incremental (or flow) theory to be valid and assum- 
ing initial shear strains yo of 0.047,/G and 0.17,/G he finds ec- 
centric buckling stresses o, of about 0.75 and 0.66 times the 
concentric ones, respectively. 

Since these eccentric buckling stresses are for these eccentrici- 
ties and for this special case of the order of magnitude of the 
concentric buckling stresses according to the deformation theory, 
he suggests that in this way the excellent confirmation of the 
latter by buckling experiments may be explained. 

The assumed eccentricities, which obviously are suggested to 
be inevitable in concentrically loaded plates, are, however, rather 
large. From Dr. Cicala’s formula the initial angle of twist of the 
mid-section is 
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Bo max v% L/xs 
According to the simple incremental theory the concentric buck 
¢ 


ling stress is equal to the elastic one—i.e., using Dr Cicaly 


notations 

Ger, = G(s/b)? 
With o.,, = 1.5¢p, G = 4,000,000 Ibs. per sq.in. and a Proportion 
limit ¢, = 30,000 Ibs. per sq.in., Eq. (2) yields b/s = 9.5. gin, 
tp = 0.5770, and assuming L/b = 10 so that L/s = 95, 4, 


assumed initial strains yo of 0.04 r,)/G and 0.1 7)/G result froy 
Eq. (1) in initial angles of twist of 0.0052 and 0.0131, respectiyely 
maz. b of the fre 
O.75k, respectively, wher 


The corresponding maximum initial deflections ,; 
edge are 0.049s = 0.3k and 0.125s = 
k is the core radius k = h/6. 

These eccentricities are large, since the other unloaded edge is 
held straight by the support. Therefore the ratios u = o,/¢, < 
0.75 and 0.66 between eccentric and concentric buckling stress 
are not extraordinary low. For a simply supported steel] column 
of rectangular cross section and slenderness ratio 100, although 
having a smaller L/s ratio than the plate, initial eccentricities; 
is the 
respectively, 


of 0.25k and 0.75k cause ratios nh = o,/ocr, (where o,, 
tangent modulus buckling stress) of 0.70 and 0.35, 
if, somewhat analogous to Dr. Cicala’s approximation, the cag 
of initial end eccentricities e is considered. 

Thus, the influence of eccentricities on the column buckling 
stress is of the same order of magnitude as for the plate, |f 
eccentricities of this amount were really inevitable in experiments 
with theoretically concentrically loaded plates and columns, it 
should be expected that with concentrically loaded columns simi 
lar reductions in buckling stress would also occur. Nevertheless, 
the experimental buckling stresses for the latter case are in good 
the (tangent modulus 
The theoretical Shanley buckling stress is some- 


agreement with concentric Engesset 
buckling stress. 
what higher than the Engesser stress, since after deflection the 
resistance increases a little because of partial elastic unloading 
of the fibers. Hence, the inevitable eccentricities obviously caus 
only a decrease from the actual concentric buckling stress to the 
Conse 


quently, the inevitable eccentricities are much smaller than 


Engesser stress, which supposes fully plastic behavior. 


assumed by Dr. Cicala, and there is no reason to expect exper 
mental buckling stresses of plates which are substantially lower 
than the concentric ones, calculated from the present writer's 
theory’ in the same way as the Engesser stress—i.e., assuming 
fully plastic behavior. 

Dr. Cicala assumes that above the proportional limit, gp, the 
tangent modulus £; is constant and half of the elastic modulus £ 
although for aluminum a better approximation would be £; = 
E/30 or E/40. It could be thought that with lower £; the ratio 
u = a-/ocr, would be less. This may be easily checked by assum 
ing E, = 0, so that a, is the yield stress. 

Assuming ocr, = op, it follows that, from comparison of e 
ternal and internal torsional moments, with a compressive stress 

o(yo + ¥) = apy 
so that the shear stress 
t = Gy = Gyl[a/(op — @)]} 

Since, according to the yield condition 
o? + 3r? = ay? 4 


the yield stress o, is reached when, from Eq. (3), 


ao? + 3G°y?[07/(o, — o)?] = op? 
or, with 
yo? = a(7p2/G?) = (a?/3)(op?/G?) (8) 
when 
(1+ w)(1 — pw)? — arty? = O (6 
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in which 


a = a/Op 


After reaching the yield stress, however, with further twisting 
r will increase, SO that from Eq. (4) the compressive stress o will 
decrease. Hence, the maximum stress (since membrane stresses 
may be neglected here) and thus the eccentric buckling stress 
7 follow directly from Eq. (6), in which, consequently, 


M = O¢/Fp (7) 


With the assumed values yo = 0 047r,/G or 0.17,/G, or from Eq. 


5), with a = 0.04 or 0.1, Eq. (6) yields 


a = 0,/o, = 0.84 or 0.91 


It follows from Eq. (2) that here, with o., = op = 30,000 
ibs. per sq-in., b/s = 11.5, so that with L/b = 10, from Eqs. (1) 
and (5), Bo maz. = 0.158a. Hence, the maximum deflection 68 maz 
of the free edge is 1.83sa. With a = 0.04 or 0.1, this yields 
eccentricities of 0.44k or 1.1k. 

On the other hand, with the same stress-strain diagram and 
assuming also ocr. = op, for a rectangular simply supported 
column the ratio « = o¢/or. for eccentricities of 0.5k and k is 
053and 0.42, respectively,? and thus much less than for the plate 
Rven for an eccentricity of 0.1% it is still 0.75. 

Hence, in this case, the influence of eccentricities is even more 
for the column than for the plate, which confirms the above 
conclusion that there is no reason to expect experimental plate 
buckling stresses much lower than the concentric ones. 

The present writer’s theory*' includes the possibilities of plas 
tic behavior according to deformation (Hencky body), as well as 
flow theory (Prandtl-Reuss body), the latter case by equating 
value e = (E/Es) — 1, where Ey is the secant modulus, to zero. 
Assuming the actual behavior to be between that described by 
deformation and flow theories, he already gave his own explana- 
tion for the agreement of the experiments with his results ac- 
cording to the deformation theory (reference 4, pages 532-533), 
on which he based his calculations from the beginning in order to 
remain on the safe side.* Eccentricity plays only a secondary 
role here, as far as it prevents unloading of the material, so that 
the plate may be assumed to remain plastic all over (reference 
4, page 536). Moreover, it helps to increase o during buckling, 
so that in most planes the active shearing stresses increase, in 
which case the plastic behavior will approximate that according 


to the deform:ition theory 
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Plastic Column Behavior 


N. J. Hoff 
Head of Department of Aeronautical Engineering and Applied 
Mechanics, Polytechnic Institute of Brooklyn, Brooklyn, N.Y. 


August 16, 1950 


I‘ A THOROUGH ANALYSIS presented by Dr. John A. Duberg and 
Thomas W. Wilder, III, under the title ‘“Column Behavio1 
in the Plastic Stress Range” published in this JouRNAL in June, 
1950, the variations in stress and strain in columns buckling in 


FORUM 743 


the plastic range and the maximum loads the columns are able to 
carry after buckling were calculated. This writer does not want 
to add to, or detract anything from, the results of these calcula- 
tions, which are undoubtedly correct and contribute considerably 
to our knowledge of column behavior. However, he wants to 
mention a few basic issues related to our way of thinking about 
columns in general. 

First, it should be pointed out that the classical definition of 
stability is not ambiguous but leads to a unique buckling load 
even in the inelastic range. According to the classical procedure, 
the buckling load is calculated as the load under which the 
stability of the system is neutral. The system consists of a given 
elastic body upon which given loads are acting. 
constant even though their magnitude is unknown. 
seen from Fig. 2 of the paper quoted that infinitesimal deflections 
are possible without a change in load only at the reduced modulus 
This is the unique buckling 


The loads are 
It can be 


load Pru derived by von Karman. 
load according to the classical theory. 

The load of real interest to the engineer is that obtainable dur 
ing some loading process rather than the theoretical neutral 
stability load. This is so because columns will seldom collapse 
under constant loads but rather during the period when they are 
loaded. To talk of the stability of a system that is in a state of 
continuous, or even discontinuous, change of loading (the latter 
being the case when dead weights are used in loading the elastic 
body) is not possible without defining first what is meant by the 
term. One possible definition of the buckling load is the maxi- 
mum load the column can support during the loading process. 
This definition has the advantage that it can be easily correlated 
to buckling tests in which the maximum load is the principal 
quantity measured. The intrinsic importance of the maximum 
load was recognized by the authors, inasmuch as they presented 
several diagrams containing it. Another proper choice of the 
buckling load is the load at which the deflections of an initially 
slightly curved column become inadmissibly large. The value of 
this load can be taken from diagrams of the type shown in Fig. 4 
of the subject paper after the meaning of the word ‘‘inadmissible”’ 
is defined in some exact manner. Finally, the authors’ definition 
is also a feasible one: The buckling load can be taken as the load 
at which an initially straight column begins to deflect if its be- 
havior is assumed to be the limit of the behavior of an initially 
slightly curved column whose initial deviations from straightness 
are made to approach zero. This definition is only one of many 
It has no particular theoreti- 


or even likely—that a 


other equally arbitrary definitions. 
cal significance, and it is not certain 
perfect column would behave as indicated by the limiting proc- 
ess. Its practical importance is also not obvious, since the 
slope of the load-deflection curve of an initially slightly curved 
column changes abruptly below the tangent modulus load if Fig 
$ is indicative of the behavior of real columns and since the 
maximum load the column can support is always greater than the 
tangent modulus load as shown in Fig. 9. 

Moreover, it must not be forgotten that the entire analysis of 
the paper is based on the assumption of an infinitely slow mono- 
tonically increasing loading process in which the inertia forces 
That this last assumption is permissible is not 
For this reason, under 


can be neglected. 
obvious to the writer of this discussion. 
the sponsorship of the Office of Naval Research of the U.S. Navy, 
he undertook the investigation of the actual motion of a column 
during the process of loading in a testing machine making use of 
the dynamic equation of motion. The results of the calculations 
as related to perfectly elastic columns were presented under the 
title ‘The Dynamics of the Buckling of Elastic Columns’”’ at the 
Annual Summer Meeting of the Applied Mechanics Division of the 
A.S.M.E. in June, 1950. Work is now in progress at the Poly- 
technic Institute of Brooklyn on the problem of the inelastic 
column. 

Finally, it may be mentioned that the definition of instability 


admits of a further generalization. A system should be con- 
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sidered stable if its equilibrium is automatically re-established 
after a disturbance. In systems that are not perfectly linear, the 
magnitude of the disturbance is a decisive factor, and in the pres- 
ence of dry friction o1 plastic deformations the original state of 
equilibrium is regained, as a rule, only approximately. For these 
reasons the writer proposed in his paper “Dynamic Criteria of 
Buckling” (Research, Engineering Structures Supplement, Butter- 
worths Scientific Publications, Ltd., London, and Academic 
Press, Inc., New York, 1949, page 121) that the stability of such 
systems be investigated by determining their motion following 
a disturbance. They should be considered stable if after a reason- 
able time they return reasonably close to their original position 
of equilibrium. The nature and the magnitude of the disturbance, 
the specification of a ‘‘reasonable”’ time, and the definition of the 
expression “‘reasonably close’ depend upon the circumstances 
and cannot be established in advance for all possible systems. 
The paper discusses the problems involved with the aid of simple 
dynamic models. 

The writer wishes to emphasize once more that it is not his 
intention to minimize the substantial advances made by Dr. 
Duberg and Mr. Wilder in column theory. He feels, however, 
that further improvements of our knowledge of column action are 
possible, and he wants to indicate the direction in which, in his 


opinion, they might be attained. 


A Note on the Calculation of Oblique Shock- 
Wave Characteristics* 


M. J. Thompson 
Professor and Chairman of the Department of Aeronautical 
Engineering; Associate Director, Defense Research Laboratory: 


The University of Texas 
July 27,1950 


5 ps NOTE IS CONCERNED WITH a simplified method for han- 
dling the numerical calculations involved in determining the 
characteristics of an oblique shock wave in supersonic flow. 
The wave is considered as being formed at the leading edge of a 
two-dimensional airfoil section having a nose angle 0, between 


its surface and the free-stream direction. The problem origi- 
nated as a part of a study of airfoil pressure distribution.! 

Using the method of analysis employed by Taylor and Mac- 
coll,? the wave angle, pressure ratio, downstream Mach Num- 
ber, and other desired flow characteristics may be determined in 
terms of the wedge angle and the free-stream Mach Number. 
The method of solution, however, is an indirect one and generally 
requires the use of elaborate chaits and a considerable amount 
of cross-plotting. By using this method it is possible to set up 
an expression for the pressure ratio in the form of a quartic, 
with algebraic coefficients that depend on the free-stream Mach 
Number 1M, the wedge angle @,, and the ratio of specific heats, 
y, for the gas. 

The present analysis follows the method used by Liepmann 
and Puckett? and others. The relation from which the wave 
angle @,, may be determined is 

* Acknowledgment is made to the Navy Bureau of Ordnance for its support 
of this work under Contract NOrd-9195, which is a part of the Bumblebee 
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+ 1\ sin 0, sin @, 

2 cos (0. — 65) (1) 
Expressions for the other characteristics may also be derived ig 
terms of the known quantities defined above and the wave angle 


computed from Eq. (1). In order to use Eq. (1), Laitones for 
example, assumes a series of values of 0, for a given M, and com 


EMBER, 


putes @,. 


Eq. (1) may be rewritten in the form 
sin* 0, + 6 sin‘ 0, + c sin? 6, + d = 0 


where 


b = — t[(vy + 1)M? + 2] tan? 6, + M,? + 2} 
{1 + tan?@, 
(vy + 1I)M? + 2)? tan?@, + 8,2 + 4} /4M,! (1 + tantg) 


(4) 
d = —1/M,' (1 + tan? 6,) (5) 
which is a cubic in sin? @, and may be solved directly by numeri. 
cal procedures such as Tartaglia’s method. Thus, any desired 
degree of accuracy can be maintained in the calculations, and, 
what is more important, 0. is obtained directly for any given 
values of M,, 0,, and y without the necessity of trial-and-error 
solutions and cross-plotting. Such a procedure is more readily 
adapted to comparison with wind-tunnel test results in which 
M, and 6, are the primary independent variables. 
It is also interesting to note that, in the special case where 
. = 0, Eqs. (3), (4), and (5) reduce to 
—(M,? + 2) 2M,2 + 1 1 
bo = : , Oo , &= (6) 
M;? M,! M;! 
Furthermore, Eq. (2) is now factorable, the significant roots 
being 


sin 0, = 1/M,, sind, = 1 (7) 


The second root, of course, has no physical significance, while 
the first corresponds to the Mach wave created by an infinitesi- 
mal point disturbance. 
The values of Eq. (6) may be introduced into Eqs. (3), (4), and 
(5) so that the latter become 
b = —[(M2 + 2)/M,?] 


2M? + 1 (y¥+D?. (y¥- 
M\* +| 4° Me 


—cos? 0,/M;,‘ (10) 


— vy sin? 6, 


While these forms do not permit the factoring of the original 
expression, Eq. (2), nor do they provide for a complete separation 
of the effects of Mach Number and wedge angle, they are be 
lieved to be somewhat simpler in form and better adapted to 
numerical calculation procedures than Eqs. (3), (4), and (5). 
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